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Dedicated to the centenary of Mark Krein 



Preface 

In the recent years, the theory of orthogonal polynomials on the real line (OPRL) 
and on the unit circle (OPUC) enjoyed the considerable development. In these lec- 
ture notes, we will explain how to construct the continuous analogs of polynomials 
orthogonal on the unit circle. It is possible to built a theory which is as rigorous and 
complete as the theory for OPUC. Spectral theories of one-dimensional Dirac and 
Schrodinger operators can then be viewed in the framework of this theory which 
establishes a solid link between an approximation theory and quantum mechanics. 

The theory is based on the ideas suggested by M.G. Krein. They were developed 
later by various authors, especially from Krein's school. In the meantime, new 
results were obtained for OPUC and OPRL and that was a motivation for us to try 
to understand their continuous analogs. We also try to give systematic exposition 
of the theory but have to refer to the literature once in a while. Also, these notes 
do not cover some aspects of the theory (e.g. continuation problems for G r classes, 
regularity of coefficients, etc.) but we give necessary references. In general, our 
objective is to give only basics of the theory by presenting complete proofs and 
filling various gaps present in the current literature. As a prerequisite for reading 
these notes, we assume that a reader is familiar with main facts from the OPUC 
theory (see, e.g. [7H [551 1ES])- The knowledge of spectral theory for Schrodinger 
and Dirac operators might also be very helpful. 

What is not covered and what is new? 

We didn't include the following subjects that are related to our topic: solution 
to the continuous analogs of Schur and Caratheodory-Toeplitz problems [45l 146] . 
We also do not discuss matrix- valued version of the theory. For the recent progress 
on more specific questions (such as continuous analog of Szego case, Rakhmanov's 
Theorem, etc.) we suggest the reader to consult the journal publications, e.g. 
[ZH [711 [TJl [T5l [16] . Also, we will deal with rather regular classes of coefficients (not 
worse than Lj 2 oc (R + )) but the general case can also be treated in the framework of 
different differential operators (see, e.g., [52 | l53 l 16]). 

In these notes, we present quite a few new results. That includes: approximation 
of continuous orthogonal system by the sequence of the discrete ones (Section 8), 
distribution of zeroes (Section 9), new criteria for A(r) € i 2 (K + ) and more on 
that case (Section 10), the continuous analog of the Strong Szego Theorem - sharp 
conditions (Section 14). We also gave complete proofs for results that were present 
in the literature without any proofs and gave alternative (hopefully, more trans- 
parent) proofs for several other statements (e.g. the continuous analog of Baxter's 
Theorem, scattering theory for Krein systems and Dirac operators). 
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Notations 

D open unit disc in C 
T - unit circle in C 
H p (il) - Hardy space in the domain fi, < p < oo 
JV(f2) - Nevanlinna class of analytic functions in fl 
B(Q) - closed unit ball in H°°(n) 
IP A - the following projection: 



Ta 



exjp(iXx) f (x)dx 



J exp(iXx)f(x)dx 



TJ2 

[0,R] 



RanTjcij] 
denote CPr± respectively 
absolute value of operator O 
the delta- function at zero 
Wiener's Banach algebra of functions 

oo 

/(A) = / f{x)exv{iXx)dx,f{x) G L l ( 



VcFo 

5(x) 
W(M) 



W + (C + ) - Banach algebra of functions 

oo 

/(A) = J f(x) C xp(iXx)dx,.f(x) e L 1 (R+) 



J ff 1 /2( R ) _ fractional Sobolev space of functions / whose Fourier transform satisfies 

oo 

J <t> \f(t)\ 2 dt < oo, < t >= Vt 2 + 1 



Co(fi) 
Xa{x) 

n A 

§p 

ln + x 
In - x 
f*9 



continuous on CI functions vanishing at the boundary of f2 
the characteristic function of the set A 
the orthogonal projection in L 2 (M + ) onto L 2 (A), A C M 
Schatten-Von Neumann class of compact operators 
= In x if x > 1 and = for < x < 1 
= In x if < x < 1 and = for x > 1 
means the convolution of / and g 

We usually use calligraphic letters to distinguish between operators and func- 
tions. For example, stands for an operator and 0(x, y) denotes the function of 
two variables. 
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1. Some classes of functions on the real line 

In this section, we recall some basic facts on positive definite functions on the 
real line. Then, we introduce certain class of functions that we will use later on. 
Let < r < oo. 

Definition 1.1. The Lebesgue- measurable function <j>(x) defined on the interval 
(— r, r) is called Hermitian if <j){—x) = <p(x) for a.e. x. 

Definition 1.2. The Lebesgue-measurable function of two variables K(x,y) is 
called Hermitian if 



K(x,y) = K(y,x), (1) 

for a.e. < x, y < r. 

Definition 1.3. The integral kernel K(x,y) is called positive definite on [0,r] if 
for any N, {%j}f = i, {%j G [0, r]), {cj}^ =1 , (cj € C), we have inequality 

N 

^ c n c m K(x n ,x m ) > 0. (2) 

n,m=l 

Consider the integral operator % in L 2 [0, r] with kernel K(x, y) £ C([0, r] 2 ), i.e. 

r 

(Xf)(x)= [ K(x,y)f(y)dy. 



Clearly, (jTJ) means %* = %. It is an easy exercise to see that if the continuous 
kernel is Hermitian, then ^ is equivalent to 3C > 0, where inequality is understood 
in the operator sense. 

Definition 1.4. A function <j>(x) is called positive definite if the integral kernel 
4>(x — y) is positive definite on M + . 

This is equivalent to <p(x — y) being positive on the whole line K. 

Notation 1. The class of continuous positive definite functions on the whole line 
is denoted by Poo- 

If dfi is finite positive measure on R, then 

exp(ixt)dfj,(t) (3) 



is positive definite. The classical result of Bochner says that the converse state- 
ment is also true. That, in a sense, is the continuous analog of the solution to 
trigonometric moment problem. 

Theorem 1.1. (Bochner, [3]) A function <f> belongs to the class Poo if and only if 
it admits the representation (0) with finite positive measure [i. The measure [i in 
this representation is unique. 

Notice that Bochner's theorem implies that all Poo functions are necessarily 
bounded on HL 
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Notation 2. Let Goo denote the class of continuous Hermitian functions g(x) 
defined on the whole line such that g(0) = and the integral kernel 

K(x, y) = g(x) + g(-y) - g(x - y) 

is positive definite on M + , i.e. on any interval [0,r],r > 0. 

Next, we obtain some rather crude estimates on g S Goo. Later, these bounds 
will be used to prove the integral representation for functions of class Goo- 
The following inequality holds 

\g(2x)\ < 8\g(x)\ (4) 

for any x. Indeed, since the kernel K(x,y) is positive definite, estimate ([2]) is true. 
Take N = 2, x\ = x, xi = 2x. If c\ = £ € R, C2 = 1, we have 

Re(g{x))£ 2 + Re(g(2x))£ + Re(g(2x)) > 

Since £ is arbitrary real, 

\Re(g(2x))\ <AReg(x) (5) 
For the same choice of £1(2), we let c\ = i£ 6 iM., ci = 1. Then, 

Re(g(x))e - [2Im(g(x)) - hn(g(2x))]t + Re(g(2x)) > 

That yields 

[lm(g(2x)) - 2Im( 3 (x))] 2 < 4Re g(x)Re(g(2x)) 
Using (JSJ), we have 

|Im(.g(2x))| < 2\1mg(x)\ +AReg(x) 

Combining estimates for the real and imaginary parts, we obtain (j4|). 
The following estimate holds true 

\g(x)\<C(l + \x\ 3 ) (6) 

Indeed, if max l6 [_ L1 ] \g{x)\ — M, then max[_ 2 ^ 2 »] \g( x )\ < 8™M by ([!]). Therefore, 
\g{x)\ < 8l log 2 l x l] +1 Af < G|x| 3 , where [.] means the integer part of the a number. 
As we will see later, the estimate (|6]) is very far from optimal. 

The following integral representation of Goo functions is an analog of Bochner's 
theorem for class Goo- 

Theorem 1.2. Function g{x) S Goo if and only if 

00 

g(x) =i(3x+ J (l + - exp(i\x)^ (7) 

— OO 

where j3 G M and positive measure a satisfies the estimate 

00 

da(X) /cA 

T+a 2 < °°- < 8 > 

— oo 

Constant f3 and measure o are uniquely defined. 
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Proof. Any functions of the form ([7|) belongs to Goo ■ Indeed, notice that the integral 
in (O converges if §§§ holds and defines the continuous function that vanishes at 
zero. Then, we have the following representation 



0O) + g(-y) - a( x - v) 



(1 — cxp(iAx))(l — exp(iAy)) 
A 2 



da(\) 



which ensures the positivity of the operator with the corresponding kernel for any 
r > 0. 

Conversely, due to ([6]), any Goo function allows Laplace transform. Consider 

00 

L(z) = z 2 J g{x) exp(ixz)dx 


This function is analytic in C + . Notice that 

00 00 

M 2 / / [9(x -y)~ SO) - g(-y)} exp(ixz - iyz)dxdy < 



L(z) + L{z) , 



21m z 







Consequently, —iL(z) is Herglotz function and has well-known integral representa- 
tion ([3], chapter 6) which gives 



L(z) — iaz — i/3 — i 



1 - Xz 



(A + z)(l + A 2 ) 



da(X) 



where 



da(\) 



< 00 



1 + A 2 

and a > 0,(3 g R. Let us take the inverse Laplace transform. Notice that 



A 2 (l- Xz) 
2 (A + z)(l + A 2 ) 

Therefore, 



1 



iXx 
1 + A 2 



— exp(iAa;) 



exp(ixz)dx, z e C + , A g E (9) 



= iaz — i(3 + z 2 / exp(iccz) 



iAx 
1 + A 2 



exp(zAa;) 



dg(A) 
A 2 



dx 



1 = — z 2 / xexp(ia;z)(ia;, z = —iz 2 / exp(ia;z)(ix, z g 



Since 



we have the formula 

g(x) = a + i/3a; + 



(10) 



iAa 



1 + A 2 



exp(iAa;) 



dajX) 



Due to normalization g(0) = 0, a = 0. Uniqueness follows from the uniqueness of 
the Herglotz function representation. □ 
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As a simple corollary one gets the following improvement of (J6j) : \g(x)\ < C(l + 
x 2 ). If cr in (JTJ) is Heaviside function and = 0, then g(x) — x 2 /2. So, the quadratic 
growth is possible. 

If the support of a is a compact, then the second derivative of g exists and 
is positive definite. That follows from the Bochner's theorem. In general, the 
derivatives of g{x) at zero can have singularities. In what follows, the reference 
measure is cr(A) = \/(2tt), and /3 = 0. It is then easy to check that g(x) — \x\/2. 
So, in this case, the second derivative in the distributional sense is delta function. 

The relation between Goo and Poo can be established by 

Lemma 1.1. If f{x) € Poo, then /(0) - f(x) E Goo- 
Proof. From Bochner's theorem, we have 



f(x) = I exp(ixt)dfi(t) 



A OO 

Take cr(A) = J t 2 d^{t), f3 = - J t(l + t 2 )' 1 d(J,(t) and use Theorem 2. □ 

-oo 

The converse is wrong. Take f(x) — ix. The corresponding kernel K(x, y) = 
0. At the same time, f(x) ^ Poo because f(x) is not bounded. The class Goo 
is convenient for description of measures generating one very important class of 
canonical differential systems, the Krein systems, which we plan to study in the 
next sections. 

Remarks and Historical Notes. 

Classes P r , G r can be introduced for finite r > (see [3], p. 190 and references 
there). Then, measures /i and a are not uniquely defined in general. The class of 
nonuniqueness for a corresponds to the continuation problems [46j . The case when 
quadratic form ^ has not more than >c negative squares is more difficult. It was 
studied in the framework of Pontryagin fl^-spaces 46; . 
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2. Factorization of integral operators 

To understand better the algebraic aspects of the theory, we will need some 
rather simple results on the factorization of integral operators. As we know from 
the linear algebra, given any matrix A — {o,ij}f t j—i with nonzero leading principal 
minors, we can always find a lower-triangular matrix X\ and an upper-triangular 
matrix X2, such that A — X1DX2, D is a diagonal matrix. The proof is simple. 
Since an ^ 0, using the first step of Gauss algorithm, we can find the lower- 
triangular matrix L\ of elementary transforms, such that L\A has the first column 
collinear to [1, 0, . . . , 0]*. Then, find an upper-triangular U\ such that the matrix 
L\AU\ has the first raw collinear to [1, 0, . . . , 0]. Notice that the first column stays 
the same. The leading principle minors of L\AXJ\ are the same as those of A. Thus, 
we can continue this process. In the end, we get a lower-triangular L = Ld ■ ■ ■ L\ 
and an upper-triangular U = U\ . . .Ud such that LAU = D. Denoting X\ = L , 
X2 = we get the desired statement. Notice that matrices Xi and X2 have 1 

on the diagonal. Taking inverse of A, we get the factorization in the reverse order. 

Now, what happens in the continuous case? First, we need to establish one 
general result about the resolvent kernels of integral operators. Fix some R > 0. 
For any < r < R, consider integral operators 

r 

Xrf(x) = J K(x,y)f(y)dy, 
b 

with a kernel K(x 1 y), continuous on [0, R] 2 , and acting in the Hilbert space L 2 [0, r]. 
Assume —1 ^ Spec(X r ) for any < r < R. Then, the resolvent kernel T r (s, t) exists 
and 



r r (s, t) + J K(s, u)T r (u, t)du = K(s, t), < s, t < r (11) 


Lemma 2.1. Let K(x,y) £ C{[0,R} 2 ) and -1 ^ Spec(X r ) for any < r < R. 
Then. 

dT r (s,t)/dr = -T r (s,r)T r (r,t), 0<s,t<r 

Function T r (s,t) is jointly continuous in s,t, r and is continuously differentiable 
in r. 

Proof. Notice that Fredholm's formula for resolvent kernel 

r r (M) = ^M ( i2) 

O r 

ensures that T r (s, t) is jointly continuous in s, t, r and is continuously differentiable 
in r, r > (Appendix, Lemma 117. lj) . 
Then, differentiate ([IT]) in r. We have 

r 

+ J K(s,u)^^du = -K(s,r)T r (r,t) (13) 
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On the other hand, multiplying both sides of 

r 

T r (s, r) + J K(s, u)T r (u, r)du = K{s, r) (14) 
o 

by — r r (r, f), we obtain an equation 

T" 

-T r (s,r)T r (r,t) - K(s,u)T r (u,r)T r (r,t)du = -K(s,r)T r (r,t) (15) 



Now we see that both dT r (s,t)/dr and — T r (s, r)T r (r, t) satisfy the same integral 
equation. Therefore, they are equal. □ 

Although the continuous kernels are very natural, we will also work with K(x, y) 
from the following class. 

Definition 2.1. Function K(x, y) belongs to the class C([0, R} 2 ), if it is continuous 
in each of the triangles: A + = {0 < x < y < R} and A_ = {0 < y < x < R} 
but might have discontinuity on the diagonal x = y if considered as the function 
on [0, R] 2 (i.e. the limits K+(x,x) and K-(x,x) might be different). 

For this case, we have an analogous statement 

Lemma 2.2. Let K{x,y) £ C([0,R} 2 ) and -1 £ Spec(X r ) for any < r < R. 

Then, 

dT r (s, t)/dr = -r r (s, r)T r (r, t), s,teA ± (16) 

Function T r (s,t) G C([0,r] 2 ). It is continuously differentiable in r for s,t € A± 
and the derivative dT r (s,t)/dr is also from C([0, r] 2 ). 



Proof. The proof is similar to the proof of Lemma 12.11 Instead of the usual 
Fredholm formula for resolvent kernel we need to use a modified one given in 
Lemma 117.21 Indeed, analysis of d r shows that it is continuously differentiable. 
The Carleman-Hilbert determinant S r (s,t) has derivative in r for s,t fixed in each 
of A±. This derivative is also from the class C*([0, r} 2 ). □ 

Notice that dT r (s,t)/dr can be regarded as a continuous function on the whole 
[0,r] 2 because the right-hand side of (fll)]) is continuous on [0,r] 2 . 

The natural analog of the lower-triangular matrix is Volterra integral operator 

X 

Lf(x) = J L{x lV )f(y)dy 
o 

acting on the Hilbert space L 2 [0,R]. We also assume that L(x,y) £ C(A_). An 
operator U is upper-triangular, if li* is lower-triangular. The product and the 
sum of two lower(upper)-triangular operators are lower (upper)-triangular as well. 
Operators I + L and I + U are both invertible and (I + — I is lower-triangular, 
(7+U) -1 — I is upper-triangular. Infact, there is the Banach algebra of lower(upper)- 
triangular operators [ 321 [30J . Assuming that we have factorization 

I + X R = (I + L){I + U) (17) 

where H is lower-triangular and U is upper-triangular, we immediately get that 
I + % R is invertible and (I + Xr)- 1 = I - Sr = (I + V+)(J + V_), where V± 



10 



SERGEY A. DENISOV 



is upper(lower)-triangular with kernels V±(x,y),0 < x,y < R and we have the 
following formula for the resolvent kernel 



R 

-V+(x,y) - J ' V + (x,u)V-(u,y)du, x<y 

( 18 ) 

-V-(x,y)- / V+(x,u)V-(u,y)du, x>y 



It should also be mentioned that 3C r , (0 < r < R) can be factorized just by using 
truncations of £ and U. 

Theorem 2.1. The integral operator %r with kernel K(x,y) G C([0, R] 2 ) admits 
factorization |J7| ) if and only if I + % r is invertible in L 2 [0, r] for any < r < R. 
In this case, 

V+(x,y) = -T y {x,y),x<y , ig . 
V-(x,y) = -T x (x,y),x>y 
where T r (x,y) denotes the resolvent kernel of I + X r . 

Proof. Indeed, assume that I + 3C r is invertible for any < r < R and T r (x,y) 
is the resolvent kernel. Define V± by (fl9|) . Now, let us check (fT8|) . Indeed, from 
Lemma 12.21 

R 



T y (x,y)- / T u (x,u)T u (u,y)du = 



—T u (x, y)du = T R (x, y),0 < x < y < R 



The case x > y can be checked in the same way. Thus, we have (|18[) . which means 
that I— Qr = (I + V+)(/ + V_). Now, (fT7|) is straightforward. Conversely, assume 
that the factorization (fT7|) exists. Then 

i + x r = n [0>r] (i + £)(! + u)n [0 , r] = n [0 , r] (i + £)n M n M (i + u)n [0 , r] 

and it is clearly invertible. □ 

Another important class of factorizations is the following one. Instead of integral 
operator on L 2 [0, R] we consider an integral operator on L 2 [—R, R] and define the 
lower-triangular operator as 

M 

tf(x)= [ L(x,y)f(y)dy 



-\x\ 

where L(x,y) is continuous in {0 < x < R, \y\ < x} and in {— R < x < 0, \y\ < \x\}. 
Introduce Cl^ = {\y\ < \x\ < R}, tl+ = {\x\ < \y\ < R} and redefine A_ = {-R < 
y < x < R}, A+ = {-R <x<y<R}. 

Similarly, we say that U is upper-triangular if U* is lower-triangular. These newly 
defined lower-triangular operators possess the same algebraic properties: sum and 
product of two lower-triangular operators is lower-triangular, (/ + XL) -1 — / exists 
and is lower-triangular. The same is true about the upper-triangular operators. In 
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general, the definition of lower (upper)-triangular operator depends on the choice of 
the so-called chain of orthoprojectors ([30], Chapter 4). 
The natural question is when can we factor the operator 



where X r is defined as 



We have 



I + Xr = (/ + £)(/ + U) 



X r f = / K(x,y)f(y)dy,0<r <R 



(20) 



Lemma 2.3. Let K(x,y) G C([-R, R} 2 ) and -1 £ Spec(X r ) for any < r < R. 
Then, 



dT r (s, i)/dr = - f r (s, r)r r (r, t) + T r (s, -r)r r (-r, t) 



, s,t £ A± 



The function T r (s,t) G C([— r, r] 2 ). is continuously differ entiable in r for s, t G 
A± and f/ie derivative dT r (s,t) / dr is also from C([—r,r] 2 ). 

Proof. We will only check the formula for derivative. The other properties can be 
checked just like in Lemma 12.21 We have 



(21) 



Differentiating in r, we get (s ^ t) 

dr r (s,t) + r ^ (aiU) ^M) du = _ ( ^ (fljr) f. r(r)t) + ^ (aj _ r) f r( _ r>t)) (22) 



r r (s,t)+ / K(s,u)t r (u,t)du = K(s,t) 



dr 



dr 



On the other hand, take ([21]) with t = ±r, multiply by —T r (±r,t) and add: 



-(r r (s, r)r r (r, t) + f r (s, -r)r r (-r, t) 



+ 



+ J K(s,u)- [-(f r (u,r)f r (r,t)+f r (u,-r)f r (-r,t) 

— r 

-(A-(s, r)f r (r, t) + -r)f r (-r, t)) 
Comparison with l|22p finishes the proof. 



du = 



□ 



Now, if we have ([20]), then (L + Xr)- 1 =I-5r = (I + V+)(I + V_), where V± 
is upper(lower)-triangular with kernels V±(x, y), —R < x,y < R. Moreover, for the 
resolvent kernel 



-V+(x,y)- J V + (x,u)V-(u,y)du, (x,y)eO + 
|»|<|«|<fl 

-V-(x,y)- J V + {x,u)V^(u,y)du, (x,y)eQ,_ 

\x\<\u\<R 

The next Theorem provides the needed factorization 



(23) 
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Theorem 2.2. The integral operator 3Cr with kernel K(x, y) G C([—R, R] 2 ) admits 
factorization h2U\) if and only if I + X r is invertible in L 2 [—r,r] for any < r < R. 
In this case, 

V+{x,y) = -f\ y \(x,y),(x,y)en + 
V-(x,y) = -T\ x \(x,y), (x,y) e Q_ 



(24) 



Proof. Calculating the right-hand side in (|2"3j) and using Lemma 12.31 we get (for 
(x,y) € &-): 

-V-(x,y)- / V + (x,u)V-(u, y)du = 



|x|<|«|<fl 
R -\x\ 



T\x\(x,y) - / T u (x,u)T u [u,y)du - I T\ u \(x,u)T\ u \(u,y)du 



\x\ 
R 



= T\x\(x,y) - / T u (x,u)T u (u, y)du + T u (x, -u)T u (-u, y) 



du 



= F R (x,y) 

The case (x, y) G fl+ can be checked similarly. So, (|23[) is true. The other state- 
ments of the Theorem can be verified following the proof of Theorem 12.11 □ 

Notice that V^(x, — x) — V+(x, — x), x ^ 0. The results of last Lemma and a 
Theorem can be easily generalized to the case when the kernel K(x, y) is allowed 
to have a discontinuity of the first kind on {(x,y) : y — —x}. We do not do that 
since the class C([—R, R] 2 ) is exactly the one we will need later on. 

Remarks and Historical Notes. 

The proofs of the results in this section are partially taken from [3D] . In [3D] , the 
general case of factorization along the chain is considered. Recently, the factoriza- 
tion problem for integral operators with less regular kernels was studied in [SH [5S] • 
Later on, we will need to use the factorization of Fredholm integral operators along 
with regularity properties of the kernels. 
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3. Continuous analogs of polynomials orthogonal on the unit circle 

In this section we start building the theory of continuous analogs of polynomials 
orthogonal on the unit circle (OPUC). For the OPUC basics, we refer the reader to 
[72 [26l EH [40] . Let H(x) be Hermitian function defined on K and H(x) £ L^O^) 
for any r > 0. In the Hilbert space L 2 [0, r], consider the following integral operator 

r 

3irf(x) = J H(x ~ t)f(t)dt, 0<x<r (25) 



This operator is called "truncated Toeplitz" operator or operator with the "dis- 
placement kernel" [551 US]- It is obvious that for any r > 0, this operator is 
self-adjoint, compact, and its lower (upper) bound decreases (increases) in r. 

Definition 3.1. Function g{x) £ Goo has an accelerant if there exists Hermitian 
function H(x) (accelerant) defined on R such that 

X 

9{x) = ^+J(x- S )H(s)ds (26) 
o 

for all 

Theorem 3.1. The function H(x) generates g(x) £ by formula i26\) if and 
only if 

I + "K r >0 

for any r > and inequality is understood in the operator sense. 

Proof. It is obvious that g(x) is Hermitian, continuous, and g(0) — 0. Consider 
any ip(x) £ C°°[0,oo) with compact support on [0, r]. For the kernel K(x,y) = 

g( x ) + g(-y) - g( x - y)> we have 



K(x, y)if' \y)ip' (x)dxdy =h+I 2 - h 



o o 



For 1 1 and I2, 



h=h = —<p(0) / g(x)ip'(x)dx 







if'(x) J g(x-y)ip'(y)dydx = -<p(Q) J (p / (x)g(x)dx+ / ip'(x) / g' (x-y)<p{y)dydx 


Using |26)) and integrating by parts, we have 



K(x,y)ip'(y)ip'(x)dxdy = J J (p(y)(p(x)H(x-y)dxdy+ J \tp(x)\ dx = ((1 + 5C r )<p, y) 
00 00 

Thus, if g S Goo, then I\ + 12 — I3 > and 1 + Ji r > for any r > 0. Conversely, 
assume that 1 + "K r > for any r > 0. Take any ip(x) £ C°°[0,oo) with compact 
support in [0,r]. It can be written as 
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where 

oo 

<p{x) = — J il>(s)ds 

X 

and <p(x) £ C°°[0,oo), ip(x) is supported on [0, r}. Therefore, 



K(x, y)ip(y)ip(x)dxdy = ((1 + !K r )(p, ip) > 



o o 



and g(x) £ G^. □ 

As a corollary from Theorem 11.21 and Theorem [3HJ we get the following formula 
for an accelerant 

X oo 

M + J(x - s)H(s)ds = l px+ J (l + ^ exp(zAx)) ^ (27) 

-co 

where (3 £ R, and 

CO 

dajX) 
1 + A 2 



< oo 



The straightforward calculation shows that the trivial case H (x) = corresponds 
to = and a (X) = A/(2tt). 
Essentially, l|27p means that 

exp(i\x)da{\) " = " 5(x) + H{x) 



exp(iXx)d(a(X) — o"o(A)) " = " H{x) 

In other words, H(x) are "moments" of a — oq. Clearly, the integrals in the last 
formulas do not have to converge in the usual sense. 

Lemma 3.1. Assume that a from £?[) is known and g has an accelerant, then the 
constant (3 is defined uniquely by the formula 

13 = -itf'(O) (28) 

where 



CC: 



*(*) = " / (l + jf^ - «P(«*)) d[a(X) - a ° W] £ C\R) (29) 

— CO 

Proof. The proof follows from the formula (|2T|) by taking the derivative. □ 

It is up to us to choose the regularity class for H(x). In these notes, we will 
consider two important cases: 

H(x) £ C[0,oo), (30) 

and 

H(x) £ if oc (R) (31) 
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Other classes of regularity (e.g., H{x) G if oc (R),p > 1) can also be treated. Al- 
though, the case p = 1 needs special consideration (see discussion in [46j 145])- 

We will start our construction with the continuous accelerants. Then the L 2 oc 
case will be treated by an approximation argument in the separate section. For 
(|30| . H(x) might have discontinuity at but the left and the right limits must 



exist and H(—0) = H(+0) due to Hermite property. Notice that the operator 3{. r 
has a kernel from the class C([0, r] 2 ), r > 0. 
Assume that we have the strict inequality 

1 + 'Kr > (32) 

for any r > 0. Then, there is the resolvent kernel T r (t, s) with nice properties (see 
Lemma [272]) such that 



T r (s,t)=T r (t,s) (33) 
IY(f, s) + I H(t- u)T r (u, s)du = H(t - s), (34) 



T 

r r (t, s) + J T r (t, u)H(u - s)du = H(t - s), < s, t < r (35) 



We emphasize that the last two identities should be understood as equalities for 
functions from C([0, r] 2 ) class. In the meantime, if H(±0) 6 R, then H(x) is 
actually continuous at and by Lemma l2. li the kernel T r (x, y) is continuous on the 
diagonal as well. 

Let us introduce the following "continuous polynomials" 

r 

P(r 7 A) = exp(iAr) — J T r {r, s) exp(i\s)ds (36) 



P« (r, A) = 1 - / r r (s, r) exp(iA(r - s))ds) (37) 
o 

Notice that function P(r, A) is of exponential type exactly r, and P* (r, A) is of 
exponential type not greater than r. 

Formulas (|36p and (|37[) can be easily explained. They are quite natural and have 
analogs in the OPUC theory. Let us consider positive finite measure (i(9) on the 
unit circle. We will denote the inner product of two functions / and g in L 2 (dfi) 
by (f,g)fj,- Let the sequence of moments be 

7T 

exTp(—in9)dii(6) , n S Z + 

Let {ej}, ej — [0, . . . , 0, 1, 0, . . . , 0]*, (j — 0, . . .) denotes the standard orthonormal 
basis. Consider the following Toeplitz matrix 



(38) 





CO 


Cl 


c n 


T — 


Cl 


c 


■ c„_i 






Cn-l • 


c 
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If D n — detTn, then one can easily show that 



Pn(z) = 



c ci 


■ Cn-l 


1 


ci c 


■ Cn-1 


z 


Cn-l C n _2 


C 


z n-\ 


Cn Cn—1 


Cl 


z n 



(39) 



are monic polynomials of degree n (i.e., the coefficient in front of z n is 1) orthogonal 
with respect to d/i. In formula (f3T>)) . exponents exp(iAr), r > 0, A 6 C play the role 
of z n , n G Z + , z G C. The exponential type r is an analog of the integer index n. If 
one introduces the truncated discrete Toeplitz operator T„ given by the matrix (f38|) . 
then P n (^) is the last component of the following vector D rl £'~^; 1 T~ 1 [l, z, . . . , z n ] 1 . 
That follows from Kramer's rule and (f3"9"| . Besides Kramer's rule, there is the 
following algebraic explanation of orthogonality relation. Consider the last com- 



ponent of the vector T n 1 [l, 



z,. 



i.e. 



z,. 



l ]*,e„). Assume that 



we have two polynomials of degree not greater than n: A(z) = a n i 
B(z) = b n z n + ... + b . Then 

(A, B) a = (T n a, b) 

where 7 n is obtained from T n by conjugating all elements, a — [ciq, 
[bo, . . . , b n ] 1 . Therefore, 



(40) 
b = 



{*>, {l-^l, z,..., z n ]\e n )) u = (z>, ([1, z, . . . , z n ]\7-'e n )) n = (T„e„T„ e n ) - S jn 

Now, consider the function f r (x) = cxjp(iXx) on an interval x G [0, r]. Then, 
P(r, A) is the value of the function (1 + Jf,,) -1 /,. at the point x = r. Because the 
spectrum of the truncated continuous Toeplitz operator !H r always contains 0, we 
need to take (1 + IK,,)^ 1 , rather than in the discrete case. This normaliza- 
tion makes the function P(r, A) "monic"- it has exp(iAr) term. In the theory of 
orthogonal polynomials on the unit circle, there is a natural procedure that maps 
any polynomial of degree n to its "reciprocal" . It is defined in the following way 



t](a n z n 



ao) 



a z ■ 



(41) 



or 



[*](A n (z)) = z n A(2-i) (42) 

Notice that the degree of the polynomial might decrease under this operation. For 
example, in the space of polynomials of degree not greater than n, [*](z n ) = 1. The 
natural analog of [*] operation for the function f r (X) of the exponential type r is 
given by the following formula 



(43) 
(44) 

In the discrete 

case, polynomials P n (z) are orthogonal with respect to d/i. In the continuous 
setting, {P(r, A)} turns out to be orthonormal family in L 2 (da) with da from the 
integral representation (|27|) . Orthogonality is understood in the usual sense- just 
like for the Fourier transform. Thus, rather than in the discrete case, continuous 
monic polynomials are already normalized. 



M(/r(A))=exp(iAr)/ r (A) 
Function P*(r, A) is a continuous analog of [*](P n (z)), i.e. 

P.(r,A) = [*](P(r,A)) 
It is clear from the formula (f3T>)) and identity T r (s, t) — T r (t,s). 
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Recall the factorization results from the previous section. Notice that since 
I + "K r > for any r > 0, the Theorem 12. H is applicable. Let us fix any R > and 
consider the factorization (| 1 7|) : I + 3€r = (I + £)(! + IX), where the lower-diagonal 
L has kernel L(x,y) and the upper-diagonal U has kernel U(x,y), < x,y < R. 
Since 5£ is Hermitian, L = U* . If I + L = (I + V-)" 1 and I + U = (I + V+) _ \ 
then 

Lemma 3.2. The following formula is true 

r 

exp(iAr) = P(r, A) + J L{r, s)P(s, X)ds, < r < R (45) 
o 



Proof. Take any R > 0. By (|19p. equation (|36)) can be written as 

P(r, X) = (I + V-) exp(iAr), <r <R 

and the both sides are regarded as functions in L 2 [0, R]. Since I + L = (I + V_) _1 , 
we have the statement of the Lemm£0. □ 

In the discrete case, the following formula is true for any n G Z + 

Z n = P n {z) + ln,jPj(z), lj,n G C 
3=0 



The Lemma 13.21 is the continuous analog of that representation. Now, we are ready 
to prove the following 

Theorem 3.2. The following map is an isometry from L 2 (R + ) into L 2 (R,da) 

oo 

: f(r) .— > (0/) (A) = / f(r)P(r, X)dr (46) 
o 

In other words, 

oo oo 

r \(0f)(X)\ 2 da(X) = / \f(r)\ 2 dr (47) 



Integral in is understood in the L 2 — sense. 

Proof. The following is true for any f,r£l 

r—t oo 

^+ / (r-t-s)H(s)ds = iP(r-t) + [ ( 1 + - «p(iA(r - ^ ^ 



l + A 2 ""J X 2 

-oo 

Multiply this equality by f"(t) (f{t) G C§°(M)) and integrate by parts. We have 
the following formula 



oo oo 



f(r)+ / H(r-t)f(t)dt = 



/(t)exp(iA(r - 



da(X) (48) 



1 Kernel L(x,y) should be regarded as a transformation kernel from one basis to another 1511 
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For f(r), 



f{r)+ J H(r-t)f(t)dt= J exp(iAr) J f(t) exp(i\t)dtda(\) (49) 

— OO — OO — O 

Consequently, we have the following analog of PU[) : 



OO oo 



\f(r)\ 2 dr+ J f(r) j f(t)H(r-t)dtdr = 

— OO —oo 

If f £ L 2 [0, r], then flSOJl implies 



f{t) exp(i\t)dt 



((1+M r )f,f) = 



f(t)exp(i\t)dt 



da(X) 



da(\) 
(50) 

(51) 



by simple approximation argument. Now, let us use Lemma 13.21 for the interval 
[0,r]: 



f(t) exp(i\t)dt = ((I + H)P(t, A), /(t)) za[0jr] = (P(t, A), (J + £,*)/(t)) La[0ir] 
Therefore, (|5T|) reads 



\\g\\ 2 = ((I + i:)- 1 (I + M r )(I + U)- 1 g,g) = 



g(t)P(t,X)dt 



da(X) (52) 



where g — (I + It)/. Since I + U is invertible in L 2 [0,r], ([52^1 holds for any 
g £ L 2 [0,r]. Since \\g\\ = \\g\\ and a number r was chosen arbitrarily, is isomctry 
on L 2 {M+). □ 

As a simple corollary of the Theorem 13.21 and the polarization identity we get 
the following Lemma. 

Lemma 3.3. For any f(r) £ L 2 (E+) 



/(r)= / P(r,X) (Of) (X)da(X) 



where the equality is understood in the L 2 (M. + ) sense. 

Remark 3.1. Notice, that is not necessarily a unitary map. In the simplest 
case H{r) = 0, P(r, A) = exp(iAr), <r(A) = A/(2tt) and the range of is H 2 (R) C 
L 2 (R,da). 

Now, let us obtain the differential system for P(r, A) and P*(r, A). In the discrete 
case, we have 



P n+ i{z) = zP n (z)-a n P*{z), P (z) = l 
P:+i(z) = P*{z)-a n zP n {z), PS{z) = l 



(53) 
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where P^(-z) = [*]P n (z), a n are the so-called Verblunsky coefficients (Geronimus 
coefficients, Schur coefficients, circle parameters, or reflection parameters). If one 
starts with arbitrary positive finite measure with infinite number of growth points, 
then the corresponding a n G P. Conversely, any sequence a n G D yields the unique 
probability measure with infinite number of growth points. 

Let us prove one property of the resolvent kernel which would yield differential 
equations for P(r, A) and P*(r, A). It holds only for the integral operators with the 
"displacement " kernel 

Lemma 3.4. IfT r (s 1 t) is the resolvent kernel for "K r , then 

T r {s,t) = T r {r-t,r-s) (54) 
Proof. The following relation holds 

where the "flip" operator 3> is defined: 3>/(x) — f(r — x) and 3i r is an integral 
operator with the displacement kernel H(x). Then, 

Writing this down in terms of the resolvent kernel gives (|54|) . □ 

As a simple corollary, we also get the following formulas for P and P* which 
were originally used by Krein 

P(r, A) = exp(iAr) f 1 ~ / Tr ( s ' °) expHAs)dsj (55) 

r 

P*(r, A) = 1 — J r r (0, s) exp(i\s)ds (56) 



An analog of the relations (|53|) is given by the following statement. 
Theorem 3.3. The following equations hold 

/ P'=i\P-AP*, P(0,A) = 1, 
\ Pi = -AP, P*(0,A) = 1 

where 

A(r)=T r (0,r) (58) 
Proof. Differentiate ([56]) and use Lemma [2?2l to get 



Pi(r,X) = -r r (0,r)exp(iAr) + J T r (0,r)T r (r,s)exp(iXs)ds = -A(r)P(r,X) 

o 

where we used the definition of P(r, A). Equation for P(r, A) can be obtained from 
the equation for P* (r, A) and (|44|) . □ 

Definition 3.2. System ((57)) is called the Krein system. 

Obviously, coefficient A(r) in the Krein system is an analog of Verblunsky coef- 
ficients. 
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Lemma 3.5. Under regularity conditions \30\) . we have A{r) € C[0, oo). We also 
have A(0) = H(-0) = H(+0). 

Proof. The continuity of A(r) follows from Lemrna [2~2"1 The equality A(0) = H(—0) 
can be obtained from (l35l) where r — > 0. □ 



Notice that H(x) G C(R) iff H(0)- real iff A(0) is real as well. 
Let us consider the model ("free") case. If <r(A) = \/(2ir) in (127p . then /3 = 0, 
s(t) = |i|/2, #"(*) = 0, r r (s, *) = 0, A(r) = 0, P(r, A) = exp(iAr), F*(r, A) = 1. 

Lemma 3.6. 77ie following is true 
1) Christoffel-Darboux formula: 



P* (r, A)P» (r, /it) = P(r, A)P(r, /it) - *(A - p) J P(s, A)P(s, /x)c2s, A,/i6C (59) 

o 

2) For AeC, 

P(r, A) = cxp(iAr)P,(r, A) (60) 

3) P* (r, A) does not have any zeroes in C + 

4) P(r, A) has zeroes in C + only 



Proof. To prove Christoffel-Darboux formula, multiply the first equation in (|57|) by 
P(r, /i). Multiply both sides of 



Pl(r,») = -A(r)P(r,ti) 
by P*(r, A) and subtract two identities. One has 



P(r,A)P(r,/i) -Pi(r,M)P*(r,A) = iXP(r, X)P(r, fx) (61) 

Write the same equation with A and /i interchanged. Take conjugate and add to 
(|6T|) . One gets ([59j. 
2) repeats (JUJ). 

Take A e C+ and /i = A. Then, ((59]) guarantees that |P*(r,A)| > 0, so P*(r,A) 
has no zeroes in C + . Assume that P*(r, A) = for some A e M. Then, by 2), 
P(r, A) = 0. However, functions P(p, A), P*(p, A) solve the problem (|57|) and then 
must vanish for all p > 0. In the meantime P(0, A) = P*(0, A) = 1. This contradic- 
tion shows that P* (r, A) has no zeroes in C+ . 

4) follows from 3) and flS}. □ 

Assume that for fixed r we know P(r, A) as the function in A. The natural 
question is whether we can find P(p, A) and A(p) for all < p < rl The answer 
happens to be positive. Notice that since P(r, A) is given, we know the values of 
r r (s, 0) for all s£ [0, r]. It easily follows from ([55]) . 

Consider 

g(t, s) = r r (0, r - *)r r (r - s, 0) - T r (t, 0)T r (0, s),0<s,t<r 
Lemma 3.7. The following formula holds true 

min(s,i) 



r r (*, s) = r r (* - s, 0) + T r (s - *,0) + / g(t -u,s- u)du, < s ^ * < r (62) 

o 

where r r (0, s) = r r (s, 0) = 0, s < for shorthand. 



CONTINUOUS ANALOGS OF POLYNOMIALS 



21 



Proof. Assume that H(x) € C^r] first. Then, T r (s,t) € C 1 (A±). Let us show 
that 

dr r ( s ,t) , ar r («,i) 

= g{s,t), s,teA± (63) 



eft 



<9,s 



We have 



-r r (ii, r)r r (r, si) = -^-r r (ii, si) = T r (r - s x ,r - ti) 
dr dr 



= -T r (r - s ll r)T r (r,r - t x ) + d Xl T r (r - s 1 ,r - ti) + d X2 T r (r - s x ,r - h) 

Taking t = r — t\, s = r — s\, we get ([M]) . 

Solving the linear first order equation, we obtain (|6^|) . Now, consider H(x) E 
C[0, r) . We can approximate H{x) by the sequence of Hermitian functions H' Jl ' (x) e 
C 1 ^, r] in C[0, r] norm. Since 1 + "K r > 0, 1 + %r 1 " 1 > 0, the corresponding kernel 
r£ (s, t) converges to T r (s, t) uniformly as elements of C([0, r] 2 ). Then, apply 
to r£ n '(s, t) and take n — > oo. □ 

Now, assume we are given r r (s,0) for all < s < r. From (I62|) . we know the 
resolvent kernel T r (s, t) for all < s, t < r and then can find H(x), say, from (|35|) . 
The accelerant -ff(x) on [0, r] defines A(x) for a; G [0,r] by construction. 

There is yet another way to find H(x), < x < r from P(r, A). It is given by the 
following 

Lemma 3.8. The following representation is true 

oo 

1 / ^ (64) 



1+ / H r (s) exp(i\s)ds, A € 



|P*(r,A)P 

— oo 

where H r (s) is Hermitian function, H r (s) G L 1 (R), and H r (s) = H(s) for \s\ < r 



Proof. Formula (|56|) . 3) in Lemma [3~6l and Levy- Wiener theorem yield existence of 
function H r G L 1 (R). Moreover, H r is continuous on K except for the points 0, ±r, 
where the left(right) limits exist. That is clear from the corresponding integral 
equation. Let us show that this H r (s) coincides with H(s) for \s\ < r. We have 



— 1 = P*(r, A) f 1 + / H r (s) exp(i\s)ds - 1 

'*(r,A) \ J J 



P, 



The left-hand side belongs H 2 (R). Therefore, 



P,(r,A)|l+ J H r (s) exp[iXs)ds J - 1 



which gives 



r r (s,0) + / iJ r (s-u)r r (w,0)dw = H r {s), < s < r 



Formula (|64|) proves that 7J r defines an integral operator !K r on L 2 [0,r] and / + 
!K r > 0. Denote its resolvent kernel by T r (s, t). The last equation shows that 
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IY(s,0) = r r (s,0),0 < s<r. Due to Lemma [33 T r (s,t) = T r (s,t) for all < 
s,t < r. So, H r (s) = H(s) for < s < r. Since H r is Hcrmitian, we obtain the 
statement of the Lemma. □ 

Thus, if we know P(r, A), we know P*(r, A) as well and can find H(x) for \x\ < r 
using the previous Lemma. 

Remarks and historical notes. 

Continuous analogs of polynomials orthogonal on the circle were introduced by 
M.G. Krein in the paper [33] but no proofs were given. We filled this gap. Lemma 
GO is in [46], see also [24], p. 100. 

If one is given the function H (r) £ C[0, R], Hermitian and such that I + !Hr > 0, 
then the Krein system can be well-defined on the interval [0, R]. In the meantime, 
the question of orthogonality with respect to some measure gives rise to certain 
continuation problem [3S] we do not want to address here. 



CONTINUOUS ANALOGS OF POLYNOMIALS 



2:S 



4. Krein systems 

In the previous section we learned that any accelerant H(r) gives rise to (|57|) . the 
system of ODE called the Krein system. But it makes sense to study this system 
per se. In this section, we will show that one can start with the Krein systems and 
then define the accelerant H(r) and measure <r(A) uniquely. 

Consider the system 

X' = VX (65) 

with X(0, A) = I, 



V = 



iX 

-A(r) 



■Mr) 




(66) 



Matrix V has very special algebraic structure and it should imply very special 
properties for the fundamental (transfer) matrix X(r). Assume first that A(r) £ 



The first obvious result is 
Lemma 4.1. We have 



detX(r) = exp(zAr) 



(67) 



Proof. Indeed, 



det X(r) — exp 



Tr V(t)dt 



exp(zAr) 



□ 



Consider the signature matrix 



J = 





-1 



(68) 



Definition 4.1. The matrix M is said to be J- contraction if M* JM < J. 

Definition 4.2. The matrix M is called J- unitary if M* JM = J 

Later, we will need the following algebraic 

Lemma 4.2. If M is J -unitary, then |detM| = 1, and Af _1 ,M* are J -unitary 
too. If M is J- contraction then M* is J- contraction also 

whose proof is given in the Appendix. 

The signature matrix J defines the corresponding indefinite metric. For general 
properties of these spaces and operators acting on them, see [3TJ [33] . 
The next very important algebraic observation is 



V*(r)J + JV(r) 



~2ImA 



1 




(69) 



for any r > 0. 



Theorem 4.1. The matrix X is J- contraction for A G C + and is J- unitary for 

AeR. 
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Proof. Consider Y = X*JX. Then, fSl) yields 



Y' 



~2ImAX* 



1 




X, Y(0, A) = J 



Therefore, for any / S C 2 , we have 

(Ff, /) = (.//,/), A el 
(17,/) < (J/,/),AeC H 
(Yf,f)>(Jf,f),\eC- 

which implies the statement of the theorem. 



(70) 
(71) 
(72) 
□ 



In this section, we define the functions P(r, A), P*(r, A) as solutions of equation 
(p3]) corresponding the Cauchy problem P(0, A) = P*(0, A) = 1. Consider also two 
functions P(r, A) and P*(r, A) such that the vector P(r, A), — P* (r, A) solves (p5|) 
and satisfies initial condition P(0, A) = 1,— P*(0, A) = —1. The simple calculation 
shows that 



X(r) 



1 


P + 


P 


p 


- p 


i 


P 


P 




' i i 


2 


P* - 


i 3 * 


p* 


fP« 


~ 2 


P* 


-i 3 * 




i -i 



The sign " - 
Notice that 



Then, 



(73) 

in the definitions of P, (r, A) was chosen for the following reason. 
JV(r)J 



i\ A(r) 
A(r) 



(74) 



Lemma 4.3. // X solves equation &65\). then JXJ solves the same equation but 
with coefficient A{r) having an opposite sign. 

Proof. Multiply ^ from the left by J. Then use ([HD and identity J 2 = I. □ 

Corollary 4.1. The vector P(r, A), P*(r, A) satisfies the same initial conditions at 
zero as P(r, A), P* (r, A) but solves system Ii65\) with A(r) having an opposite sign. 
This system is called the dual Krein system. 

Instead of dealing with the transfer matrix X(r,X) which solves (|65p . we will 
first study functions P, P*,P, and P*. Below, we list some simple properties. 



Lemma 4.4. 



1. For any A £ C and r > 
P(r, A)P* (r, A) + P* (r, A)P(r, A) = 2 exp(iAr) 



(75) 



2. All statements of Lemma lS. 6\ are true for both P(r, A), P*(r, A) andP(r, A), P*(r, A). 

3. For A € C + , we have 



Re 



Pr x (r,A)p;(r-,A) > |P.(r,A)| 



(76) 



and the last inequality is equality for real A. 



Proof. Part 1) follows from (f67|) and (173]) . 

To show 2), notice that Christoffel-Darboux formula is the direct conse- 
quence of the differential equations (jST)) . Then, ([60)) holds because the functions 



exp(iAr)P* (r, A), exp(iAr)P(r, A) solve the same Cauchy problem as P(r, A), P* (r, A) 
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do. The statements about the zeroes of P, P, can be proved in the same way as it 
was done in Lemma \TM The analogous results for P, P* follow from Corollarv l4.ll 
To prove 3), notice that by Lemma [4.21 X* is J-contraction for A £ C + and 
J-unitary for real A. Writing XJX* < J in terms of P, P*, P, P, we get 



< 



1 

-1 



= J 



1 PP + PP PP ± - PP* 

2 [ P*P- P,p -(P*P* + P*P*) 

Element (2, 2) gives 

P*P*+PP>2 (77) 
which implies 3). For real A, we get equality because X* is J-unitary. □ 

Lemma 4.5. If A(r) is real and A = 0, the exact solution can be obtained, i.e. 



X(r,0) = 



cosh 



fA(t)dtj sinh (- jA(t)dtj 



sinh ( - J A(t)dtj cosh ^- / 



Proof. The proof is a direct calculation. 

Lemma 4.6. The following estimate is true if X € 



a 



exp 



\A(s)\ds 



< \P*{r,\)\ <exp 



\A{s)\ds 



Proof. The second inequality easily follows from the differential equations for P 
and P*. The first one is then immediate from (1771). □ 



Now, that we studied the general properties of system (|55|) . let us show that for 
any A(r) G C[0, oo), there is the unique accelerant H(r) £ C[0, oo) that generates 
it. 

The following result says that solutions of Krein system are indeed continuous 
polynomials. 



Lemma 4.7. For any r > 0, we have the following formulas 



T 

P(r, A) = exp(iAr) — J A(r, s) exp(i\s)ds 



(78) 



P*(r, A) = 1 - / j4(r, s) exp(iA(r - s))ds 



(79) 



where function A{r, s) is continuous in s and r : < s < r < oo. 
Proof. Consider Q = exp(— i\r)P. We have the following equations for P and Q 
Q' = -cxp(-ar)AP, Q(0,A) = 1 



Pi = -cxp{i\r)AQ, P»(0,A) = 1 
The corresponding integral equations are 



Q(r, A) = 1 - / cxp{-iXs)A{s)P« (s, A)ds 
o 



(80) 



(81) 
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r 

P*(r, A) = 1 - J exp(i\s)A(s)Q(s, X)ds 



Let us find the solutions to (|57|) in the following form 



(82) 



P(r, A) = exp(iAr) — / A(r, s) exp(i\s)ds, P* (r, A) = 1 — / B(r, s) exp(iA(r — s))ds 
o o 

(83) 

where A and P are continuous function. Then, part 2) of the Lemma 13.61 yields 



B(r, s) = A(r, s). Plug (JS3J) into JgTJ to get the equation for A(r, s) 



A(r,t) = A(r - t) - / A(s) A(s, r - t)d 



(84) 



Fix any positive R. In the triangle = {0 < t < r < R}, consider the operator 



[Of](r,t) = / A(s)f(s,r-t)ds 



(85) 



Let us shows that O is Volterra in C(Ar). That would allow us to solve (|84|) 
uniquely. 

The following inequalities hold true 



(2n) /l ( r , t ) 



(2n-l) 



M) 



< 



< 



n! n! 

ll^ll^ll/llcCAajCr-t)"- 1 ^ 



(86) 



(87) 



n!(n- 1)! 

Let us prove them by induction. For n — 0, 1, the estimates are obvious. Assume 
that (|86|) is true for n. Then, we get 

II^I| 2 C^H]II/IIC(A H ) 



< 



< 



||A|| 2 ™ +1 



(n!) 2 

c[o,B]\\f\\c(A R ) 
n\(n + 1)! 

Assuming that ([57]) is true for n we obtain 
II 41 ' 2n 

(2n) / (r,t) 



-(r-t) n / [s - (r - t)] n ds < 

r-t 



(r - i) n i 



n / n+1 



Jl^ll 2 cV]li/ilc(A R ) 



< 



n!(n- 1)! 
\\A\\% , R] \\f\\c(A R) 



(r - t) n t" 



(n!) 2 

Therefore, ||0 (,l) || -> as n — > oo and O is Volterra. Notice that the estimates 
obtained above prove convergence of the series obtained by the iteration of (|84f . 
Therefore, the solution A(r, t) is continuous inO <t < r < oo. The corresponding 
Q and P* solve integral equations {EH and (|82]), So, P and P* from (Jj8| and [[79]) 
solve the Krein system. □ 
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Remark 4.1. From (l84|) . we have an identity 



A(r,0) = A(r) 

Theorem 4.2. For any Krein system |57[ ) with 

A(r) eC[0,oo), (89) 

there is the unique accelerant H{x) which generates it and satisfies \30\) . Con- 
versely, any accelerant satisfying \30\) gives rise to the Krein system for which (8ty) 
is true 

Proof. The converse statement follows from the construction done in previous sec- 
tion. 

Now, let us find an accelerant that generates the given Krein system. The clue 
is given by Lemma [3.81 Consider the function P^ 1 ^, A)P (r, A). We know that 
P*(r, A) does not have zeroes in C+. By Lemma 14771 both P, and P* are continuous 
polynomials, i.e. 

r r 

P* (r, A) = 1 — J A(r, r - s) exp(i\s)ds, P* (r, A) = 1 - J A(r, r - s) exp(i\s)ds 
o o 
Therefore, Levy- Wiener theorem yields the representation 



P' 1 ^, A)P* (r, A) = 1 + 2 / H r (s) exp(isA)ds, (90) 



with H r (s) € i 1 (K + ). H r (s) is continuous on [0, r] and [r, oo] but the right and 
the left limits at r are not necessarily the same. From differential equations for P* 
and P* and Lemma l4~4l (parts (2) and (3)), we have 

V(r,A)P.(r,A)]= 2 ^ A ( ; rA) ,Ae R (91) 
Consequently 



d 
dr 



P-\r 2 , A)P*(r 2 , A) - P" 1 ^, A)P*(n, A) = J ^f*^ da, < n < r 2 

ri 

and that implies 

H ri (s) = H r2 (s) = H(s) (92) 

for < s < r±. Therefore, the function H(s) is well-defined and continuous on 
[0,oo). 

Therefore, if we let H r (-s) = H r (s), H(-s) = H(s),s > 0, then H(s) is Her- 
mitian and satisfies ([30]) . Now, let us show that it actually generates the Krein 
system with given coefficient A(r). 

Indeed, from Lemma 14.41 part 3, we have 

oo 

' = Re 



\P*(r,X)\ 



p- L (r,\)P*(r,\) =1+ / H r (s)exv(is\)ds, (93) 
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Notice that the last identity yields 



\h{x)\ 2 dx + J J H r (x-y)h(y)h(x)dydx>0 (94) 



for any h £ C^°(0,oo) and the inequality is strict for any nontrivial h. Indeed, 
one needs to rewrite (|94p in terms of Fourier transform. Then, (|9"2")l shows that 
"K r + I > for any r > and H does generate the Krein system with coefficient 
A^\r) that satisfies ([89]). Now, let us prove that A^(r) = A(r). Denote the 
solutions of §7$ with coefficient A {l \r) by pW and P* (1) . But P* (1) (r, A) = P,(r, A) 
for any r > 0, A £ C. Indeed, from Lemma applied to Pp^ and (|9"5)) , we get 



1 

1 



|P*(r,A)| 



i 1 



P* (1) (r, A) | 



2 



Then, by Lemma [17.41 from Appendix, we get P* (r, A) = P*(r, A). Therefore, 
pW(r, A) = P(r, A) for all r > and A^)(r) = A(r). □ 

Remark 4.2. Notice that the values of accelerant on [0, R] depends solely on the 
values of A(r) on [0, R] and vice versa. 

The application of Levy- Wiener theorem to P7 1 (r, A)P*(r, A) shows that 

2H r (+0) = -A{r,r)+A(r,r) 

and therefore 



H(+0) = lim P r (+0) = [-A(0, 0) + A(0, 0)]/2 = A(0) 

r— »0 

where we used f|88j> . Therefore, if A(Q) G M, then H(x) is continuous at and 
P(0) € R. 

Theorem 14.21 establishes a one-to-one correspondence between continuous A(r), 
defined on R + , and continuous accelerants for which (|32p is true. But what happens 
to a map {H(x) — > -A(r)} if (l32|) fails at a finite point? For OPUC, if the measure 
has only k growth points, then Du-i ^ 0, Dk — 0. The corresponding \a,j\ < l,j = 
0, . . . , k — 1, |afc| = 1. For the Krein system, the situation is similar. Assume that 
1 + 3i r > for all r < R and ker(7 + !Kr) 7^ 0. Following argument given above, 
one can construct A(r) £ C[0,R). Vice versa, given A(r) £ C[0,R), we can define 
H £ C[0, R) such that ([32]) holds up to R. But as long as ker(7 + JC fl ) ^ 0, A(r) 
blows up as r approaches R from the left. More precisely, this process is governed 
by a pair of simple (and clearly very crude) estimates 

\A(r)\ < ||r r (t,0)|| <7[Oir ] < ||(/ + ^r)- 1 ||c[o J r]l| J ff||a[o,r] (95) 

and 

\\(I + X r )-i\\ c[0;r] <l + Cr||A|| C[0>r] (l + r||A|| C[0 , r . ] )e X p[Cr||A|| c[0;r] ] (96) 

The first estimate easily follows from with s = 0. It shows that A(r) can not 
blow up unless ((32|) fails at a finite point (here we also assume that H{x) £ C[0, R]). 
One can get (|96ll from the following inequalities. 

||(/ + J{ r )- 1 || c[0 . r] = ||7-r r || c[0 . r] < l + r max \T r (s,t)\ 

0<s,t<r 
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To estimate the last maximum, we follow the proof of Lemma 14.71 (estimates ([86]) 
and ([87]) ). This gives us the following bound 

max |r,.(0,s)| < C||A|| 0[ o lT .]e3tp[Cr||A||(7 [ o, r ]] 

sG[0,r] 

because A(r, s) from Lemma 14771 is actually equal to T r (r, s) — T r (r — s, 0) by (|54l) . 
Using Lemma l3.7[ we obtain an estimate 

max |r r (s,t)| <C||A|| C [o )r ](l + r||A||cro,r])exp[<7r||A||c.[o jr ]] 

9<s,£<r L 

which yields (|96p . Obviously the left-hand side of ([96)) is non-decreasing in r. If it 
blows up at a finite time (i.e. ([32]) fails at a finite time), then A blows up at the 
same point as well. Comparing the OPUC and Krein systems, we see that infinity 
(for Krein systems) plays the role of T = <9B for OPUC. 
Remarks and historical notes. 

The problem of constructing accelerant from the Krein system with locally inte- 
grable coefficient A(r) was solved by Rybalko in [61]. In [61], estimates on A(r,t) 
from the Lemma 14.71 are a bit stronger than what we obtain. Some generalizations 
of Krein systems were considered by L. Sakhnovich in [65]. 
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5. Accelerant and A(r) are from L 2 oc (R + ) class 

In this section, we will show that the accelerant from L 2 oc (IR) class generates 
the Krein system with A(r) E Lf oc (M. + ) and, conversely, the Krein system with 
A{r) g L 2 oc (]R + ) generates the accelerant H(x) £ L 2 oc (R). Moreover, this highly 
nonlinear map is homcomorphism in L 2 oc , i.e., in L 2 [0, R] for any R > 0. Wc 
will prove that all statements from the previous two sections find their analogs for 
Lf oc -case. 

First, consider an accelerant H(x) € L 2 [0,r] for any r. Then, / + 3i r > and 
T r (x,y) is well-defined as L 2 ([0,r] 2 )- function. We also have 

r 

T r (t, s) + J H(t - u)T r (u, s)du = H(t - s) 
o 

Fix any s G [0, r] in the last equation. Then, H(t — s) is continuous in s as L 2 [0, rj- 
function in t. Therefore, T r (t, s) is continuous in s in L 2 [0, r] norm with respect to 
the first coordinate. It is also Hermitian function so the same is true for s and t 
interchanged. Let g r {s) = r,.(s, 0) for s S [0, r) and g r (s) = for s S [r, R]. Due to 
the formula T r (s, 0) = (/ + Jf r ) _1 i?, we have the continuity of g r in r with respect 
to L 2 [0. R] norm. 

The special displacement structure of the kernel of !H r allows to generalize The- 
orems O and O 

Theorem 5.1. For Hermitian H(x) S L 2 [—R,R], the operator I + "Kb, admits 
factorization p7\) if and only if I + "K r > for any r S (0, R\. In this case, 



V+(x,y) = -T y (x,y),x<y , g7 > 
V-(x,y) = -T x (x,y),x>y 

where T r (x 1 y) denotes the resolvent kernel of I + "K r - 

Proof. Now, assume that I + K r >0 for any r € (0, R]. Define V± by d97]). Notice 
that operators V± are well-defined. It follows from the representation 

xx R 

[V-f](x) = -f T x (x,y)f(y)dy = - J T x (s,0)f(x - s)ds = - J g x (s)f(x-s)ds 



which shows that V_ is actually bounded from L 2 [0,R] to L°°[0, R]. Analogous 
formula is true for V+. These operators also have Hilbert-Schmidt and Volterra 
properties. Approximate H(x) by H^(x) S C[—R,R] in L 2 [—R,R] norm and 
apply Theorem 12.11 For each n, formula (TT5)) is true. Moreover, V±^ (x, y) — > 

V±{ x , y)i r^ 1 "* (x, y) — > rfl(x, as n — > oo and convergence is in L 2 ([0, R] 2 ). Taking 
n — > oo, wc get (fT5|) for G and V±. That implies the needed factorization. The 
converse statement is simple and repeats the argument in Theorem 12. II □ 



The analog of the Theorem 12. 21 can also be easily proved in the same way giving 

Theorem 5.2. For Hermitian H(x) € L 2 [— 2R, 2R], the operator I + "Kp admits 
factorization H20\) if and only if I + !K r is invertible in L 2 [—r, r] for any < r < R. 
In this case, 

V+(x,y) = -t lyl (x,y), (x,y) € (l+ , gg s 



V-(x,y) = -T\ x \(x,y), (x,y) G O- 
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All functions that are used in the definition of continuous polynomials are now 
well-defined as elements of L 2 [0,r]. Indeed, F r (r,s) and IY(0, s) are both from 
L 2 [0, r] and we can consider the corresponding continuous polynomials P(r, A), P*(r, A). 
Moreover, Lemma \3. 2 1 and Theorem 13.21 are true for H £ L 2 oc (IR) as well. Indeed, 
their proofs were based on the factorization of Fredholm operators (the analog of 
which we just proved). 

Now, let us define the coefficient A(r) and show that A(r) £ L 2 oc (R + ) and that 
equations ((57)) are true. 

Consider an accelerant H{x) £ L 2 [—R,R] and approximate it with H^ n \x) £ 
C[—R,R] in L 2 {— R, R] norm. Then, each H^ n \x) generates the Krein system on 
[0,R] with A<™)(r) £ C[Q,R). We have 



A(«)(r) = H {n) (r) - J H [n) {r - u)T { r n) (u,0)du -► H(r) - J H(r - u)T r (u,0)dv 
o o 
in £ 2 [0, R}. The conjugate of last function will be denoted by A(r), i.e. 



A{r) = H(r) - / H(r - u)T r (u, 0)du 



o 



Notice that the second term is continuous function in r. Therefore, all singularities 
of A(r) and H{r) coincide. In particular, A(r) £ L 2 [0,R]. For any A £ C, the 
polynomials P< n ) (r, A) , Pi™ } (r, A) converge to P(r, A) and P*(r, A) uniformly in r £ 
[0, R\. Notice that equations (|5T[) for P^ and P^ are equivalent to the system of 
integral equations 



P (n) (r, A) = 1 + / i\P (n) (s, A) - AW (s)P* (n) (s, A) 



P* (r, A) = 1 - / A (Tl) (s)P (n) ( S7 A)ds 



Taking the limit n — > oo, we see that P and P* satisfy the corresponding equations 
that are equivalent to (|57[) . The proofs of Lemma 13. 61 and Lemma 13751 work for L 2 oc 
case. Lemma 13.71 can be shown by an approximation argument. 

Now, following the arguments from the Section 01 we consider A(r) £ i 2 oc (R+). 
We have 

Theorem 5.3. For any Krein system |57[ ) with 

A(r) £ L 2 oc (K + ) (99) 
there is the unique accelerant which generates it and satisfies 

H(x) £ i 2 oc (R) (100) 

Conversely, any accelerant satisfying il00\) gives rise to Krein system with A(r) £ 
i 2 oc (M + ). This map is a homeomorphism. 

Proof. As we just showed, the L 2 oc accelerant does generates the Krein system with 
L 2 oc coefficient. Now, let us start with A(r) £ Lf oc (R + ). Notice that all arguments 
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from the proof of Theorem l4.2l are valid as long as we have an analog of Lemma 14771 
Therefore, we just need 

Lemma 5.1. If A(r) S Lj 2 oc (K + ), and P(r, A) , P*(r, A) are defined as solutions to 
the Krein system, then the representations |7ff| ], |7ff| ) ZioW wii/i A(r, s) € L 2 [0, r]. 

Proof. We need to consider the operator O given by (|85|) and show that it is Volterra 
in the space of functions f(r,y),0<y<r<R such that 

1/2 

I/O", y)| 2 rfy 



|oo. 2 = sup 



< oo 



Let us prove by induction that 



(n)f 



M) 



< 



t,L 2 [0,r] 



\A(s)\ds 



l/IU,2/(n!) 



(101) 



For n — 0, the statement is elementary. Assume that this estimate is true for n. 
Then, for n + 1, we have 



Q(n+1) r 



(r,t)= / A(a) 



(s, r — t)ds 



Application of Minkowski inequality and induction assumption yields 



Q(n+l) j 



M) 



t,L 2 [0,r] 



< ll/lloo^n!)- 1 \A(s)\ 



\A(u)\du 



|/IU,2[(n + l)! 



\A{s)\ds 



n+l 



Thus, we have (| 101(1 and 

P(M)lki 2 [(Vr] < miU2[o, r ] exp [||A|| L i[ , r ]] 

This estimate finishes the proof of the Lemma. 

We are left with proving 
Lemma 5.2. The following estimates are true for any R > 

\\H{x)\\ L 2 { _ RM] < C\\A(r)\\ L 2 [0tR] exp(C\\A\\ L i [0tR] ) 

\\A(r) -H$)\\ L - [0iItl < \\H\\i 2[0jR] \\(I + X R )-%, 2 
and the map A(r) — > H(x) is homeomorphism in Lf oc . 

Proof. From ([90]) and (l9T|) . we have 



□ 



(102) 
(103) 



n 2 0R,A) 



H R (s) e~xp(i\s)ds 
o o 

R R 

A(s) exp(i\s)ds — 

o Lo 



ds 



A(s) exp(iAs) 
n 2 (*,A) 

A(u) exp(i\u)du 



A(s)P( Sl X) 



ds 
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Therefore, using \P(s, A)| = |P*(s,A)| and Lemma l4~6l we get 

R 

l|Bfl||a< II^IU=[o,fl]exp(2||A|| z , 1[0)fl] ) + J ||i4|| i3[0ia] |il(s)|ejcp(2||i4|| il[0ja] )ds 

o 

which yields (|102p since H(x) = Hr{x) for \x\ < R. This argument also shows that 
H(x) depends on A(r) continuously in the Lf oc . 
To prove (|103[) . we use to write 



T r {t, 0) + / H(t - u)T r {u, 0)du = H(t) (104) 
o 

Taking t = r, 

r 

A(rj + J H(r-u)T r {u,0)du = H(r) (105) 
o 

The formula CEl implies IY(t,0) = (/ + "K r )~ l H and, therefore, \\T r (t, 0)|| L 2 [0 ^ < 
IKZ + Jfr)- 1 !^^!!^!!^^]- That yi eld s CHI- Also, the map H(x) — > A(r) is 
continuous in Lf oc . □ 

□ 

We want to mention here that slight modification of the arguments allows to 
prove that the map H{x) — ► A(r) is homeomorphism in Lf oc (R) — Lf oc (R + ) for 
any p > 1. For one direction, one just have to iterate equation (|104|) sufficiently 
many times to achieve the necessary gain in regularity and then plug it in (|105l) . 
The other direction is straightforward. 
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6. Continuous analogs of Wall polynomials and Schur function. 
Bernstein-Szego approximation 



X(r,X) 



(106) 



Let us introduce the continuous analogs of the so-called Wall polynomials. Con- 
sider the functions 2t(r, A), 58 (r, A), 21* (r, A), 58* (r, A) defined by 

21* (r, A) 58* (r, A) 
58(r,A) 2l(r,A) 

were X(r, A) is the transfer matrix given in (|65j) . By analogy with OPUC theory, 
it makes sense to call 21 and 58 the continuous Wall polynomials. They can be 
rewritten in the following way 

*(r, A) = Ml^M), ^ (r , ,) = ^A)+P(r,A) 

®(r,A) = P *>*- P *> X \x.{r,X) = ^f^X) {m) 

Lemma 6.1. For continuous Wall polynomials, the following identities are true 

1) For A g K, 

| 21 1 2 - 1 58 | 2 = 1,| 21 1 = | 21* |, 1 58 | = 1 58* |,2LB* = S21 (108) 

2) for AeC, 

2l(r,A)2t*(r,A)-58(r,A)58*(r,A) = exp(iAr), (109) 

3) For X g C+, 



1211 



|58| 2 >l,|2l| 



58* 



> 1 



(110) 
□ 



Proof. The proof follows directly from Theorem [JT] or Lemma 

Notice that pTO) implies 58(r, A) 2P 1 (r, A) g B(C+). 

Theorem 6.1. The ratio 58(r, A) 2l _1 (r, A) converges to f(A) g B(C+). T/iis con- 
vergence is uniform over all compacts in C + . 

Proof. Take r\ < r 2 . Consider the Krein system on the interval [r\, oo). Denote the 
transfer matrix from r% to r 2 by X(r\, r 2 , A). Then, in our notations, AT (0, r, A) = 
X(r, A). If we introduce 



X(n,r,X) 



(111) 



a*(r, A) b*(r, A) 
b(r,A) a(r,A) 

then a and b are Wall polynomials for the same Krein system considered on the 
interval [ri,oo). Obviously, Lemma I6TT1 will hold for these functions as well. 
The semigroup relation 

X(0,r 2 )=X(r 1 ,r 2 )-X(0,ri) 

yields 

2l(r 2 ) = b(r 2 )58*(r 1 ) + a(r 2 )2l(r 1 ) 
®(r 2 ) = b(r 2 )2l*(r 1 ) + a(r 2 )58(r 1 ) 



(112) 



For the function 58 21 



5B(r 2 ) = 58(r 1 ) + (bo- 1 )2l*(ri) 
2l(r 2 ) - 2l(r 1 ) + (ba-i)58*(r 1 ) 



(113) 
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and 

g(r 2 ) g(ri) _ (ba~ 1 )(2l(r 1 )Sl.(r 1 )- i B(r 1 )S4r 1 )) 
2l(r 2 ) a(n) 2l(r 1 )[a(r 1 ) + (ba- 1 )Q5*(r 1 )] 

(bQcxpfzAri) 
2t(r 1 )[2l(r 1 ) + (ba- 1 )?8 t (r- 1 )] 1 > 

where we have used (|109[) . 

We have bo -1 G B(C+). Therefore, from (fTT0|) . we get 

|2l(r 1 )[2l(r 1 ) + (ba- 1 ) ( B >t (r 1 )]| > (| 2l(n)| - | «B,(n)|) | »(n)| > (115) 



>l^)i 2 -]^(n)l 2 >i, AeC+ (ii6) 



Then, 



< 2 exp(— ri Im A) 



a(ra) 2l(n) 

That means Q3(r, A) 2l _ (r, A) converges to a certain function f(A) as p — > oo. This 
convergence is uniform in any compact in C + . □ 

We will call f(A) the Schur function corresponding to [0, oo). Notice that f(A) e 
55 (C + ). The following formula is the consequence of (| 1 1 3|) if one takes r 2 — > oo 

g(p,A)+f,(A)g.(p,A) 
HJ a(p,A)+f„(A)B.(p,A) 1 J 

where f p (z) is Schur's function for the same Krein but on the interval [p, oo). 

Not every function from B(C + ) is Schur's function of some Krein system. The 
characterization of that special subclass will be given later but now we just want 
to mention that 

f(A) -> (118) 

if ImA — > +oo. It easily follows from the formula (|117[) with any fixed p and 
relations 

f p (A) G B(C+)Mp, A) — > 1, Q3(p, A) — > 0, 21* (p, A) -> 0, <8* (p, A) - 
as Im A — > +oo. 

In the previous section, we constructed an accelerant from the given Krein sys- 
tem. Then, the measure a and the constant /3 can be found from the formula (|2T[) . 
But there is more direct way to find these data. The next Theorem develops an 
analog of the Weyl-Titchmarsh theory [49] for Krein systems. 

Theorem 6.2. The ratio P*(r, X)P~ 1 {r, A) converges to the function F(X) uni- 
formly in any compact in C + as r — ► oo. This function F(X) has the positive real 
part in C + and allows the following representation 

oo 

F(X)/2 = -0 + i J — L+^— A,® (119) 

— oo 

where da and j3 coincide with those from the formula \27j) . Moreover, the sequence 
of measures 

^= (2^W ^ MA) (120) 
in the weak-{*) sense (analog of Bernstein-Szego approximation). 
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Proof. From (flDT]) . 

P* = 2t + Q3,P» = 2l-<8,P = 2l» +P„P = 2t*-«8* 

Therefore, 

p- 1 (r,A)P st (r,A) = (521 - 93) (21 + S)" 1 = - — ^—j^- 

That shows convergence of P~ (r, A)P* (r, A) to the function 

p(A) = (l-f)(l+f)- 1 (121) 

Again, this convergence is uniform for compacts in C + . Function P(A) has positive 
real part in C + because f(A) G B(C + ). Therefore, P(A)/2 admits the following 
integral representation [4 



1 + At 



F(X)/2 = -il3-ia\ + i / ^ ^7- ^dait) (122) 



where a, (3 £ R, a > and non-decreasing function <r is such that 

oo 

da(t) 



i 2 + l 



< oo 



Recall the way this formula is obtained. If we map A € C + onto z € D by 
conformal mapping z = (A — i)(A + i) , then the function g(z) — 2~ 1 F(i(z + 
1)(1 — z) _1 ) is the Herglotz function in ID. It has a canonical representation, which 
can be written as follows 

g(z) = -0+ [ P^dr(0 



T 

If a = da{l}, the mass at point 1, then we have formula (|122p with (l+t 2 )~ 1 da(t) = 
dr[(t - i)(t + i)- 1 ]. 

Next, our goal is to show that a — 0, and /? and da coincide with those from 
(|2T)) in section 3. 

From [7], p. 630, we have 

iF(iri) 
a = lim — : 

77— >+oo irj 

But (jll8p implies lim^^ +oc F(irj) — 1 so a = 0. In other words, measure dr has 
no mass at £ = 1. 

Notice now that each function P~ 1 (r, A)P*(r, A)/2 has positive real part as well 
and admits the same representation (|122|) with a r — 0, j3 r £ R, and absolutely 
continuous measure doy, which is the transplantation of some dr r . From (|76|) . 
we have a' r (X) = (27r) _1 |P*(r, A)|~ 2 . Measures da r are analogs of the so-called 
Bernstein-Szego approximations for OPUC. They converge weakly to da. Indeed, 
the convergence of P*~ 1 (r, A)P*(r, A)/2 to F(X) in C + implies convergence of the 
corresponding functions within the unit disc D. By Stone- Weierstrass theorem, that 
yields weak convergence of measures dr r to dr, (3 r to (5 and so the weak convergence 
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of da r to da. Since dr has no mass at £ = 1, the family of measures (1 + t 2 ) 1 da r {t) 
is tight, i.e. for any e > 0, there is T(e),R(e) > such that 

da r (t) 



1 + t 2 



< e 



(123) 



\t\>T(e) 



if r > i?(e). 

Let us show now that a coincides with the measure from (|27|) . Using ([9]) and 
we obtain 



+ / H r {x) exp(i\x)dx 



itx 
1 + t 2 



exp(-itx) 



da r (t) 



e-x.p(i\x)dx 



From (|10p . we get 

t 



+ / (t - s)H r {s)ds = i/3 r t + I ( 1 + _^ - 2 - exp(its) 



its 



da r (s) 



Fix t and take r oo in the last equation. The formula ([52")) and tightness p23p 
yield 



^ + / (t - s)H(s)ds = if3t+ I [ 1 + ^-^2 ~ ex P(^ s ) 



its 



da{s) 



Since the integral representation of Goo functions is unique (Theorem 1 1.2|) . we get 
the statement of the Theorem. □ 

We also have the following important 
Corollary 6.1. For any fix) G L 2 [0,p] ; the following identity is true 



f(x) exp{i\x)dx 



da(X) 



f(x) e~xjp(i\x)dx 



dX 



27r|P(p,A)|2 



(124) 



Proof. We have \P(r, X)\ — \P*{r, A)| for real A. Then, the Plancherel theorem for 
the Fourier integrals and (|53")l yield that the right hand side of (|124|) is equal to 



\f(x)\ 2 dx+ / H p (x-u)f(u)f(x)dudx 





p p 



\f{x)\ 2 dx + J J H{x-u)f(u)f{x)dudx (125) 



where (|92|) is used to get the last equality. Then, the formula (f50|) shows that the 
l.h.s. of (fT24l) is equal to r.h.s. of (p5|) . □ 
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The formula similar to (|117|) is true for Weyl-Titchmarsh function as well. From 
(|117j) and relation between F and f, we get 

F(x) = P* (P, A) - Hp, A) + F P (X) (P, {p, A) + P(p, A) ) 
lJ P4p,A)+P(p,A) + F p (A)(n( P ,A)-P(p,A)) 1 > 

Remarks and historical notes. The Weyl-Titchmarsh theory for Krein sys- 
tems was developed to some extent in [5Tj . 
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7. Dual system. Some simple considerations 

The dual Krein system is obtained by changing the sign of the coefficient A(r) 
(see Corollary 14. Due to Corollary 14. 11 functions P(r, A), P*(r, A) are continuous 
orthogonal polynomials for the dual system. They are usually called the dual 
continuous orthogonal polynomials. 

The dual Krein system can be characterized by the dual accelerant. Let us call 
it H. The relation between accelerant and dual accelerant is very simple. 

Lemma 7.1. For the dual accelerant H, we have 

X 

H(x) + H{x) + 2 J H(x - s)H(s)ds = 0, x £ R (127) 
o 

Proof. We have 



P,(r,A) (P*(r,\) 



1 



P*(r,A) \P*(r,X) 

Substitute (|90[) to the seconds factor and analogous formula to the first one gives 



which implies 

X 

H r {x) + H r (x)+2 J H r (x - t)H r (t)dt = 
o 

Then, use j92|) to get (fT27| . □ 

Clearly, the last Theorem allows one to find H from H by solving Volterra 
equation. The algebraic explanation to (|127[) is as follows. Consider two operators 



[Af]{x) = f{x) + 2 J H(x - u)f(u)du, [Af]{x) = f(x) + 2 J H(x - u)f{u)du 



acting in L 2 [0, r]. They both have positive real parts: 

ReA = I + "K r ,RcA = I + 'k r (128) 

and the formula (|127[) is equivalent to 

A = A~ X (129) 

These identities arise naturally from the solution to continuous Caratheodory- 
Toeplitz problem [45] . The relation between T r (x,y) and the dual resolvent kernel 
T r (x,y) is also quite simple and can be obtained from (|128|) and (|129[) . 

Consider the dual Weyl-Titchmarsh function F(X) — linv^oo P*(r, \)P~ 1 (r, A). 
Then, by Theorem E2 

F = F- 1 ,?=-f 

Now, let us study how the parameters of the Krein system change upon some 
simple transformations of the coefficient A{r). 
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Lemma 7.2 (Shift). Let A(r) be coefficient o/(57p. Then A® (r) = A(r) exp(irt),t E 
M. corresponds to 

o-W (A) = cr(A + t), H {t) (x) = exp(-itx)H (x), (x, y) = T r (x, y) exp(-it(x - y)) 

Proof. Let pairs {P, P*}, {P^\ P*^} be solutions of Krein system with coefficient 
A and A", respectively. Introduce Q = exp(— i\r)P, Q^> = exp(— i\r)P^ . We 
have (see formula ([50)1 ): 

Q'= -iexp(-£Ar)P*, Q(0,A) = 1 
Pi= -Aexp(iAr)Q, P*(0,A) = 1 

and 

f QW'= -^exp(-i(A + i)r)Pi t) , QW(0,A) = 1 
\ Pi* V = -Aexp(t(A + t)r)QW, pJ 4) (0,A) = l 

Therefore, Pi 4) (r,A) = P*(r,A + t). By TheoremEl (2n)- 1 \PP(r, X)\~ 2 dX 
da^\ we get the shift in the measure. For dual system, we have the same results. 
Then, by j90j), 



(130) 
(131) 



OO OO 

1 + 2 J H^ t] (s) exp(isX)ds = 1 + 2 J H r (s) exp(is(A + t))ds, (132) 

o o 

and we have the needed formula for ffw. The way resolvent kernel changes is easy 
to obtain from (O or i35i □ 



P, (t) (r, A) 
P* (t) (r, A) 



It is an easy exercise to show directly that if H(x) is an accelerant, then H{x) exp(— itx) 
is an accelerant as well. Coefficient /?w from formula (|2T|) changes in a more intri- 
cate way. It can be recovered by noticing that (A) = P(A + t) and by integral 
representations for both Weyl-Titchmarsh functions. We then get 



— oo 

Lemma [721 has its analog in the OPUC theory, rather than the following Lemma 
Lemma 7.3 (Dilation). For any 7 > 0, coefficient jA(-yr) corresponds to 

0"(7)( A ) = 70-(7 -1 A),7 > 0,P( 7 )(x) = jH(jx),T {lhr/ ^(x,y) = jT r (jx,jy) 
Proof. Under the change of variables p — yr, system ([57)1 changes as follows 

( dP(yr, A) = , a) - 7 A( 7 r)P* ( 7 r, A), P(0,A) = 1, 

■> \ (133) 
^ = -7^(7r)P(7r,A), P*(0, A) = 1 

which proves the Lemma following the same lines as in the proof of Lemma l7.21 □ 

Again, if H(x) is an accelerant, then one directly checks that 7^(72;) is also an 
accelerant. Since P( 7 )(A) = P(A/ 7 ), we have 
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Lemma 7.4 (Conjugation). The coefficient A(r) corresponds to the measure a : 
er(7) = <j(— J) for any Borel set I, accelerant H(x), resolvent kernel T r (x,y), and 
coefficient —(3. 



Proof. Take conjugation of ([57]) . The pair {P(r, — A), P*(r 7 — A)} solves Krein sys- 
tem with parameter A and coefficient A{r). New Weyl-Titchmarsh function is equal 
to F(— A). The integral representation for F(X) yields the value of (3. □ 

Consider the case when A(r) is real. Then, a = a and H(x) is also real. It 
is continuous on the whole line provided that A(r) S C[0, oo). In (|122[) . constant 
(3 = 0. Later on, we will consider this case in greater details. 

The next calculation will be important to understand the scattering problem for 
Krein system and Dirac operators. Consider Krein system on the interval [0, R] 
with coefficient A^ R \r) = A(R - r) for r E [0, R}. For r > R, we let A^fr) = 0. 

Lemma 7.5 (Mirror symmetry). The Schur function of Krein system with the 
efficient A^ R ' (r) is equal to 



coe ,„ 

/C*>(A) = 



2t*(iJ,-A) 



Proof. Consider the matrix AT that solves Krein system 

X' = VX, X(0, A) = / 

and V is given bv ([66]) . At the same time, matrix Y(r, A) = JX(R—r } —\)X~ 1 (R,—\)J 
solves the following system 

Y' = V(R — r)Y, F(0) = I 
Therefore, new Wall polynomials are 

2t (fl) (i?,A) = exp(Ui?)2U(i?,-A) 

& R \R,\) = exp(iA-R) —A) 
That finishes the proof. □ 

Remarks and historical notes. The dual systems were studied before, e.g. 
|46j . They also appear in the solution to various continuation problems and in 
continuous Caratheodory-Toeplitz problem. 
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8. SZEGO DISTANCE AND KREIN SYSTEMS 

As any orthogonal system, {P(r, A)} has the reproducing kernel. Consider the 
scale S p = 3 3 [ 0iP iI/ 2 (K) of Paley- Wiener spaces. Recall that for any p > 0, this space 

consists of functions /(A) that can be represented as 

p 

/(A) = J f(x)exp(iXx)dx 
o 

with f(x) e L 2 [0,p}. 

Lemma 8.1. The following function K„{X' ', A): 



.P,(p, X)P*(p,X>)-P(p,X)P(p, A') 



K p (X', A) = / Pfo A')P(x, A)da: = f VP ' ; ^_ -, ; ^' ' (134) 
o 

is the reproducing kernel in S p space, i.e. 

f(X') = (f(X),K p (X',X)), A'GC (135) 
where the inner product is defined as follows 

oo oo 

dA 



(/i,/ 2 )= y /i(A)/ a (A)«fer(A) = J /i(A)/ 2 (A) — — - — (136) 

— oo — oo 

Proof. The second equality in (|134[) is formula ([55)) . For any fixed A', P(x, A') G 
L 2 [0,p], so the kernel if p (A',A) itself is an element of S„. Formula (|135[) follows 
from Plancherel-type identity (|4"T|) . Indeed, by (|4"5)) . we have 

p 

/(A') = f f 1 (x)P(x,X)dx 



with some /i(a;) G L 2 [0,p]. Now, using (|47|) . we get 

p 

(f(X),K p (X',X))= f h{x)P{x,X')dx = J{X') 



The second equality in (|136|) is the contents of Corollary [6Tj □ 
The reproducing kernel property (|135[) yields 

K P {X 1 ,X 2 ) = (K p (X u X),K p {X2,X)) 
Together with Cauchy inequality, that implies 

7(A')f < \\f\\l H M, d a)K P (X',X') (137) 

and the equality holds if and only if /(A) = jK^X' , A), I7I = 1. 

Lemma 8.2. The following identity is true 

00 

1 



K P {X',X') / e 5 p n /(A')=i 



|/(A)|^a(A)=m 2 (A') (138) 



for any A' G C. The minimizer is unique and is given by /(A) = K (A', X')K p {X' , A). 
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Proof. Divide (JT33 byK p (X',X') /(A') 
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Now, the natural question to ask is what happens if p — > oo? Since S P1 C S P2 
for p\ < p2, the minimum m p (X!) decreases. Now, can we characterize the case 
when it decreases to zero? To do that, we need the following classical result (see, 
e.g. [21] . page 84). We give its proof in Appendix (Theorem 1 17. 2 j) . 

Assume that da is a positive measure on R such that 



da(X) 
1 + X 2 



< oo 



Consider the linear manifold X of functions /(A), having the following representa- 
tion 

/(A) = J exp(iXx)f(x)dx, < r*i < r 2 

ri 

where f(x) € C^ri,^] and is zero outside [ri,ra] Q [0, oo). Notice that each 
/(A) £ L 2 (da). Denote the closure of X in L 2 (da) by X. 

Theorem 8.1. The linear manifold X is not dense in L 2 (da) if and only if 

oo 

lno-'(A) 



1 + A 2 

Moreover, the following formula is always true 



dX > — oo 



(139) 



Dist 



1 



A-A ' 



X 



1 



\/2ImAo 



exp 



ImAo f ln(27rcr'(A)) 
J |A-A | 2 



dX 



, A e 



(140) 



Next, we will apply this Theorem to the Krein systems. Let da be the measure 
generated by some Krein system (the measure from (|27p). Recall the definition of 
S r and notice that for each finite r, S r C L 2 (da) (see ([5T])). Denote the closure of 
U r> o<SV in L 2 {da) by 5. 



Lemma 8.3. //dcr is generated by some Krein system, then X = S. 

Proof. It is clear that X C S. On the other hand, any function f E S r can be 
approximated in L 2 (da) by a sequence of functions from X. That easily follows 



from (I5T 



Lemma 8.4. The following formula 
'2ImA 



Dist 



,S 



L 2 (da) 



inf m r (Xo) = m oo (A ) 

r>0 



A — Ao 
is true /or any Ao G C + . 

Proof. Denote the l.h.s. by I\ and the r.h.s. by I-i. We have 



□ 



(141) 



2 Im A 

A-A 



/(A) 



2i Im Ao A — Ao 



A-Ao 



A- A 
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and so 



Since 



j 1=Dtot ('2*^0,^5, (142 ) 



A — Ao A — Ao 



L 2 (da) 



oo 

ot \ f r-\ ^ c \ ^ 2iImA A-A 2ilmA 

21mAo / exp(iAa:) exp(— iAo^jaa; = — = — , =— = 1 = — (143) 

J A — Ao A — Ao A — Ao 

o 

any function 

?^_*Z*o/ (A ),/ (A)e S 
A — Ao A — Ao 

can be approximated in L 2 (da) by the sequence of functions f rn ( A) S SV„ , / rri (Ao ) = 1. 
Therefore, /i > 1-2.. Assume that /(A)- arbitrary function from some S r and 
/(Ao) = 1. Then 

/(A) ^ ^ - ^(A) (144) 
A — Ao A — Ao 

with g(X) G S. Indeed, the function 

/(A) - 2am _ A ° 
A-A 

belongs to H 2 (<C + ) and has zero at A = Ao. So, by Paley- Wiener Theorem, (I144p 
holds with g G H 2 (C + ). Let us show now that g E S. The first formula in (I143|) 
suggests 



. , , s A - A 

gix) = 



-/(A) - 2 , Im ^° — =-rr / exp(ix(A - A ))cte 

1 - cxp(in(A - A )) 



r n (A) 



The Paley- Wiener Theorem yields that the first term belongs to S p where p = 
max(r, n). One can easily see that ||r n ||2,cr — > as n — > oo. Thus, the formula (|144|) 
holds. Due to fl42]) . Ji < J 2 and therefore h = h- □ 

The next result describes the continuous analog of the Szego case in OPUC 
theory. If any of the conditions bellow is satisfied, we will say that da G (Szego). 

Theorem 8.2. (The Szego case) The following statements are equivalent 

(a) The operator from the Theorem \3.2\ is not unitary. 

(b) Inequality 

oo 

W(A) 



1 + A 2 



■dX > -oo (145) 



holds. 

(c) 



sup |/(A )| = (Ao) < oo 

r>0,/(A)eS r ,||/(A)|| 2 , <T = l 

for at least one (and then for all) Ao G C + . 

(d) P(r, A ) G L 2 (R + ) for at least one (and then for all) A G C+. 

(e) liminff— >oc |P*(r, Ao)| < oo for at least one (and then for all) Xq G C + . 
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Proof, (a) and (b) are equivalent. Indeed, by Lemma 13.21 the range of coin- 
cides with S. From Lemma 18. 3[ we have X — S and then we only need to use 
Theorem HU 

The formula (|59|) with A = \i = Ao shows that if (d) holds at some Ao then (e) is 
true as well at the same point. The converse is also true. 

The identity (|138|) and the formula for reproducing kernel yield 

oo 

m" 2 (A )= sup |/(Ao)| 2 = ifoc(Ao,Ao)= f \P(x,\ )\ 2 dx 

r>0,/(A)eSV,||/(A)|| 2 ,„=l J 

and that proves equivalence of (c) and (d) for fixed Ao- 

Let us show that (c) is satisfied with some Ao if and only if (a) holds. Assume 
that the operator is unitary. That means its range is the whole L 2 (da) and 
(A — Ao) -1 G S, for any Ao G C + . Due to Lemma l8~4l (c) fails. Conversely, assume 
(c) fails for some A € C+. Then, by LemmaH (A - A ) _1 G S. The Theorem 
now implies 

oo 

In er'(A) 

d\ — —oo 



|A-Ao| 2 

— oo 

but that means (b) fails and therefore (a) fails too. Notice that both (a) and (b) 
do not depend on parameter Ao- Therefore, if any of (c), (d), or (e) holds for some 
A then it holds for all A G C+. □ 

It is not in general true that |P*(r, A)|,A G C + is even bounded as r — > oo 
under the conditions of the Theorem 18.21 That is due to continuous nature of the 
problem. Moreover, it is possible that |P*(r, A)| has a limit, but lim^oo P* (r, A) 
does not exist inspite of the fact that the corresponding A(r) — ► at infinity and 
A(r) G i p (M + ) for any p > 2. This phenomena was observed for the first time by 
Teplyaev (see [7H[7S]). That can be explained as follows: in the discrete case, the 
orthonormal polynomials are usually normalized such that they have the positive 
leading coefficient. For Krein systems, normalization is quite different: P*(r, A) are 
normalized to be equal to 1 at infinity, the point on the boundary of C + . Therefore, 
the argument of P* (r, A) is not stabilized and that leads to the ambiguity in the 
definition of linv^oo P» (r, A). 

Consider some Ao G C + . If conditions in Theorem 18. 21 are satisfied, then there is 
a sequence r n G [n, n + 1] — * oo such that P(r n , Ao) — * 0. Take the outer function 



n(A) = -^exp 

V Z7T 



3C 



1 f (1 + sA) In a'(s) , 
-ds 



2-Ki J (A-s)(l + s 2 ) 



(146) 



that satisfies |II(A)| = [2tkt' {X)]- 1 ' 2 for a.e. A G K. Notice that [(A + i)II(A)] 1 G 
P 2 (C + ) and is outer. 

Lemma 8.5. If da G(Szego) and r n — > oo is such that P[r n , Ao) — > for some 
Ao G C + . then the following convergence \P*(r n , X)\ — > |n(A) | takes place uniformly 
inC+. 

Proof. From the Theorem l6.2l we know that the sequence h n (X) = [(A + i)P» (r n ,X)\ 
is bounded in H 2 (C + ), i.e. ||/i„(A)||i2( ffi ) is bounded. Assume that h(X) is any 



-l 
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L 2 (R)-weak limit point of this sequence. Then, h(X) G H 2 (C + ) and the conver- 
gence is uniform in C + over the same subsequence 7ifc. From (I76p . 

Re [P,- 1 (r nk ,X)P4r nk ,X)] > |P*(r„ fe , A)|" 2 (147) 
Taking k — > oo, we get 

| (A + i)| 2 |/i(A)| 2 < ReP(A) 
for A G C + , where P(A) is the Weyl-Titchmarsh function. From (|119[) . we get 

(A 2 + l)\h{\)\ 2 < 2na'(X) 

for a.e. Ael. Therefore, from the multiplicative representation of (A + i)h(X) 

|(A + i)ft(A)| < |n- 1 (A)|,AeC+ 

At the same time, at Ao, we have \(Xo + i)h(Xo)\ = |Il _1 (Ao)|. It follows from (|134p . 
(fI55j) . (|T30|) . |141|), ([TiH]) . Therefore, (A + i)h(X) is an outer function different 
from II -1 (A) only by a unimodular constant factor. Thus, for any subsequence rife, 
\P*(r nh , A) | — > |II(A)|. That means we actually have convergence over the whole 
r n . □ 

It is known that all outer functions II 7 (A), satisfying |II 7 (A)|~ 2 = 27rcr'(A) a.e 
on R, have the following representation 



n 7( A ) = ~^ exp 



3C 



1 f (1 + sA) In <t'(s) , 



, 7 e[0,27r) (148) 



i.e. they can be parameterized by the angle 7. The function 11(A) = IIo(A) satisfies 
the following normalization condition: Tl(i) > 0. There are some quite interesting 
examples |75] when P* (r n , A) — > II 7 (A) and the constant 7 depends on the choice 
of subsequence r n . In the meantime, the following is true [64] 

Lemma 8.6. Assume that A(r) is real-valued and r n is such that P(r n ,Xo) — > 
for at least some Xq G C + . Then, P*(r n , A) — > 11(A) uniformly in X G C + . 

Proof. Following the proof of the previous Lemma, we get convergence of (r nk , A) 
to some II 7 (A) uniformly for A G C + . Taking A = i, we get 7 = 0. Indeed, 
P*(0,i) = 1, P*(r, i) is real and has no zeroes for r > 0. Therefore, it must be 
positive for all r > 0. Thus II 7 (i) > and 7 = 0. □ 

In the OPUC theory, we can not directly characterize the set of moments such 
that the corresponding measure belongs to Szego class. The same is true for the 
Krein systems: we are not aware of the characterization of the Szego case in terms 
of accelerant. In the meantime, 

Lemma 8.7. If da G (Szego), then for the dual system we also have da G (Szego). 
Proof. We know that for any Krein system, 

PM) F(A) 

as r — * 00 and F(X) has positive real part in C + . Therefore, if condition (e) of the 
Theorem 18.21 is satisfied for the original Krein system, it must be satisfied for the 
dual one as well. □ 
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Next, let us show that we have weighted P 2 -convergence for P*(r, A) for A€l. 
First, we need the following auxiliary result 

Lemma 8.8. Assume that the Szego case holds. Let r n be a sequence such that 
P(r n ,i) — > 0. Then 

J \P(r n ,\)\\ 1 
rSL ./ A 2 + l ^ = 2 



Proof. Indeed, we have 

K r (i, A) = i 



P* (r, A)P* (r, i) - P(r, A)P(r, i) 



A + i 



or 



P*{r,i) = i 



P(r,A) 



P(r,t)+#r(t,A) 



X + i v ' ' A+i 
From Lemma 18751 we know that |P*(r n ,i)| — > |II(i)|. Then, since P(r n ,i) — > 



(149) 



lim 



P*(^n,A) 



A + i 



2, do 



|n(z)|- 2 lim |i^ rn (i,i)| = i 



□ 



The following result establishes an L 2 (da, R) asymptotics of P*(r, A). It will be 
used later to prove existence of wave operators for Dirac equation. 

Lemma 8.9. Assume that da G(Szego) and r n —> oo is such that P{r n , i) — ► and 
P*(r„,A) -> n 7 (A) /or A G C+, ft € [0,2rr)J. Then, 



A 2 + 1 



P*(r„,A) 



n 7 (A) 



- 1 



dX -> 



(150) 



as r„ ^ oo. 

Proof. The left-hand side of (|150[) is equal to 



1 



A 2 + 1 



P*(r„,A) 



n 7 (A) 



dX- 



dX 



A 2 + l 



2 Re 



P.(r„,A) 

(A 2 + i)n 7 (A) 



By the Cauchy formula, 



P*(r n ,A) 7rP*(r„,i) 

uA = — ► 77 



(A 2 + i)n 7 (A)"" n 7 (i) 

— oo 

as r n — * cxd. The first term of (|151p can be written as 

oo oo 

1 



dX (151) 



(152) 



2tt 



' \P*(r n ,X)\ 2 do-(X) -2tt 



A 2 + 1 



A 2 + l 



|P i ,(r n ,A)| 2 do- s (A), 



where d<r s (A) is the singular component of cicr(A). From the Lemma 18.81 we infer 

oo 

^J—\P*(r,X)\ 2 da(\)^~ (153) 
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Bearing in mind (Pol]) . (fT52]) . and (fT53]) . we have (fT50|) . □ 
Remark. We also proved 

oo 

/ 7I^yl^(^A)| 2 da s (A)->0 (154) 

— oo 

as v n — > oo . 

We want to finish this section with the following observation that relates the 
regularity of 1 11(A) | at +ioo to some approximation problem. Consider the function 
/ Ao (A) = (2ImA ) 1 / 2 (A- A )- 1 ,AeM. Since 

o 

(2ImA ) 1/2 (A- Ac,)" 1 = i (2ImA ) 1 / 2 / e- tX ° x e tXx dx 



the function f\ (A) has frequency concentrating near zero as Im Ao — > +oo and the 
constant L 2 norm. How regular the distance Dist(/A (A), S) behaves as ImAo — > 
+oo depends on the regularity of II (A) at infinity. Indeed, from Lemma 18.41 

Dist(/ Ao (A),5) = |n(A )r 1 

That infinitesimal phenomena is not present in the discrete case. 

Remarks and historical notes. The approximation results of this section 
can be interpreted in the framework of prediction theory for stationary Gaussian 
processes with continuous time |36j . The original paper by Krein contained some 
inaccuracies in the formulation of the Theorem 18 . 21 and the same mistake was made 
in some later papers. The correct statement was given later by Teplyaev [74\ I75j. 
Some sufficient conditions for the Szego case were given in the series of papers 
[541, H~5] , There is no known criteria in terms of A(r) for the Szego case to hold. In 
the meantime, if one assumes some regularity of A(r), say, A(r) <G i°°(M + ) then 
da € (Szego) if and only if A(r) € H~ 1 (M. + ) (see [15]). It is probably impossible 
to give reasonable characterization of the Szego case in terms of A(r) without 
any apriori assumptions. For example, one can construct a sequence of compactly 
supported iW(r) with growing H~ 1 (R + ) norms but such that the corresponding 

sequence (oo,i) is bounded. This can be achieved by a simple modification of 
Teplyaev's example [75] . 
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9. SCHUR'S ALGORITHM AND APPROXIMATION OF CONTINUOUS ORTHOGONAL 

SYSTEM BY DISCRETE ONES 

It is well-known [53], that any function f(z) 6 B(H>) can be expanded into 
the continued fraction (the so-called Schur's algorithm). This expansion can be 
obtained by the iteration of 

= zf+^+a = = 

1 + fl„z/ n+ i(z) 

By doing so, we obtain the one-to-one correspondence between B(JD>) and all se- 
quences {a n } such that |a„| < 1. The Geronimus theorem asserts that these 
so-called Schur parameters are actually equal to Verblunsky parameters for the 
measure dr in the representation 

\^f\ = I I^MO (156) 

In this section we will study Schur's function associated to Krein system. The 
Schur function f(A) associated to the Krein system was introduced in the Theo- 
rem HHJ Functions f p (A) in the representation 1|117|) are Schur's functions that 
correspond to the same Krein systems but on the interval [/?, oo). This is the same 
as if we would take A p (r) = A(r + p) with r G K + . 

Lemma 9.1. For any fixed A G C + , the Schur functions f r (A) are continuously 
differentiable in r and satisfy the following equation 

^M = - iX f r (X)+A(r)-A(r)f r (X) (157) 

Proof. The smoothness of f,.(A) in r follows immediately from (| 1 1 7[) . Taking deriv- 
ative of (I117P in r at r = 0, we get 

df r (X) 



= -iAf (A)+A(0)-A(0)tf(A) (158) 

dr r=o 

Now, (|157|) follows from the definition of f r (A). □ 

It is very important to keep in mind that the initial condition for (|157[) is f (A) = 
f(A) and it is not independent of the coefficient A{r). Now, let us compare the 
continuous and discrete Schur algorithms. Consider the following Mobius transform 

T 7 {Z) = — — 
l+JZ 

The inverse to r 7 (z) is equal to r_ 7 (z). 

For any 7 G B, r 7 is a conformal map of ID onto D that takes to 7. Another 
important property of 17 is the preservation of the pseudohyperbolic distance on 



P(Z1, z 2) 



zi ~ z 2 



,p(t 7 (zi),t 7 (z 2 )) = p(z 1 ,z 2 ) (159) 



1 - Z\Z 2 

Given any function f(z) from B(U>) the Schur algorithm can also be defined as 
follows. The n-th Schur's iterate f n {z) is defined by the relation 

f(z) = S , Z)O0 ,... )O „_ 1 (/„) = T ao o ZT ai o ... ZT an _ t o (zf n ) (160) 

and 

fn{z) = S~l o ^ ani (/) = z" 1 r_ Qii _ 1 o z" 1 T_ a „_ 2 o ... z-V_ O0 o / (161) 
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Notice that for z € T the map S z ,a ,...,a n _ 1 is the composition of rotations and 
Mobius transforms. 

The differential equation (|157|) is a continuous analog of (|155p . It is Riccati 
equation and the Cauchy problem for solving it from the right to the left happens 
to be well-posed for suitable initial data: 

Lemma 9.2. For any A(r) € L l [Q, R] and any f(0) = f € B, A e C + , there is the 
unique solution f r (A) to Cauchy problem for equation ftl57\) with initial condition 

Proof. Consider Y(r, A) = 2/2(7*)) , solution to the following Cauchy problem 

y[ = -i\y x + A{r)y 2 , y 1 (R,X) = f Q , 



y' 2 = A(r) yi , y 2 (R,\) = l 

and solve it for r £ [0,R]. Simple calculations show that 

Y(r, A) = X B(r) (0, R-r, X)Y(R) = X\r, R, X)Y(R) (162) 

where Y(R) = (/o, 1)* and Xb is the transfer matrix for the Krein system with 
coefficient B(r) = A(R — r). We have A e C+ so Xb is J-contraction by the 
Theorem 14. II Thus, we have |j/i(r, A)| < 1 2/2 A) | . In particular, y2{r,X) ^ since 
otherwise Y = on [0, R\. Consider f(r, A) = yi(r, A)?/ 2 _1 (r, A). The straightforward 
calculation shows that f(r, A) is solution to our Cauchy problem and uniqueness 
follows from the general theory of ODE. □ 

In analogy with discrete case, we denote the solution of this Cauchy problem at 
zero by Sx,a,rQo) and now we have f(A) = S\,A,r(f r (A)), the direct analog of (|160p . 
The formula (|162[) shows that S\,A,r allows the following representation 

St,(r,A)z + <B(r,A) 

£>\,aA z ) = 7^7 — n — r~™i — vT ^ lb6 > 
<B*(r, \)z + 2l(r, A) 

Notice that the inverse to S in discrete case is not contraction anymore and we 
have the same problem in the continuous setting. 

As we mentioned earlier, the class of Schur functions f(A) in the continuous case 
can not be all B(C + ). For example, (|118[) must hold. Next, we will describe the 
subclass of B(C + ) in which A(r) have the meaning of intrinsic parameters of the 
function /(A) just like {a n } are intrinsic parameters of f(z) € £>(D). 

Let a function C(x) S L^ oc (K + ) be given. For any R > 0, consider the operator 
Qr acting in L 2 [0, R] by the following formula 

X 

G R f(x) = [C(x- u)f{u)du (164) 



We start with the definition. 

Definition 9.1. The function s(A) € S(C + ) if the following is true: 

(1) There is a function C(x) £ L^ oc (R + ) such that for any R > there is a 
function $r(A) G H°°(C + ): 

R 

s(A) = / C(a;)exp(iAa;)da; + exp(iA J R)$ jR (A) (165) 
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(2) For any R > 0, 

WGrWlho.h] < 1 (166) 

It is an easy exercise to see that C(x) and $i?(A) are both uniquely defined 
for any s(A) G S'(C + ). Conversely, if the function C(x) is given, then there is at 
most one s(A) G S(C + ) having C(x) as a function in the formula (|165p . Indeed, 
assume that there are two s^'(A), s^(X) having the same C(x) in (|165[) . Then, 
(s^(X) - s^(X))/(\ + i) G H 2 (C+). At the same time, 

..■»(A) .^(A) = <(A) <(A) 

X + i v ; A+z 

for any R > 0. Since (^(A) - $^(A))/(A + i) G H 2 (C+) as well, we have 

- s (i)( A )- s (2)(A)- 



0> 



[0,-R] 



= 



A + i 

for any R > 0. So, a W(A) = s (2) (A). 

Consider s(A) G S'(C + ). We have s(A) € H°°(C+). In the space H 2 (C+), denote 
the operator of multiplication by this function by S. Also, consider the operator C 
acting on L 2 [0, oo) and given by the formula 

X 

e/(ar) = J C{x- u)f{u)du (167) 
o 

Since Iljo^CIIjo.ii;] = 6^ and ||Ci?|| < 1, operator C is well-defined on I/ 2 [0,oo) and 
is contraction, i.e. ||C||2,2 < 1- 

Lemma 9.3. One has the following inclusion: S(C + ) C B(<C + ). The operators C 
and § are unitary equivalent. 

Proof. Indeed, s(A) G iJ°°(C + ) and it is known that 

IN| ff °o (c+) = sup \\sg\\ H 2 {c+) = ||S|| 

Il9llff2(c+) =1 

The condition ()166[) is equivalent to the estimate ||3 3 [o,.r]S3 3 [o,.r] I < 1- Then, ||S|| < 
1, or s(A) G B(C + ). The unitary equivalence of § and C follows from the unitary 
equivalence of operators Il^fljCII^ij] and J'[o,_R]§J'[o,fl] via the Fourier transform. 

□ 

Recall the definition of the accelerant: given Hcrmitian H(x) E i 2 oc (IR), we say 
that it is an accelerant if the operator 

/ + Vi R > 

for any R > and IKr is given by ([2"5]) . Given any C(x) G L 2 oc (IR + ), consider the 
function H{x) G L 2 oc (M + ) which is the solution to 



H[x) + C{x) + J C(x - u)H{u)du = (168) 
b 

The direct iteration of the equation proves existence and uniqueness of this H(x). 
Let H(-x) = H(x) for x > 0. 
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Lemma 9.4. The function H(x) is an accelerant if and only if C(x) is such that 
U66\) holds for any R > 0. 



Proof. For any R > 0, consider the Caley transform of Cr: 

i + Ur = (i - e R )(i + e R )~ l (169) 

Since Cr is a Volterra operator, the Caley transform does exist. Moreover, Ur is a 
Volterra operator with the kernel given exactly by H: 

X 

U R f(x) = 2 J H(x - s)f(s)ds (170) 
o 

This is an easy corollary from (|168p . Clearly, Re(I + Ur) = I + "Kr- Now, the 
equivalence of / + "Kr > and (| 1 66[1 is a simple algebraic fact. □ 

Now, we can easily characterize the class of all C(x) that generate s(A) G S(C + ). 

Theorem 9.1. The Schur functions of Krein systems with A(r) G Lj 2 oc (R + ) are in 
one-to-one correspondence with functions s(A) G S(C + ). For each s(A) G S(C + ), 
the coefficient A{r) of the associated Krein system plays the role of the Schur pa- 
rameter. 

Proof. Assume s(A) G S(C + ) and C(x) is the corresponding function. Denote by 
H(x) the accelerant corresponding to C(x), i.e. 



H{x) + C(x) + J C(x - u)H{u)du = (171) 
o 

Then, by Lemma \9. 41 H(x) is an accelerant that generates the Krein system with 
coefficient A(x) G L 2 oc (R + ). Consider the corresponding Schur function f(A). Let 
us show that f(A) = s(A). For each R > 0, we have (|1 1 7[) 

g(fl,A) + f fl (A)g > (fi,A) 
H ; a(iZ,A) + f B (A)B.(i2,A) 1 } 

for any R. Clearly, by (|114|) with r% — oo, r-i = R 

_ «B(fl,A) exp(»Ai?) f R (X) 

K j 2l(i?,A) 2l(fl ) A)(2l(i?,A) + f fl (A)Q5 <f (i?,A)) 1 ' 

By pigjl . the right-hand side of ([175)) is equal to exp(iXR)$ R (X) with <J>r(A) G 
H°°(C + ). Consider the function <8(R, A) ST^i?, A). Due to Levy- Wiener Theorem, 

oo R 

Q5(-R A) f f 

' = / Cr(x) exp(iXx)dx = / Cr(x) exp(iXx)dx+ 
w(R, X) J J 





+ exp(iAi?) J C R (x + R)exp(iXx)dx, C R (x) G L 1 (R+)nL 2 ( 
o 

and the last term can be written as exp(iAi?)$i?(A) with < I ) r(A) G H°°(C + ). We 
can also write 

2l(i?,A) \ P,(R,X) i I P*(i?,A) i 
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From PJJ), ([92]) and (fT7T|) . we infer C R {x) = C(x) for x £ [0,i?]. Since R is 
arbitrary positive, f(A) G S'(C + ) and f(A) = s(A). 

Now, assume that the Krein system with the coefficient A(r) is given. Repeating 
the arguments above, one has f(A) G S(C + ). □ 

Remark 9.1. It follows from (I168j) that C(x) on the interval [0, R] depends only on 
the values of A(r) on the same interval. This is because an accelerant has analogous 
property. 

Notice that s(A) £ S(<C + ) implies certain regularity at infinity, a boundary point 
of C + . For instance, s(iy) — > as y — > +oo. So, there are plenty of functions in 
B(C + ) that do not belong to S(C + ). 

The same arguments immediately yield 

Remark 9.2. A(r) G C[0,oo) iff H(x) G C[0,oo) iff C(x) G C[0,oo) in the corre- 
sponding representation for f(A). 

Consider this case for the rest of the section. Let C(x,r) be function associated 
to / r (A) by formula (|165[) . 

Lemma 9.5. (Continuous analog of Geronimus theorem). The following relation 
holds true: C(0,r) = —A(r). 

Proof. It is enough to prove the statement for r = 0. From (|171| . we have C(0) = 
-H(+0). Then, by Lemma \5M C(0) = -A(0). □ 

Notice that the value of C(r) at zero gives the main term of asymptotics of f (A) as 
A = iy,y — > +oo. For instance, if A(r) G C^O, oo), then H(x),C(x,r) G C^O, oo) 
as well (see (|84|) and ()171|) ) and Lemma [9751 yields f(r,iy) = — A(r)/y + d{y~ 1 ). In 
general, for A(r) £ C[0,oo), we have asymptotics in the mean (see Lemma 117.51 in 
Appendix) : 



lim — / sf(r, is)ds = — A(r) 

y^oc y J 





Anyway, the number C(0,r) = — A(r) is an intrinsic parameter of the function 
f(r, A). Therefore, Lemma 19.51 can be regarded as the continuous analog of the 
celebrated Geronimus theorem which says that the Schur parameters of the function 
from B(H) coincide with the Verblunsky parameters of the associated sequence of 
orthogonal polynomials. 

Lemma 9.6. Assume A(r) G C^O, oo), then C{x,r) is continuously differ entiable 
in x and r and satisfies the following nonlinear integro- differential equation 

dC(x,r) _ dC(x,r) 
dr dx 



A(r) C(x-u,r)C(u,r)du (174) 



o 

Proof. Let us prove smoothness in fir — {r > 0, x £ [0, T]} for any T > 0. 
Apply (I84|) to the Krein system on [r, oo). We have the corresponding r-dependent 
function A^ (T, x) £ C 1 ^). Notice that A^(T,x) =rf ) (T-x,0) and 

T 

T { t \x 1 0)+ f H^(x-u)T^ ) (u,0)du = H^(x) 
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To prove (I174|) for, say, < r < R and < x < T, we can consider new 
A\ equal to A on [0, R + T], smooth on K + with compact support. Then, new 
d(x,r) = C(x,r) for < r < R,0 < x < T. Moreover, &(x,r) G L 1 (R+) in x 
and the formula (|165|) holds with R = oo. Then, if one substitutes (|165[) to (|157jl . 
equation (|1T4[) pops up. □ 

Plug A(r) = — C(0,r) into this equation and solve the first order PDE with 
boundary condition C(x,0) regarded as known. Then (I174[) becomes a nonlinear 
integral equation which one tries to solve by iterations. Since A(r) — — C(0,r), 
that gives us a solution to inverse problem since C(x, 0) can be read off the spectral 
data, say da. 

Next, let us focus on the differential equations (|157[) and (|57|) . Looking at the 
formula (|160[) . one might guess that the map S\,A,r should also be represented as a 
combination of Mobius transforms and certain multiplications. This is indeed the 
case. The following result gives approximation of continuous orthogonal polynomi- 
als by the sequence of properly scaled polynomials orthogonal on the unit circle. 
These discrete polynomials are given in terms of Schur parameters that depend 
upon the step of discretization h. 

Theorem 9.2. Let A(r) € C[0, oo). Fix any r > and consider the sequence of 
Verblunsky coefficients 



<$> = hA(t 1 ),a\ n ' = hA(t 2 ), = hA(t n ) 



Ah) 



(h) 



and 



= 0, j > n, (h = r/n,tj = jh) 

where h is chosen so small that all of these coefficients are less than one in absolute 
value. Consider the discrete transfer matrix generated by these coefficients 

M (0, k, z) = W{a k )Z . . . W(a Q )Z 

with 



(175) 



Z(z) = 



z 
1 



,W(af) = 



1 

— at 



1 



Then, we have 



X(r, A) = lim M(0, n, exp(i\h)) 

h—>0 

and the convergence is uniform over A in compacts in C. 



(176) 



Proof. Consider small h > 0. Then, from the definition and properties of the 
multiplicative integral [llj . we have 



X(r, A) = / exp[V{t)dt] 



lim 



(1 + hV(t n )) . . . (1 + hV(h)) 



(177) 



lim 



W(a^l 1 )Z(w)W(a^l 2 )Z(w) . . . W(a^)Z(w) 



(178) 



where w = exp(i\h). Now, the statement of the Theorem is an elementary corollary 
from ([TTg)) . □ 



The next Corollary follows directly from the Theorem. 
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Corollary 9.1. For continuous orthogonal polynomials, 



P(r, A) P(r, A) 
P*(r,A) P*(r,A) 



lim 



P n ( W ) P n (w) 



where w = exp(iA/i), polynomials Pk,P£ are monic orthogonal polynomials gen- 
erated by the prescribed Verblunsky parameters and Pj~ , P£ are ditaZ to them. The 
convergence is uniform in A /rom any compact in C. 

Corollary 9.2. Let ^4(r) 6 C[0, oo). For the map S\ t A,r{z), we have 

Sj r i, r (z) = limS („) (h ) (z), (179) 

h^O ' '■■■> a n-l 

w — exp(i\h) and the convergence is again uniform over the compacts in C. 

Proof. (|179jl follows from the formula (|163p . Indeed, in discrete setting, there is a 
formula analogous to (|163p (see [40], formula (4.19)) 

A n (z) + zB*(z)f 



>,.*(/) = 



B n (z) + zA* n (z)f 



We also have 



21* (r, A) 55* (r, A) 
58(r,A) a(r,A) 



= lim 



-«;A„(iu) 



(180) 



and j4„ , P„ are the standard Wall polynomial^. On the other hand, by Lemma [4T3l 
JXJ is the transfer matrix for Krein systems with coefficient —A. Comparing the 
corresponding formulas to (|163[) . we get the statement of the Corollary □ 

Remark 9.3. One might wonder why the formula (|179p contains the sign minus in 
front of Al The answer to this question is contained in the definition of continuous 
Schur function and map S\ } a,t- Indeed, we defined f as 

X 2 i(r,X) 



f(A) 



r->oo X22(r, A) 



In the discrete case, the Schur function is defined as 

A n {z) 



lim 



The Schur function with opposite sign corresponds to the dual system with coeffi- 
cient of opposite sign. So, having (|181[) in mind (see the footnote below), we see 
why the opposite sign was picked up. 

If we view an operation 5 , z , ao ,...,a„_ 1 {w) introduced in (|160|) as a map of w € D 
to D with parameters z, do, . . . , a„_i € B, then 

• For z G T, it preserves the pseudohyperbolic distance. This is simply 
because both multiplication by z and the Mobius transform preserve this 
distance. 



We want to emphasize some abuse in notations in the definition of continuous Wall polyno- 
mials. If one wants to be consistent with discrete case, then the choice must be made according 
to I I180I I so that for the transfer matrix X in Krein system: 



X 



—A 



-A, 
B 



(181) 



In the meantime, we want to keep our notations to be consistent later on with terminology accepted 
in the scattering theory. 
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• For z € B, it acts as a contraction, i.e. 

p{S z<ao ,... tan _ 1 {w\), S z<a0l ... <an _ 1 (w 2 )) < p{wi,w 2 ) 

The contractive property follows solely from the contractive property of 
multiplication by z. 
Analogous properties for the map S\ a,t is given in the following 

Lemma 9.7. For any AgR, the map S\^A,r preserves the pseudohyperbolic metric 
and for A S C + it is contraction. 

Proof. For continuous A, the proof follows immediately from the properties of dis- 
crete map S and Corollarv l9.2l Approximating A € L 1 [0, r] by continuous functions, 
we get the statement of the Lemma in general case. □ 

Now, we can really regard the map S\^A,r as a combination of Mobius transforms 
and rotations. In the particular case A = 0, we have the following representation 

S ,aA z ) = ®o(AA,z) 
where we use notation $g (F, G, z) for the continuous continued fraction invented 
by Puig Adam [58] and later developed by Wall [76]. In this case, equation (|157p 
takes the form of the Riccati-Stiltjes equation 

®M = A(r) -A(rjf?(0) (182) 

We do not get deeper into this subject and refer the interested reader to the original 
papers. 

One should notice that there are many ways to approximate Krein system by the 
sequence of OPUC's. For example, one can take the following system of Verblunsky 
coefficients: 

hA(tx), 0, hA(t 2 ), 0, . . . , hA{t n ), 0,0,... (183) 
Then, the only difference will be a different scaling, e.g. 

X(r,X) = limM(0,2n,exp(zA/i/2)) 

h— >0 

There are at least two other ways to approximate the Krein system with the 
sequence of OPUC. They are discussed below and we will make use of them later 
on. 

Previously, we started with finite differences approximation to a system of ODE. 
That produced the approximation of the related analytic functions. Now we start 
with an accelerant and approximate it first. 

Theorem 9.3. Assume that we are given an accelerant H(x) € C[0,oo). Fix any 
R > and let h — R/n. Consider the Toeplitz matrices 



7j 



hH(-h) . . . hH(-jh) 



hH(jh) hH((j-l)h) ... 1 

For h small enough, 7 n > and it generates the Schur coefficients {aj^YjZ.Q such 
that 

lim sup \A(jh) - h- l a\ h) \ ^0 (184) 

h ^°8<]h<R 

where 6 is a small fixed number. 
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Proof. Indeed, we know that I + "Kr > 0, H(-x) = H(x), and H(x) G C[0,oo) 
Thus T„ > for n large enough and it generates the Schur parameters {a^ ^}»- = q 



_( h ) _ det 7j 



detTj_i 3 



(185) 



Consider the resolvent equation 



T R (x, 0) + / iJ(a; - w)r fl (M, 0)dw = H(x), 
o 

Its /i-step discretization leads to the system of linear algebraic equations with 
the matrix T„. If one takes any 6 < r < R, then the matrix T[ r /j-i] +1 is the 
discretization of the operator / + "K r but with the step of discretization (relative to 
the length of the interval [0,r]) slightly bigger than that for [0,i?]. Nevertheless, it 
tends to zero as h — > and this is why we need to keep r > S > 0. It allows us to 
use the following argument. We have 



T r (x, 0) + j H(x - u)T r (u, 0)du = H(x), 
o 

The discretization with the step h gives 

1 hH(-h) ... hH(-jh) 

hH(jh) hH((j-l)h) ... 1 

and j — [rh^ 1 ] + 1. Application of the standard arguments that use Hadamard's 
Lemma on the determinants yields 





" 7 (r,A)( ) " 




' H(+0) ' 




. l {r ' h \]) _ 




_ H(jh) _ 



lim sup 

/i^O fc=0 ,...j 



S r (kh,0) 7 ( r , fe)(fc) 







(186) 



where S r (x,y) and 8 r are introduced in Lemma 117.21 From this Lemma, we also 
know that 

1^,0) = ^} 



The convergence in (|186j) is uniform in r as long as S < r < R. Since A(r) = r r (0, r), 
we have 

sup |Zp)- 7 (r ' fe) 0')H0 

j:S<]h<R 

as h — > 0. At the same time, Kramer's rule gives us the following 

1 hH(-h) ... hH(-(j-l)h) H(+0) 



7 



(r,h) 



(?) 



1 



det 7j 



det 7j 
■det 



■ det 



hH(jh) hH((j - l)h) . . . hH(h) H(jh) 
hH(-h) ... hH(-(j-l)h) H{+Q)-h- 1 



hH(jh) hH((j-l)h) 



hH(h) 
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where the last formula follows from (|185j) . We have det 7j — > S r uniformly in r as 
h — > and 5 r is defined in Lemma ri7.2l Therefore, (| 1 84[) follows. □ 

The next Theorem is technical but we will need it later in the proof of the Strong 
Szego Theorem. It will give an approximation of Krein system through its measure 
although we need to take da of a very special kind. Assume that da is purely 
absolutely continuous with density 



27rcr'(A) = exp (^J 



l(x)e iXx dx 



with l(x)- Hermitian, continuous on R with compact support within [-R, R]. This 
measure da will generate the Krein system with continuous H(x) and A(r). For 
H(x), we have an expansion (see (|27p) 

* l](x) [I * ... * 1} (x) 



H(x) = l(x) 



21 



x 



Theorem 9.4. For large n, consider h = Rn 
a. c. measure dfi n on T with the density given by the formula 



k\ 

j = jh,j = 



(187) 
n and the 



H' n {9) = exp 



hl(xj)z^ 



j = -n 



Let a,j be the associated Verblunsky parameters. 



The 



lim sup \A(jh) 

6<jh<R 



(188) 



where 5 is any small fixed number. 
Proof. The 0-th moment of the measure is equal to 

h k 



„(o) _ 



= i + «(0) + ... + — 



/.■! 



l(x h ) ■ . . . ■ l(x jk ) + 



(189) 



which can be written as 

l + hH(Q) + o(h) 

as it follows from (I187j) , approximation of the integral by the Riemann sum, and 
simple estimates on the tail of the series (|189[) . The same is true about the higher 
moments, i.e. 

= hH{-kh) + 5(h), \k\ < n 
Moreover, \h~ 1 o(h)\ — > as h — > uniformly in \k\ < n. Application of the same 
arguments that proved Theorem 19.31 completes the proof. □ 



Remarks and historical notes. 

The continuous analogs of Schur and Caratheodory-Toeplitz problems were consid- 
ered in [45j . The corresponding classes of analytic contractions were introduced in 
the same paper. Equation (|158[) is rather standard in the theory of inverse prob- 
lems. For Schrodinger operators, the equation analogous to (|174[) was obtained 
and studied in |67j , [28] . The explicit approximation of Krein system by sequence 
of scaled OPUC's is new to the best of our knowledge. For discretization of con- 
tinuous Toeplitz operators, see [53] Chapter 8. It is a very good exercise to take 
A(r) — const on [0, R] and explicitly compute polynomials that correspond to dis- 
cretization with step h = R/n. 
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10. Zeroes of P(r, A) 

In this section, we study zeroes of the function P(r, A). It is convenient for us to 
write P(r,X) = exp(iAr)<3(r, A) and consider Q(r, A) instead. From (|55[) . we know 
that 

r 

Q(r,A) = l- / T r (s,0)exp(-iXs)ds 



o 

Assume that A(r) is not identically zero. By the Hadamard Theorem, function 
Q(r, A) has the following factorization 

oo 

Q(r, A) = CX m exp (a A) JJ (1 - A/A„) exp(A/A„), |A X | < |A 2 | < 
Since P(r, A) has zeroes in C + only, m — 0. Also, 

C = Q(r ' 0) ' Q= Q(r,0) 

Clearly, |Ai(r)| — > +oo as r — * 0. For each r > 0, zeroes A n (r) accumulate at 
infinity in a very regular way. For instance, if A{x) is smooth on x G [0, r] and 
A{r) ^ 0, then A„ has the following trivial asymptotics at infinity (see Lemma 1 17. 61 
in Appendix) 

A„ = A° +o(l),n-> oo, 
A° = x n + iy„, ^ + yl = \A{r)\ 2 exp(2r?/ n ), 
x n = r^ 1 [tt/2 + nn - Arg(A(r))} ,neZ 

i.e. the zeroes are accumulating evenly near the graph of the logarithm. Moreover, 
as r — > oo, the graph of logarithm is getting closer to the real axis and the spacing 
between consecutive zeroes decreases □• 

In the meantime, an interesting question is distribution of zeroes for finite r and 
A inside the compacts in C + . The Fejer Theorem for polynomials (f n (z) orthogonal 
on T with respect to measure d\i says that all zeroes of each (f n (z) are inside the 
convex hull of supp(d/i). Let us prove similar statements for P(r, A). Assume that 
supp(dcr) has a gap, say (a, b). We want to show that X n stay away from (a, b). Let 
M r be introduced by the following formula 

M r = sup|P(r, A) | 



LemmaEUand \P{r, A)| = |P»(r, A)| for A G R yield 



M r < exp 



\A(s)\ds 



,o 



Theorem 10.1. (Continuous analog of Fejer Theorem). Let (a, b) n supp(da) 
Then, P(r, A) has no zeroes in f2 r given by 

7 2yda(X) 2 

il r : A = x + iy, / jt rj- — j < M 

J (X-xy + y 1 

— oo 



''More on the asymptotics of {A n } can be found in recent preprint 1341 . 
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Proof. Assume that Ao = x+iy is a zero of P(r, A). Consider /(A) = P(r, A)/ (A — Ao) G 
By (H36j>, 

OO OO 

l/(A)| 2 ^ 



2vr|P(r, A)| 2 



That can be rewritten as 



1 

Thus, a simple estimate follows 



l^,A)| 2 
(A — x) 2 + y- 



\f{X)\ 2 da{X) 



■da(A) 



1 < Mi 



2yda(X) 
(A - x) 2 + y 2 



□ 



Clearly, f2 r contains a domain in C + contiguous to (a, b). 



Remark 10.1. One can modify this proof in the following way. For simplicity, 
assume (a, b) = (—1,1). Introduce 



N r 



\P(r,\)f 
A 2 + 1 

Then P(r, A) has no zeroes in the following set 



da(X) 



Q' r : A = x + iy,y < 



2N r sup 



A 2 + l 



AGsupp( CT ) (A - x) 2 + y 2 



For a large class of coefficients A(r), N r is bounded in r. This is because 2N r = 
Tr Im Gi (r, r) where G is the resolvent kernel for corresponding Dirac operator 
which will be introduced later. For example, A £ L°°(IR + ) is sufficient for N r to be 
bounded in r > 0. 

The next Theorem yields yet another result on the distribution of X n . 

Theorem 10.2. If Z\ is a zero of P{r, A), then there is no any other zero of P(r, A) 
in Oi 

Qi : X e C+, |A - zi \ < Dist(A,supp(do-)) 



Proof. Assume z\ is a zero of P(r, A). By the variational principle, function /o(A) = 
K r (zi,X)/K r (z 1 ,zi) minimizes ||/(A)||2, CT in the set of all /(A) G S r ,\f(zi)\ = 1. 
Since P(r,zi) = 0, 

(z 1 ~z 1 )P4r,X) 
(A - zi)P*(r, zi) 

In the meantime, if P(r, A) has a zero zj 6 Oi, then P*(r, I2) — and the function 

/i(A) = Z -f^fo(X) 

A — Z2 

belongs to S r , /i(zi) = 1, but ||/i||2,<r < ||/o||2,<r, a contradiction. □ 
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Notice that this Theorem makes no assumptions on coefficient A(r). It also 
implies that the isosceles triangle with base (a, b) and angles 7r/6, 7r/6, 2ir/3 can 
contain only finite number of zeroes. 

The next result is the continuous analog of the Widom's theorem on the zeroes 
of OPUC. It says that the zeroes of P(r, A) can not accumulate in the compact 
of C + provided that the support of da is not the whole R. The proof is a rather 
simple modification of proof for the discrete case. 

For any compact K C C + , define Njc(r) as the number of zeroes of P(r, A) in 
K. Fix any R > 0. We have elementary estimates 

max \P(r,\)\ < C(A,R,K), min \P(r,iy)\ > 1/2 

if y is large enough. Therefore, by Jensen's formula ([60], Theorem 15.18), we know 
that Nk{t) is bounded for r G [0, R]. 

Theorem 10.3. (Continuous analog of Widom's theorem). Assume that the mea- 
sure da of the Krein system is such that suppler) ^ R. Then, we have 

supiVif(r) < oo (190) 

Proof. Fix K and Krein system with the measure da. Cover K by disjoint cubes 
Cj with side e. We choose e small enough to satisfy the following conditions. For 
any cube Cj, consider £ G Cj and a map 4>^(X) — (A — The reflected cube 

Cj = {z,z G Cj} will be mapped to a set Dj ^ and the support of the measure 
supp(dcr) to a set Fj£, a proper subset of some circle. For each j, consider 

Dj = U ieCt D M , Fj = U 6eC] F 3 , 6 (191) 

We now require e to be so small that for each j we have: Dj and Fj are disjoint, 
C \ Fj is connected. We can always satisfy these conditions because the function 
0£ (A) is jointly continuous and supp(dcr) has a gap in it. 

Fix this e. Assume (|190|) is wrong. Clearly, among all cubes Cj there will be al 
least one, call it Cj' , such that for any k we can find r so that P(r, A) has n zeroes 
in Cji and n > k. Denote these zeroes by Xj,j = 1, . . . ,n. Fix this cube. Let 
Dji and Fj' be the corresponding sets defined by (|191[) . Let m < n be some fixed 
number to be specified later. 

By the variational principle, the function 

f m = ^r(Ai.A) Ai — Ai P,(r,A) 
Jo1 ' KriX^Xi) A-Ai P*(r,A x ) 

minimizes ||/||2,<t hi the set / G S r , |/(Ai)| = 1. We can write 

(A — A2) ... (A — A m+ i) 



fo (A) -.9(A) 



(Ai — A2) . . . (Ai — A m+ i) 



Notice that g(Xi) — 1. We will find a polynomial Q(X) satisfying the following 
properties: degQ < m, Q(Xi) — 1, and 



|Q(A)j < 2- 1 



(A-A 2 )...(A-A m+ i) 



(Ai — A 2 ) . . . (Ai — A m+ i) 



(192) 
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for any A £ supp (da) . That would give us a contradiction since for the function 
/(A) = g(X)Q(X) we have: /(Ai) = 1, 



J |/(A)| 2 da<2- 



l/o(A)| 2 da 



and 



/(A) - /o(A)Q(A) 



(A - A2) ... (A - A m+ i) 



(Ai — A2) ■ • • (Ai — A m +i) 
by Paley- Wiener Theorem. Thus we have a contradiction with the variational 
principle. 

To find Q(A), we first take a map z — (A — Ai) -1 . It sends Ai to infinity, the 
support of da will be mapped to \ lt and the set Cf (reflection of Cy with 
respect to K) will go to a compact Dy^. 

We now use Widom's Lemma (see Appendix, Lemma I17.7|) for two compacts 
Dj> and Fy . They are disjoint and C \ F.y is connected so the Lemma is applicable 
and the number m can be chosen so that for any points Zj £ D.j 
can find a monic polynomial Q(z) of degree to such that: 



1. 



, m, we 



Q(z) 



1 

< - 
- 2 



{z- z\)...{z- z m ] 

for all z £ F. In particular, for the points Zj = (Aj+i — Ai) -1 
, m there is Q(z) so that we have (|193p for any z £ Fy t Xi 



(193) 



£ Dy M C D,j = 

1, . . . ,m there is Q(z) so that we have (|193j) for any z £ Fy t Xi Q F . Translating 
it back to the A variable, we have (fT92]) where Q(A) = (A - Ai) m Q((A - Ai) -1 ). 
Clearly, degQ < to, Q(Ai) = 1. □ 

Remarks and historical notes. The asymptotics of zeroes for the exponential 
functions of the special type (e.g., P(r, A)) is a classical question. The problem 
here, of course, is how this asymptotics depends on the regularity of the function 
in representation (function T r (0,t) in our case). We do not consider this problem 
here, interested reader can check [34] for related results. 

The results from this section are new. We addressed only some of the basic 
questions about the distribution of zeroes. Clearly, there are many questions left 
open. 
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11. The case A(r) € L 2 (R+) 

Now, let us study an important class of Krein systems: A(r) G L 2 (R + ). 

Theorem 11.1. If A(r) G L 2 [0,oo), then da G (Szego). Moreover, 

P,(r, A) -» n o (A), P, (r, A) -» fi^(A), 2l(r, A) -» 21(A), 03(r, A) -> 03(A) 

as r — > oo uniformly in ImA > e, £ > 0. Function 03(A) G iV(C + ), f(A) 
03(A) 2l _i (A), ST 1 (A) is an outer function from B(C + ), 



ln|2l(A)|dA = 7T / \A{r)\ 2 dr 



(194) 



Proof. From 
P*(r,A) = 1 



, we have 

r s 

exp(iXs)A(s)ds + J exp(i\s)A(s) J exp(—i\t)A(t)P*(t, \)dtds 
o o 

(195) 



7 

P* (r, A) = 1 - J exp(iXs)A(s)ds 

r 

+ J P*(t,\) 



exp(—iXt)A(t) / exp(iXs)A(s)ds 



(196) 



dt 



Use Cauchy-Schwarz and Young inequalities to get 



cxp(— iXt)A(t) J exp(i\s)A(s)ds 



< T ^f,AGC+ 
1m A 



L 1 [0,< 



From Gronwall-Belmann inequality, 

|P*(r,A)| < C(A)ex P ([ImA]- 1 ||A||2) ,A G C+, C(X) < 1 + \\A\\ 2 [Im A]" 1 / 2 

Recall the function (|146|) and Lemma [831 From (|196|) , we have P*(r, A) — > n a (A) 
as r -> oo, A G C+. Similarly, P*(r, A) -> S /3 (A). Thus, 2l(r, A) -> 21(A), 03(r, A) -> 
03(A) as r -> oo. Moreover, from Lemma |6"TTI |2t(A)| 2 > 1 + |03(A)| 2 , ImA > 0. 
Consequently, ST 1 (A) G P(C+), 21(A) G N(C+), 03(A) G N(C+). We also have 



21(A) = (II Q (A)+IL3(A))/2 



n Q (A) 



n^A) 
n Q (A) 



Clearly, II Q (A) and II/^A) are outer from 7V(C + ). Due to Theorem 16.2 

ft/3 (A) 



P(A) 



n Q (A) 



Since P(A) has positive real part, the function 1 + IL3 (A)II a 1 (A) is outer from 
N(C+). Consequently, 21(A) is outer from N(C + ) and 2l _1 (A) is outer from P(C+). 
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Thus, 



2t _1 (A) = exp 



1_ f (1 + sA) ln|g(a) 

TTi J (S- A)(l + S 2 ) 



■ds 



,7€[0,27r), (197) 



For 2l(r, A), we have 

2t(r,A) = 1+ / »(f,A) 



exp(-iXt)A(t) / exp(iXs)A(s)ds 



ds 



(198) 



Iterating this identity, estimating the lower order terms, and taking r — > oo, one 
has 



= i + 

Clearly, 



exp(yt)A(t) / exp(— ys)A(s)<is 



+ <3(y 2 ),y^+^ 



exp(yt)A(t) / exp(— ys)A(s)ds 



y-icv 



du> 



where A(oj) is the Fourier transform of the function A(t)-Xjn+ (t) and |^4(o;)| 2 S L 1 (M). 
Thus, 

OO 



a(i») = l + y- 1 J \A(s)\ 2 ds + o(y" 1 ), y -» + 
o 

From (|197[) . we have 

OO 

y /■ In | 21(A) | 



|2l(iy)| = exp 



A 2 + 



dX 



(199) 



(200) 



Since | 21(A) | > 1 for a.e. A e R, relations (fTM)) and |(2DD} imply In | 21(A) | € L^K) 
and ifTM]) . □ 



Corollary 11.1. For a.e. A G K, we ftave 

|2i(A)| 2 = l + |<B(A)| 2 
Proof. The equation is equivalent to 

|n a (A) + n / ,(A)| 2 _ |n a (A)-n>(A)| 2 

4 4 



(201) 



or 



Re F(X) 



1 



27rcr'(A) 



|n Q (A)| 2 

and the last identity is elementary and follows, e.g., from the integral representa- 
tions for both functions F(A) and II a (A). □ 



Corollary 11.2. If A(r) € L 2 (R+), then 

OO OO 

2tt J \A(r)\ 2 dr = - J ln[l - | f(A)| 2 ]o!A 



(202) 
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Proof. From IpOl]) . we have 1 - |f(A)| 2 = |2t(A)|- 2 for a.e. A £ I. Now, (12021) 
follows from (fl94l . □ 

Notice that ([2021 implies 

f(A) € L 2 (R) (203) 
So, f(A) € H 2 {C+) n S(C+) and (compare with (11651) ): 



f(A) = / C(x)exp(iAx)dx,A € C+ (204) 
o 

Corollary 11.3. If A(r) E L 2 (R+), then ln[27rc/(A)] € i 1 (R) + L 2 (M). 
Proof. We have for a.e. A S K 

[2™' (A)]" 1 = |n Q (A)| 2 = | 21(A) + <B(A)| 2 = | 2l(A)| 2 • |1 + f(A)| 2 
and therefore 

- 1ii[2trj' ( A)] = 2 In | 21(A) | + 2 In 1 1 + f( A) | 

By (|194| . the first term is from L 1 (R). As about the second term, 

ln + |l + f(A)|<ln[l + |f(A)|]<|f(A)|eL 2 (M) 

by (|203|) . For the negative part of the logarithm, we use an elementary estimate 
1 1 + f | > 1 — | f | which yields (check the definition of ln~ ) 

| In" |1 + f(A)|| < ln(l - | f(A)l)- 1 = - ln(l - | f(A)| 2 ) + ln(l + | f(A)|) 

for |1 + f(A)| < 1. The second term is again in L 2 (R). The first one is in L X (K) due 
to Corollary E21 □ 

Corollary 11.4. For the constants a and 7 from the multiplicative representations 
o/IT Q (A) (formula (fl48l) ) and 21(A) (formula (fl97l) ) we have 

1 ? sln[27rcr'(s)l , 1 7 sln|2l(s)| , 
a = — / - as, 7 — — / - — as 

27T J 1 + S 2 ' ' 7T J l + s 2 

— OO — OO 

Proof. From (fTW|> and pM)) . we have 

OO OO 

_ 1 Bm r wwi, _ 1 I .NiWl t 

7T J (s 2 + 1J 2 )(1 + S 2 ) 7T J 1 + S 2 

— OO —OO 

because ln|2l(s)| e L 1 (M). 

Take r — ► 00 in (|196[) . One has 

n Q (A) = l + 0([ImA] 1 / 2 ]) 

as ImA — » +00. Then, from (| 148[> . we have 

1 [ s(l-y 2 )+iy(l + s 2 ) : ro 

a= urn — / — — ttttz — — m\Ziro (s)\ds 

y^2nj (s 2 + y 2 )(l + s 2 ) 1 WJ 

Using Corollary 1 11.3( we get the needed formula for a. □ 
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Clearly, the Corollary implies the following integral representations 

oo 

1 f ln|2t(s) 



21(A) = cxp 



n a (A) = cxp 



-ds 

tti J s — A 

— oo 

oo 

1 f ln[27rcr'(s)] 



2ni 



A 



ds 



(205) 



(206) 



As we know from the discussion of the general Szego case, the function (A + 
?') _1 n^ 1 (A) 6 H 2 (C + ). For our situation, much more is true 

Theorem 11.2. If A(r) G L 2 {R+), then 

OO 

n- 1 (A) = l+ (\(x)e lxX dx = l + j(X)^(X)eH 2 (C+) (207) 



where 



a{E) + / \ 7 (x)\ 2 dx = / \A(r)\ 2 dr 



(208) 



where E- support of da s , the singular component of da. 

Proof. From Lemma 18.91 as r — > oo: 

P,(r,A) n Q (A)- Xg c(A) 
A + i X + i 

in L 2 (R, da), where E c - the complement of E. On the other hand, for any r > 
0, A e K, we have (by ([57])): 



That implies 



r 

P* (r, A) = 1 - J A(t)P(t, X)dt 



U a (X)-XE^X) 1 i(A) 



(209) 



X+i X+i X+i 

where the generalized Fourier transform 

oo 

A(X) = / A(r)P(r,X)dr e L 2 (M,da) 



by Theorem O Thus 
and 



A(X) = 1 — II a (A) • Xb c (A) 

OO oc 

| |,4(r)| 2 dr = (7{i5} + -L | 



(210) 



n Q (A) 



Now, we have 



n Q (A) 



dX < 



1 - 



x + i V n Q (A) 



dX 



e h 2 (c a 
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Therefore, by elementary Lemma [17.91 in Appendix, 1 — H a 1 (X) <E H 2 (C + ). The 
Paley- Wiener Theorem now implies (|207l) and (1208]) . □ 

An interesting corollary from this result is that the total variation over the whole 
line of the singular part of the measure da is finite. 

After proving the representation for n Q (A), the following result is quite natural. 
By Levy- Wiener Theorem, 

oo 



= 1 + / lr (x)e iAx dx 



P.(r, A) 

where j r (x) € L 1 (M+) n L 2 (R+). 

Lemma 11.1. Assume A(r) e L 2 (R+) and E = 0. Then, j r (x) -> j(x) tnL 2 (R+). 

Proof. From (|209|) and Corollary 16. 11 we get 

r oo oo 

J\A(x)\ 2 dx= J \l-P«( r ,\)\ 2 da(\) = ± J |1 - P,(r, A)| 2 , p (211) 



|PM))p 



1 

2t7 



P*(r, A) 



d\ = \j r (x)\ 2 dx 







Then, since P*(r, A) — > n a (A) for A S C + , we also have that 7 r (a;) — ► 7(2;) weakly in 
i 2 (R+). But since E = 0, we also get ||7r(x)|| 2 -> ||7(x)|| 2 from (j2Hgj) . Therefore, 
7 r (x) -> 7(x) in L 2 (R+). □ 

If the singular part of the measure da is not trivial, then we have only the bound: 
117,112 < ||A|| 2 . 

Now, let us characterize the class of Schur coefficients corresponding to A(r) E 
L 2 (R+). We start with 

Theorem 11.3. For any A(r) £ L 2 [0,R) and any R > 0, we have 

oo R 

\n\%(R,\)\d\ = n J \A{r)\ 2 dr (212) 

-oo 

and 

R oo 

2tt J \A{r)\ 2 dr = - J ln[l - | Q3(P, A) ST^P, A)| 2 ]rfA (213) 



Proof. The proof repeats those of Theorem 1 11. II and Corollary [Tl~2l □ 

Let us introduce a certain subclass of P(C + ). Consider a function f(A) € P(C + ) 
such that for its boundary value: 

ln(l - |f(A)| 2 )dA > -oo (214) 
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Then, a simple estimate ln(l - | f(A)|| 2 ) < -| f(A)| 2 yields that f(A) G L 2 (R). Thus 
f G H 2 (C+) and 

OO 

f(A) = J C(x)e lXx dx 
o 

with C(x) G L 2 (R+). 

Definition 11.1. We say that f(A) G S (C+) if f (A) G B(C+) and (pTi]) holds. 

The set 5o(C + ) can be regarded as the metric space [71] with the distance given 
by the formula 



ln[l-p 2 (f(A), (A))] d\ 



(215) 



where pseudohyperbolic distance p(-,-) is defined by (1159p . It turns out that the 
resulting metric space is complete (see Lemma 1.5, Theorem 1.6, Corollary 1.9 in 
|71j . The geometry of this space is studied in the same paper). 
Let us write p s (f) = p s (0, f) for short-hand. 

Lemma 11.2. If in Krein system f r (A) G 5o(C + ) for some r, then f(A) G Sq(C + ) 
and 

p 2 (f(A)) =^(®(r,A)2l- 1 (r,A))+p 2 (f r (A)) (216) 
In particular, this is true for any A(r) G L 2 (R + ). 

Proof. Let us prove this lemma using certain "orthogonality" argument. For A £ K, 
we have 

21* (r, A) = e lAr 2t(r, A), 05* (r, A) = e iA, '25(r, A) 
Therefore, (|117p yields 



I f(A)| s 



«8(r,A)a (r,A) + exp(iAr) f r (A) 



1 + 05(r, A) 2t _1 (r, A) exp(zAr) f r (A) 



, A g 



The formula 

In [ 1 



z + w 



1 + zw 
and Theorem 111.31 give 



ln(l - \z\ 2 ) + ln(l - |w| 2 ) - 21n|l + zw\ 



p 2 (f(A)) = p 2 (<8(r, A) STV, A)) + p?(f r (A)) 

OO 

+2 / ln|l + ( 8(r,A)a" 1 (r,A)exp(iAr)f I ,(A)|dA 



p 2 (<B(r,A)2l- 1 (r,A)) + p 2 (f I ,(A)) + 2 / In |1 + 58, (r, A) ST^r, A) f r (A)|dA 



The last integral is zero because 

95, (r, A) VT\r, A), ST^r, A), f p (A) G B(C+), 
05*(r,A) = o(l) as A G C+, |A| -> oo, 05* (r, iy) = d{y~ 1/2 ) as y -> +oo 
f r (A) G H 2 (C + ), so f r (iy) = ^y- 1 / 2 ) 
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functions fB*(r, A), f r (A) G L 2 (M) and the mean-value formula (Lemma Il7.8| Ap- 
pendix) is applicable. □ 

The result above is sometimes called "the layer stripping". That is because 
p 2 (f) is equal to the sum of the terms that correspond to different intervals of the 
coordinate r. The formula (|202[1 is called the non-linear Plancherel Theorem. The 
both results are well-known in the theory of orthogonal polynomials. 

It is an important observation that any function from 5*0 (C + ) gives rise to a 
certain Krein system. Indeed, we can show So(C + ) C 5(C + ), where the class 
5(C + ) was introduced in the Dcfinition l9.H and the Theorem 19 . 1 1 applies . To prove 
this inclusion, notice that C(x) generates operator C/j (see (| 164[) ) for any R > 0. 
By a standard argument, HC^H^o^] < 1 for all < R < oo. Then, for any R > 0, 

R 

f(X) = J C(x)e iXx dx + e aR $ R (X),\ G C+ 
o 

with 

oo 

$ R (A) = J C(x + R)e tXx dx 
o 

Since <&r(\) € H 2 (C + ) by Paley- Wiener Theorem and both 

R 

/(A), / C{x)e tXx dx G H oc {C+) 



we get $ H (A) G L°°(R). Therefore, $ H (A) € H°°(C + ) and /(A) G S(C+). 
In [71], the analog of the following Theorem was established. 

Theorem 11.4. For Krein system, A(r) G L 2 (R + ) iff the corresponding Schur 
junction /(A) G S^C" 1 "). 

Proof. If A(r) G L 2 (R+), then f(A) G 5 (C+) follows from the Corollary [TT21 
Assume now that f(A) G <So(C + ). The observation made right before the Theorem 
says there is the corresponding Krein system with A(r) G £ 2 OC [0, oo). Let us show 
that actually A(r) G L 2 (R+). Indeed, fix any R > 0. Then, 

f(A)a(fl,A)-fB(fl,A) 
Tfll j 21* (iJ, A) - f (A) 05* (R, A) 1 J 



and 



Therefore, 



f fl (A)| 



f(A) - *B(i?, A) 2t _1 (i?,A) 



1 — f (A) gj(i2, A) 2t _A A) 



,AG 



Ps (f R (X)) = p s (f(A), A) St" 1 ^, A)) < oo 
and ffj(A) G 5o(C + ). Now, the Lemma 111.21 is applicable together with Theo- 
rem QX1 

R 

2t: J \A(r)\ 2 dr = p 2 s (<B(R,\)%- 1 (R,\))<p 2 s (ft<C (218) 
o 

uniformly in R which means A(r) G L 2 (M. + ). □ 
This result has an interesting corollary. 
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Corollary 11.5. (i) If f(A) is the Schur function corresponding to Af(r) £ 
L 2 (R + ) and g(X) £ B(C + ) then gf generates the Krein system with A g f £ 



and \\A gj \\ 2 < ||A f || 2 . 



(ii) The set of measures da that correspond to Adair) £ L 



is convex. 



Proof. The first statement is obvious due to Theorem 1 11. 41 The second one follows 
from the following calculations. 



1-1/1' = 



1 + F'- |1 + F| 2 

If Fj and f ^ correspond to daj , j =0,1 and F t and f t - to icier i + (1 — t)da ,t £ [0,1], 
then 

<|1 + Fi| 2 



4(<RcFi + (l-t) RcF ) 
ll + tPi + a-t^ol 2 



= (i-lAr)- 



l + tF 1 + (l-t)F \ 2 



+ (l^l/o| 2 )- 



(l-i)|l + ^o| 2 



l + tF! + (l-i)F | 2 -j=o,i 



> minfl-l/,! 2 ) 



due to convexity of \z\ , z £ C. Therefore 



ln(l - |/ t (A)| 2 )dA > / ln(l-|/ (A)| 2 )dA+ / ln(l - |/i(A)| 2 )dA 



and ^ (A) S S , (C H 



□ 



Remark. Clearly, the last estimate is not optimal. As about the first statement, 
notice that multiplication f A by any inner function does not change L 2 norm of the 
coefficient A(r). 

The Theorem 111.21 leads to the following natural question: what can be the 
singular component of da if A(r) £ L 2 (R + )? The answer is given by the following 
result which can be regarded as another criteria for A(r) £ L 2 (M + ). In particular, 
it says that the singular component can be any singular measure finite over M. + . 

Theorem 11.5. Let da be a nonnegative measure on R with decomposition da = 
da s + a'(X)dX, where a s (R) < oo and ln[27rcr'(A)] £ L 1 (M) + L 2 (R). Assume also 
that 



exp 



1 

2tt7 



ln[27r</(*)] 
t-X 



dt 



Then da generates the Krein system with A(r) £ L 2 ( 
Proof. We have 



27rcr'(A) = 1 + 7 + 7+ |7| 2 
Take H(x) Hermitian such that 



- 1 = J 7(1) exp(iXx)dx = 7(A) £ H 2 (C+) (219) 

and (|208[) holds true. 

(220) 



H(x) = (2tt)- 



+ / exp(iAa;)(icr s (A) 



(221) 



7(2;) + / j(x + u)j(u)du 


for x > and define (3 by (j28|) . Then, it is not difficult to check that the formula 
H27]) holds. Moreover, H(x) £ L 2 oc {R) is an accelerant and generates some A(r) £ 



functions. 



That follows from ([50]) . (|220p . and the uniqueness theorem for analytic 
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Let us show that A(r) € L 2 (R + ). Indeed, take da n such that it has only finite 
number of jumps and 

exp(i\x)da n (X) — > / exp(iAx)d(T s (A) 



uniformly in a; S [0, i?] with any fixed i? > 0. 
Also, take purely a.c. d[i e so that 

2tt/4(A) = |l+7e(A)| 2 ,7 £ (A) =7(A + *e),e> 0, 

Take dcr„ !C = da n + d/j, e . Then, the corresponding accelerants H n ,e(x) — ► i?(a;) and 
^n.efV) — > in L 2 [0, R] for any fixed R > as long as n — > oo, e — > 0. 

Let us show that Theorem 111.41 yields A„ i£ € L 2 (R + ). Indeed, fix n and e. Then, 
we can write dcr nje as a convex combination 



da n .f = t 



da n dX 
~ + 2tt 



+ (1-*) 



// e (A) t 



dA 



1-t (l-t)27T_ 

and t G (0,1). If we can show that both measures in this convex combination 
generate square summable A(r), then the second claim of the Corollary QTT5] finishes 
the argument. We will apply Theorem lll.4l to 

da n dX 
~T + 2t7 

and 

The first measure gives rise to 

da n (t) 



F(X) = 1 - 2ii _1 / 



t — A 



This representation follows from the formula (|122[) and F(iy) — ► 1 as y — > +oo. 
Since 

i-ifi"-^ 

and (i<T n has only finite number of jumps 

l-|f| 2 = l + 0(|A|- 2 ) 

as |A| — * oo. The local singularities of ln(l — | f | 2 ) are integrable and thus we have 
ln(l - |f | 2 ) e L*(R) and the Theorem [TOl can be applied. 

Let us show that Theorem II 1 .41 can also be applied to the measure (|222l) as long 
as t is chosen properly. Indeed, we can always take t so small that the density u(X) 
of this measure is strictly positive. We have 



FA =l + i 7- -T- r ds 

Jr n(l-t){\- s) 

Since % is infinitely smooth, F is continuous up to the boundary and ln(l — | f | 2 ) 
is locally integrable by (|223[) . Therefore, we are left with showing that | f(A)| — * 
if A e R, A -> oo and that | f(A)| € L 2 (R). Since %(s) € L 2 (R) and it decays at 
infinity, the simple properties of Hilbert transform imply F(X) — 1 € i 2 (R) and 
F(A) — l^OasA^oo. Since 

£ 1 ~ F 
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we have ln(l - | f | 2 ) g L 1 (R) and Theorem [l 1 . 41 applies . 

The Krein system with coefficient A n ^ has II £ (A)-function with inverse 

RJ 1 (X)=1 + %(X) 



Then, notice that (|208|) gives a bound ||^4 n ,e||2 < C uniformly in n and e and 
A{r) S L 2 (R+) because A n<e {r) -> A(r) in L 2 [0,P] with any i? > 0. We get (208)) 
and the function 7(2;) coincides with the one introduced in Theorem 1 11. 21 □ 



Later on, we will need the following bound which sharpens the Lemma 18.9 
Lemma 11.3. If A(r) E L 2 (R+), then 



1 



A 2 + 1 



n Q (A) 

— CO 

Proof. First, we notice that 



dX 



1 



A 2 + l 



|P*(r, A)| 2 da s (X) < C / \A(s)\ 2 ds 



Therefore, 



00 

n a (A)-P,(r, A) = - J A{s)P{s,X)ds,XeC+ 

r 

|n a (A)-P(r,A)| 2 < (J\A(s)\ 2 d S \ ■ (J\P(s,X)\ 2 ds 



< 



1 



21mA 



\A(s)\ 2 d S \ ■ [|II a (A)| 2 -|P*(r,A)| s 



by Cauchy-Schwarz and analog of ([59| . For the last expression, we use 

|n a | 2 H^*l 2 <|na-P*h(|n Q | + |p*|) 

and then 

00 

|n Q (A) — P*(r, A) I < C(A) /" |A(s)| 2 cis,IniA > 

r 

The last estimate also implies 

\P( S ,X)\ 2 ds = (21mA)- 1 (|II Q (A)| 2 - |P(r,A)| 2 ) 



(224) 



(225) 



< C(A)|n Q (A) — P*(r, A)| < C(A) / |A(s)| 2 ds 
For P(r, A), we have 



(226) 



P 2 (r,A) = -2 P'(s,X)P(s,X)ds = -2 / iXP 2 (s, A) - A(s)P(s, A)P*(s, A) ds, A e C+ 
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Therefore, for A € 



\P(r,\)\ 2 <C(A) / \P( Sl X)\ 2 d s + C(X) 



\A{s)\ 2 ds- / \P{s,X)\ 2 ds 



1/2 



< C(A) / \A{s)\ 2 ds 

r 

where the last inequality follows from (|226|) . 

Now, let us improve estimates from Lemma [8781 From (I149p . we have 



(227) 
(228) 



\P*(r,iW 



l P -^ A )'-MA)= |n( ^ )|2 - |P(M)|2 + lP(M)l 2 



A 2 + 1 



A 2 + l 



dcr(A) 



2 Re 



P(r, A)P(r, i) 
A + z 



AT r (i, A)dcr(A) 



For real A, \P(r,X)\ = |P*(r,A)|. Therefore, 

IP (V A"ll 2 1 r 

' Ar(A) = -+2[|P*(r,i)| 2 -|P(r,i) 

— oo 

Using the representation 



Re 



iP(r, i) 



P(r,A) 
A + i 



(229) 



K r (i,X)da(X) 
(230) 



# r (i,A) = Xoo(z,A)- / P(s,X)P(s,i)ds 



and the property of reproducing kernel, we get 



P(r,A) 



-K r (i,X)da(X) 



P(r,i) f P(r,A) 
~~ 2i J X + i 



P{s 1 X)P{s,i)dsda(X) 



The last integral can be bounded by Cauchy-Schwarz and Theorem I3.2I as follows 



P(r,A) 
A + i 



P(s,\)P(s,i)dsda(X) 



<C(A) 



\P{s,i)\ 2 ds 



- 1/2 


oo 


<C(X) 









1/2 



where we used (|226[) for last inequality. Estimate (|228p and (|230[) yield 



\A(s)\ 2 ds 



(231) 



Now, repeating the proof of Lemma [8?9l with estimates (|225[) and (|23ip . we obtain 
(12241). □ 



The estimates in the last Lemma are not sharp but good enough for us. 

Remarks and historical notes. The results in this section are partially new. 
For the Helmholtz equation, analog of Lemma 111.21 was obtained in 71J where 
the nonlinear Fourier transform was introduced. In our case, this transform is 
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given by the map J : A(r) € L 2 {WL+) — ► f(A) e 5 (C+). In [71\, the space 
So(C + ) is studied in detail as well as properties of the nonlinear Fourier transform. 
For example, its homeomorphic property is proved by means of weak convergence 
argument. See also [73) . For the Schrodinger operators and Jacobi matrices, the 
analysis is more involved [12]. The recent paper [42], contains analysis of L 2 (R) 
potentials for Schrodinger operators. In the OPUC theory, many of these results 
were well-known for quite a long time. The paper 16J studies the relation between 
the decay of the tail ||A||^ 2 r % and the Hausdorff dimension of the support of da s . 
In conclusion, we want to say that the case of square summablc coefficient is studied 
pretty well by now. Perhaps, the only problem left open is the following nonlinear 
(non-commutative) analog of the Carlesson Theorem in the Fourier analysis. Prove 
(or disprove) that solution of the ODE: P^(r, A) = A(r) exp(i\r)P*(r, A), P*(0, A) 
1 has a limit at infinity for a.e. A 6 R. We assume here, of course, that A(r) € 
L 2 (R + ). This is a deep and difficult problem whose analog for OPUC case is also 
open for quite a long time. We mention the paper [10] for some recent closely 
related results in this direction. 
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12. Continuous analog of the Baxter theorem. The case 

A(r) G L 1 (M+) 

In this section, we assume that A(r) and H(x) are both from regularity class 
L 2 oc (K + ). Our goal is to prove the following analog of Baxter's theorem in the 
OPUC theory. The proof is an adaptation of the one for the discrete case ([66], 
Chapter 5). 

Theorem 12.1. For any Krein system, A(r) € L 1 (R + ) n L 2 oc (lR + ) if and only if 
the accelerant H(r) G L (M) n Lf oc (M. + ) and the Hopf-Wiener operator 

oo 

(I + Xoo)f = f(x) + J H(x - y)f(y)dy 
o 

is strictly positive on L 2 (R + ). The last condition is equivalent to 

l + p(A) > 0,A > 



where p(A) 6 W(M) is the Fourier transform of H(x). Moreover, the measure da 
is purely absolutely continuous, has continuous derivative and 

exp(-2||,4||i) < 27rcr'(A) = l + p(A) < exp(2||A||i) (232) 

Proof. Assume that we are given A(r) S L 1 (M + ). We are then in the Szego case. 
Indeed, an elementary application of Gronwall-Bellman inequality to (|80p yields 
the uniform convergence of P* (r, A) and P* (r, A) to II Q (A) and n Q (A) in A € C+ . 
Both iIo,(A) and n a (A) are continuous in C + . Moreover, 

|P*(r,A)| < cxpdlAHi), |P*(r,A)| < exp(||A||i), A G C+,r > 
\P*{r t \)\ > exp^PHx), |P(r, A)| > exp( — 1| A|| i), A G C+,r > 

and then 

exp(-p||i) < |n (A)| < exp(||A|| 1 ),exp(-||A||i) < |fi a (A)| < cxpO^id) 

The Theorem 16 . 2 1 savs that measures (27r) _1 |P t (r, A)|~ 2 <iA converge to dtr(A) in the 
weak-(*) sense. Therefore, da(X) is purely a.c. and its continuous density allows 
an estimate 

exp(-2||A||i) < 27rcr'(A) = |n Q (A)T 2 < exp(2||A||i) (233) 

Now, let us show that H{x) G L 1 (R + ). Indeed, let 

2/i (r, A) = P, (r, A) - 1, y 2 {r, A) = P(r, A) - exp(iAr) 

For each r > 0, y 1{2 ) G W+(C+) and \\yi\\ w+ = hj2\\\v + by ([36]) and ([37]). From 
([57]) . we get 



yi (r)=a(r)- / A(s)y 2 (s)ds (234) 

where 

a ( r ) = - I A( s ) exp(is\)ds G W + 
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and (|234|) is considered as an integral equation for functions with values in W+. 
Taking the norm of the both sides in ()234|) . we get 

r r 

\\yi(r)\\w + < J \A(s)\ds + J \A(a)\-\\ yi {a)\\ w+ da 
o o 

The Gronwall-Bellman inequality yields convergence of y\ (r) to some y\ in the W+ 
norm (as r — ► oo) and 

\\yi\\w + < PHiexpdlAlli) 

Therefore, n a (A) = 1 + y\. In the same way, we have ila(A) = 1 + y\ and F(X) — 
U&/IL a = (1 + yi)/(l + yi ) = l + h, where 

h = h - V - e W+ (235) 

i + yi 

since the spectrum of 1 + y% does not contain zero. Now, ([50)) implies 

oo 

h = 2 J H{x) exp(i\x)dx (236) 
o 

and so H(x) £ L 1 (R) since heW + . We also have 2ira' (X) = 1 + p(A). 

Now, assume that we are given H(x) € L 1 (R) and 7 + 3-Coo > 0. Then, the 
equivalence of I + 3-Coo > and an estimate 1 + p(A) > follows from the simple 
identity 

oo 

((7 + JC 00 )/ ! /) = 2tt / (l + P (-A))|/(A)| 2 dA 



Last identity shows that H generates a Krein system with A(r) £ Lf oc (M. + ). To- 
gether with formula (|50|) and a simple approximation argument (like in the proof 
of Lemma 1 1 7 . 1 01 in Appendix), it also imply da = (27r) _1 (l + p(X))dX. 
We need a simple 

Lemma 12.1. < 1, then A(r) e L 1 (IR + ) and \\A\\i < ||-ff||i/(l - ||77||i). 

Proof. From (|34|l . we have 



T r (t, 0) + / F(t - u)r r (u, 0)du = H(t) (237) 
o 

Iterating this identity, we have for 

r r r 

T r (t, 0) = H(t)-f H{t-ux)H{u x )dui+ J H(t-m) J H(u 1 -u 2 )H(u 2 )du 2 du 1 -. . . 



(238) 

where the series converges in L x [0,r]. Then, |r r (t, 0)| < g(t), where 
g = h + h*h + h*h*h+ h(t) = \H(t)\ 
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Notice that g(t) does not depend on r and \\g\\i < ||i?||i/(l— ||i?||i). So, by taking 
t = r in ([237]) . we get 



A(r) + / H(r- u)T r (u,0)du = H{r) (239) 



\A(r)\ < J \H(r - u)\ ■ \g{u)\du + \H(r)\ (240) 
o 

The Young inequality finishes the proof. □ 

To finish the proof, we will apply this Lemma to the interval [R, oo), where R 
is so large that the accelerant for the Krein system considered on [R,oo) has small 
L 1 norm. To prove the existence of such R, we need to use Baxter's Lemma (see 
Appendix, Corollary 117. 2[) to the operator / + "K^ acting in the space L 1 (R + ). 
Since H(x) £ L : (R) and 1 + p(A) > 0, this Lemma is applicable and gives 

||(/ + Jt r )" 1 || i i[o ir ],i,i[o, r ] < C, (r > r ); \\T r (0,x)x[o, r ](x) -r(x)||i -> 0, (r -> oo) 

(241) 

where 

r(ar) = (I + 'K ao )- 1 H{x) e L 1 (R+) 

is solution to the Wiener-Hopf equation. Then, we are immediately in the Szego 
case since 

r oo 

P*(r,A) = 1 - / T r (0,s)exp{i\s)ds — ► II Q (A) = 1 — / T(s) exp(i\s)ds 



and convergence is uniform in C+. 
For F(X), we use (fTT9| and 



X-t i 2 + l 



1 + ^-dt 



2rr 



F(A)/2 = -i(3 + i 
Since F(ioo) = 1, we have 

oo 

F(A) = l + i | ^PW* 

— OO 

and the integral is understood in v. p. sense. That immediately implies 

oo 

F(\) = 1 + 2 J H(x) exp(i\x)dx 


From F = F^ 1 we have H(x) G L 1 (M) for the dual accelerant and 

||f r (0,a:)x[o,r](a!) - f (z)||i - 0, (r - oo); f (as) = (I + Woo)" 1 ^) G i 1 (R+) 
So, we also have 

|| B(r, A) - <8(A)|| W+ - 0, || 2l(r, A) - 2l(A)|| w+ - (242) 
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where 

oo ^ OO ^ 

21(A) = 1 - J F(x) + r(x) exv{i\x)dx, 05(A) = J ~ F(x) exp(zAx)^ 



From (fI73|) . 

exp(Ui?) f fl (A) = [f (A) 2l(i?, A) - B(J2, A)] • [2l(fl, A) + f fl (A) 23* (fl, A)] (243) 

Consider the first factor in the right-hand side. It can be written as 

a _1 (A) [23(A) 2l(i?, A) - 23(i?, A) 21(A)] 

Since 21(A) — 1 E W+ and 21(A) has no zeroes in C + , relations (|242[) imply this 
factor goes to zero in W + norm as R — > oo. On the other hand, this factor can be 
written as 

23(i?,A)" 



2l(i?, A) 



f(A)- 



(244) 



2l(i?, A) 

and (|173[) says that the Fourier coefficient of the second factor in (|244[) is equal to 
on [0,R]. Thus, 

f(A) 2t(-R, A) - 23(i?, A) = exp(i\R)Q R (\) 
where ||Q.r(A)||w + — > 0. Now, 

Q R (\)%(R, A) 



and ||/i?.(A)||v7 + -> since || 2l(i?, A)|| w+ < C, || 23*(.R, A)||w+ < C 1 uniformly in R. 
Since F fi = (1 - f R )(l + f a ) _i , we also have jjl - F R \\ W+ -» or ||lffl(a:)||i -» 
where H R (x)— accelerant corresponding to the interval [R, oo). The Lemma 112. II 
now yields A(r) S L x [i?, oo) as long as \\H R \\ < 1. □ 

Let us obtain the formula for a in the representation for IIq^A), just like we did 
for the square summable A(r). Since II a (A) = l + yi(A) and yi(A) G W + , we have a 
trivial asymptotics: LT Q (A) = 1 + o(l) as A — > oo, A G C+. From the multiplicative 
representation for II Q , we get 

i_ /■ s(i_y^±Ml±£ 2 ) 

y "*<i 27T 7 (s 2 + y 2 )(l + s 2 ) 
and simple estimates yield 



01 = J im _ oZ I — ~ a 2 T „2Vi a, C 2N — M(s)ds, /i(s) = ln(l + p(s)) G 



3C 



1 /" Sfl(s) 

— v.p. / ^ds 



2tt 7 1 

— oo 

Thus II Q (A) allows the same representation (|206[) . Notice also that in contrast with 
square summable case, function 21(A) does not allow asymptotical formula (]199|) . 
It should also be mentioned that the Theorem does not provide a quantitative 
estimate on ||A||i in terms of and, say, |j(l + p)~ 1 \\ 00 - On the other hand, it 

is easy to bound ||(1 + p) _1 ||oo in terms of ||A||i. 
Later on, we will need the following result 

Lemma 12.2. Assume that conditions of Baxter's Theorem hold. Then A(r) G 
C (R+) iffH{x) G C (R+) iffC{x) G C (R+). 
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Proof. Since H(x),C(x) £ L 1 (R+), we have H(x) £ C* (M+) iff C(x) € C (K 
because of (|T68| . Assume € C* (K + ). Then, ([240]) implies A(r) £ C (M H 

Now, let A(r) G L 1 (R+) n C (M + ). Consider (81]). Since A(r,f) = T r (t,r) 
r r (0,r-t), 



T r (0,t) = A(t) - I A(s)-T s (0,t)ds 
t 

Take R large and iterate this identity for Ar = {R < t < r}. Since A(r) £ 
L 1 (R + ) n Co(M + ), we will get convergence. To be more precise, if 

j R = sup |r r (0, t) |, a R = max\A(t)\ 

then 

oo 

1r <ur + jr / \A(s)\ds 



and ■jr — ► as i? — * oo. Since 1^(0, i) — * in L 1 (M + ), we have || || i,=^[JE,oo) — » 
as R — > oo. The same is true about T(t). Application of formulas (|235[) and (|236|) 
finishes the proof. □ 

Remarks and historical notes. For an excellent exposition of the proof for 
the Baxter theorem in OPUC case, see [55]. In our case, some modifications were 
needed. Apparently, the first proof of the Baxter theorem for continuous case was 
given in [52] • See also (20] [47] (one has to pay attention to some inaccuracies in 
statements regarding the regularity of coefficients generated by summable acceler- 
ants). 
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13. DlRAC SYSTEMS 

In this section, we relate Krein systems to the well-known object in mathematical 
physics: one-dimensional Dirac operator. Consider the Krein system, given by (|57p 
and assume some regularity conditions, e.g. a(r),b(r) G Lf oc (R + ). Let A G C and 



<p(r, A) = [P(?r, A) + P, (2r, A)] , 



exp(— iXr) 
2i 



[P(2r,A)-P*(2r,A)j 



These functions are of the exponential type r and are not from H 2 (<C + ) anymore. 
They should be regarded as analogs of trigonometric polynomials (or Laurent poly- 
nomials). For the free case, i.e. A(r) = 0, one has ip(r,X) = cos(rA), ip(r, A) = 
sin(rA). If A € M, 

<p(r, A) = Re £(r, A), ip(r, A) = Im £(r, A), £(r, A) = exp(-iAr)P(2r, A) (245) 

Define £(r,A) for r < by £(-r,A) = £(r, A). Let a(r) = 2Re A(2r), b(r) = 
2ImA(2r). Consider the following Dirac operator 



D 



' fl ' 




. h _ 







' fi ' 







(246) 



—b(r) d/dr — a(r) 
-d/dr — a(r) b(r) 

where the Hilbert space is f\,fa G L 2 (R + ) x L 2 (R + ) and operator is made self- 
adjoint by imposing condition /2(C)) = 0. Another way to write D is as follows 

d 
dr 



where potential Q{r) is 



and 



Q(r) 



—b(r) —a(r) 
—a(r) b(r) 



1 

-1 



(247) 
(248) 



This form of Dirac operator is called canonical. Any Dirac operator can be 
reduced to this form by a suitable change of variables [49], p. 48-50. Just like the 
Krein system, the Dirac operator in the canonical form has a lot of structure due 
to a special choice of potential. 

Rather than, say, Schrodinger operator, D can always be defined as the closure 
of the naturally chosen minimal operator (HO], Theorem 7.1, p. 493 or [77], p. 99). 
We will start with the following 



Lemma 13.1. Functions tp and ip are generalized eigenf unctions ofD, i.e. 



D 



,ip(0,X) = 1, ip(0, A) = 



(249) 



More generally, the fundamental solution Xd for the system J^6| ) can be expressed 
via the fundamental solution for Krein system in the following way 



X d (r, A) = exp(-iAr)f/ X(2r, \)Uq 
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Uq 



1 1 

—i i 



Proof. The proof is a straightforward calculation. □ 

Next, we will show that the system {E(x, A)} is an orthogonal system in L 2 (M, da). 
There are many ways to see that but we prefer an algebraic one, based on the proper 
factorization of certain integral operator. Just like in Section 3, we start with the 
following consideration. Let H(x) be an accelerant and H(x) £ Lf oc [0, oo). For any 
r > 0, consider the integral operator 3i r : 

r 

Krfix) = [ H(x - u)f{u)du (250) 



given in L 2 [—r,r]. Since "K r is the translation of 'Kir defined on [0, 2r], we have 
I + "K r > for any r > and the conditions of the Theorem 12.21 in section 2 are 
satisfied. Moreover, we can express one resolvent via the other one, i.e. 

t r (x,y) =T 2r (r + x,r + y), \x\,\y\<r (251) 

Let us consider some R > and the factorization I + ^Kr — (I + L)(I + VL) where 
the lower-diagonal (in the sense of Theorems 12.21 and 15. 2p operator £ has kernel 
L(x,y),\y\ <\x\ < R. Since % R - Hermitian, we get U = &*. If I + t = {1 + V-)" 1 
and / + U = (I + V+)-\ then 

Lemma 13.2. The following is true 

e*p(i\x) = £(x, A) + / L(x, u)E(u, X)du, x e [-R, R] (252) 



-|x| 

Proof. We have the following formula for E(x, A) if x > 

2x 

£(x, A) = P(2x, A) exp(— iXx) = I exp(2iAa;) — / T2 X (2x, s) exp(is\)ds I exp(— iXx) 



o 



(253) 



Doing the change of variables s — x = t in the integral, we obtain 

X 

£(x,A) = exp(iAx) - / T 2x (2x,x + t)exp(iXt)dt = (254) 



= exp(iAir) - J T 2x (x - t, 0) exp[iXt)dt (255) 

— X 

We also have 

£.{-x,X) = exp(-iXx) - / T 2x (0,x + i)exp(iXt)dt,x > (256) 
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Therefore, from (|25T|) and the Theorem El we get E(x,X) = (I + V_) exp(iXx). 
The Lemma then follow from I + L = (I + V_) _1 . □ 

Clearly, this Lemma is an analog of Lemma [3721 but for the different chain. Now, 
we are ready to relate Krein systems to Dirac operators. But first we need the 
following 

Definition 13.1. We say that a non-decreasing function tJd(A), A € K is the spec- 
tral measure for Dirac operator ()246|) . if the following is true (see [50], Chapter 8): 
for any f\ G L 2 (R + ) and /2 € L 2 (R + ), the operator 5F given by 

oo 

P7] (A) = J (/i(rMr, A) + / 2 (r)^(r, A))dr (257) 
o 

is unitary onto L 2 (R, dcr^). 

The next Theorem establishes a further link between the Krein systems and 
Dirac operators 

Theorem 13.1. The measure da^X) = 2do~(X) is the spectral measure for Dirac 
operator. Moreover, the mapping 

oo 

f(x) e L 2 (R) — > [W/] (A) = j f(x)£(x,\)dx 

— oo 

is unitary onto L 2 (R, da). 

Proof. Let us first show that W is an isometry map. Indeed, let f{x) 6 L 2 (—R, R). 
From ([50]) . we get 

oo 

((/+&*)/,/) = I 

— oo 

In the meantime, from Lemma 113.21 
it 

J J{t)eM^t)dt = ((I + i)E(t, X),f(t)) L 2 [ _ RtR] = (£(i, A), (I + i*)f(t)) L *[-R,R) 

-R 

Therefore, (|258|) gives 

2 

dcr(X) = ((/ + ^)- 1 (/ + ^.)(/ + U)- 1 g , 5 ) = ||<?|| 2 

where g = (I + £*)/ = (I + II)/. Since J + IX is invertible on L 2 [—R, R] and R was 
chosen arbitrarily, we learn that W is isometry. Now, let us show that W is also 
unitary. Indeed, for any R > and f(x) G L 2 [— R, R], we get (Lemma 113.21) 
R 

J f(x)exp(iXx)dx = ([(I + £)£](!, A), /(a:)) = (£(x, A), [(I + U)f](x)) 

-R 

Functions of that kind are dense in L 2 (R, da(X)) because the span of characteristic 
functions of the intervals {[a, b)} are dense and each of these characteristic functions 



/(*) exp(iXt)dt 



da(X) 



(258) 



g{t)E(t,X)dt 
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can be approximated by the Fourier transform of finitely supported L 2 (R) function 
due to the regularity condition ([8|). Therefore, the range of W is the whole of 
L 2 (R, dcr(X)). That means W is unitary. 

Now, we can conclude the proof of the Theorem. Take any /i(r) € L 2 (R + ) 
and / 2 (r) G L 2 (R+). Let fi(-x) = fx(x), f 2 (-x) = -f 2 (x), x > 0. Consider 
/(x) = — if2(x) on R. Function y> is even, ip is odd. So, 



[W/](A) = / f(x)£(x,X)dx = / (h - if 2 )(<p + i4)dx 



{h<P + f*l>)da: = 2[5f\ (A) 



That proves 3" is unitary mapping to L 2 (R, 2d<r(\)). Since the range of W is the 
whole L 2 (R, rfcr(A)), we can write W~ 1 g = f(x) = fi(x) — if 2 (x), where f\{x) — 
[f(x) + f(-x)}/2, f 2 {x) = -\f{x) - f(-x)]/(2i) and g is arbitrary from L 2 {R,da). 
Then, 25"/ = g and 5" is unitary. □ 

It is easy to show that the spectral measure for Dirac operator is uniquely defined 
(see Lemma ri7.10p in Appendix. 

As usual, the following representation can be easily obtained from the Theorem: 



f(x,X)ip(y,X) 



ip(x,X)ip(y,X) 
%/}(x,X)il>(y,X) 



d(Jd{X) = 



S(x-y) 
S(x-y) 



(259) 



and this identity should be understood in the weak-[L 2 (R + )] 2 sense (i.e. it is true 
after multiplication by L 2 functions and integration in x and y). 

In case dao(X) = dX/(2Tr) discussed above, one has a(r) — b(r) — 0, £(x, A) = 
exp(iAx), dad(X) = dX/ir. The map W is then the standard Fourier transform. 

The Dirac operator plays the role of the so-called CMV matrix for polynomials 
orthogonal on the unit circle. Many results about the Krein systems and functions 
P(r, A) can be viewed from that perspective. 

It is also quite helpful to introduce the auxiliary dissipative operator. For any 
R > 0, consider the operator Dr on L 2 [0,R] x L 2 [0,R] given by the differential 
system (|246j) and the boundary conditions: 



/ 2 (0) = 0,/ 1 ( J R)+i/ 2 ( J R)=0 



(260) 



The domain of definition for that operator consists in functions from W /1 ' 2 [0, R] x 
W 1,2 [0,R} satisfying (|260[) . It is an elementary calculation to show that T)r is 
dissipative since 

Im(D fl /,/)-2|/ 2 (i?)| 2 >0 

for all / in the domain of D. This operator has compact resolvent, an integral 
operator that can be written explicitly in terms of the solutions to the corresponding 
equation. The formula is as follows and can be easily checked 



{'DR-X)- 1 (Ji,f a ) t = 



<p(r,X) 
i>(r,X) 



{-Z 12 {a,X)f 1 (a) + Z u (a,X)f 2 (a))da 
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where 
and 



lll]x) ] J (^(s,A)/i( S ) - Z 21 (s,\)f 2 (s))ds 



( yi (r, A), 2, 2 (r, A))* = X d (r, X)X^(R, A)(l,«)* 



(261) 



p(r,A) yi(r,A) 
_ V>(^A) 2/2(7-, A) 

Since (2/1,2/2)' is solution to the Cauchy problem, it always exists. Therefore, the 
kernel of resolvent has a pole at point A if and only if 



= det 



ip(R,X) 1 
ip(R,X) i 



= iip(R, X) - ip(R, X) = i exp(-iXR)P(2R, X) 



Thus, the spectrum of this operator is discrete and coincides with the zeroes of 
P{2R 1 X). Since the dissipative operator has the spectrum in C + , this is another, 
operator-theoretic explanation to the fact that all zeroes of P(r, A) are in the upper 
half-plane. The zeroes of P*(2i?, A), on the other hand, are naturally characterized 
by the spectrum of the operator DJj, adjoint to Dr. The symmetry of zeroes for 
P(r, A) and P*(r, X) is now a consequence of a simple fact in operator theory. 

One can say even more, infact 



Lemma 13.3. The following representation is true 

R 

P*(2P, A) = exp - / [a(s) + ib(s)} exp(2iXs)ds 



det; 



(262) 



where Do,k denotes the operator with a(r) — b(r) — 0. The regularized determinant 
is understood as the regularized determinants of operator I + (Dq r — A) _1 Q (see 
[68], p. 106). 

Proof. A simple calculation shows that the spectrum of Dq r is empty (infact, is 
equal to infinity). Therefore, (Dq r — A) -1 is always well-defined. We have 



<Pt R -\)- 1 (Juh) t 



cos(rA) 
sin(rA) 



[i exp(iXs) fi(s) + exp(iXs)f2(s)]ds 



exp(zrA) 
— i exp(zrA) 



[i cos(sA)/i(s) + i sin(sA) f 2 {s)\ds 



(263) 



Simple calculations show that (Dq r — A) _1 Q € § 2 and the regularized determinant 
exists. 

Then, we use the following trick (see [55], p. 75). Introduce the so-called "cou- 
pling constant" fi£C and the potentials = [iQ. Then, consider the correspond- 
ing functions P* (2P, A, fi) and 

R 



/(A, m) = exp 



— /x / [a(s) + ib(s)] exp(2iXs)ds 



det 2 (/ + M(DS, fl - A)- X Q) 
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It is easy to see that these functions have the same zeroes. That follows from the 
properties of regularized determinants and relation between spectrum of Dq R and 
zeroes of P*(2R,X) discussed above. For fixed A, the function P* (2R, A, fi) is of 
exponential type in fi. Therefore, it can be factored 



P*(2P, A, fi) = exp(ci^ + c ) ]^[ (1 m/M") exp(/i//i„) 



(264) 



where the constants co,ci and zeroes \x n all depend on A. Since P*(2P, A,0) = 1, 
we get co = 0. Taking logarithm of both sides in (|264p . and comparing the Taylor 
coefficients in front of fx, we get 

2R R 



exp(i\s)A(s)ds 



\a(s) + ib(s)] exp(2iXs)ds 



On the other hand, for det2, we have the following factorization result (see [55] . 
Theorem 9.2, part (a)) 

oo 

det 2 (I + KV*q,r - A) _1 Q) = [J (1 fi/Hn) exp(/VM«) 

n=l 

Comparing these two expansions, we get the statement of the Lemma (take fj, = 
1). □ 

The determinantal representations are usually very useful in practice. They pro- 
vide the natural factorization of entire functions of interest. In case potential A(r) is 
small at infinity (say A(y) £ L p (R + ),p < 00), one can get asymptotical expansion 
of any order by using the further regularization of det 2 involving det3,det4, .... 
Formulas for the kernel K(x, y) show that it has discontinuity on the diagonal. 
Therefore, one could have used the Carleman-Hilbert determinant instead of det 2 
regularization. 

Since we have determinantal formula for P* (2P, A) for finite R, we might hope 
to get analogous result for H a (X) in case A(r) £ L 2 (R + ) by just taking R — * 00. 

Theorem 13.2. If A(r) £ L 2 (R+), then (D - A)" 1 ^ £ S 2 and 



II a (A) = exp 



a(s) + ib(s)] e-xjp(2i\s)ds 



det 2 ,AeC+ (265) 



Here Do denotes the free Dime operator, i.e. D with a{r) — b(r) = 0. 

Proof. The integral operator (Do — A) _1 Q has the following kernel Kq(x, y, A): 



-J Xy cos(Xx) [a(y) + ib(y)] e lXy cos(Xx) [b(y) - ia{y)} 



K (x,y,X) -- 
if y > x > and 
K a (x,y,X) = 



if < y < x. Since A(r) £ L 2 (K+), we have (D - A)" 1 ^ £ § 2 for any A £ C+ and 
the regularized determinant exists. Now, fix A £ C + . The function LT(2P, A) — » 



-e tXy sin(Ax) \a{y) + ib{y)] e lXy sin(Ax) [b(y) - ia(y)\ 



ie' lXx [cos(Xy)b(y) + sin(Xy)a(y)} ie lXx [sin(Xy)b(y) — cos(Xy)a(y)] 
-e lXx [cos(Xy)b(y) + sm(Xy)a(y)] e lXx [sm(Xy)b(y) — a{y) cos(Ay)] 
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lie (A) as R — * oo (see Theorem On the other hand, one can easily check 

that det 2 (/+ {T>* R - \)~ 1 Q) -► det 2 (I + Q(T> - \)~ 1 Q) as well. □ 

Now, let us study the wave operators for D. The following result establishes a 
connection between the stationary and non-stationary scattering approaches. 



Theorem 13.3. Ifa(x),b(x) £ L 2 



Q±f = lim e""~ i 

t — >=foo 



then the wave operators 

tV -itTtn 



10 f 



exist. The limit is understood in the strong sense, f — (/i, f%) t G [L 2 (K + )] Z 



Proof. The free evolution of Do is given in Lemma 117.111 from Appendix. It is 
actually a shift after some unitary transformations. Since each of the operators 
is unitary, it suffices to check the existence of strong limit for vectors 



/ = (/, 0)', where the scalar function f(x) G C 3O (IR + ). The existence of strong limit 
for vectors with zero as the first coordinate can be proved in the same way. Due to 
linearity, that is enough to conclude the convergence for all Cg 30 ^" 1 ") vector- valued 
functions that give rise to subspace dense in L 2 (R + ) x L 2 (M + ). From Lemma 117.111 
we have 

f{x +t ) + f{x-t) l i r f(\x-t\) 

-i(f(x-t)-f(x + t))\ 2[-if(\x-t\) 

where the last formula holds for t large enough because the support of / is finite 
(say, in the interval (0, a)). Consider the function 



-itD 



J 2 



i?(*,A) = — 



f(\x — t\)ip(x, X)dx — i J f(\x — t\)ip(x,X)dx 

.0 

To prove Theorem ll3.3[ it suffices to show that the limit of $(t, A) in L 2 ' 2a (R) exists 
as t — > oo. The following relations are true 



WI/.A) — / f(\x-t\y Xx P(2x,\)dx= i J f(\s\)e- tXs P42t + 2s,X)ds 

—a 

a 

(2t + 2s, X)ds = P* (2t + 2a, A) / /(r) cos(Ar)dr 



f(\r\)e 



dr 



+ J A(2t + 2s)P{2t + 2s,X) I f{\ T \)e~ lXT dTds. 

— a 

Since f(x) is smooth, 

s 

J f(\T\)e— x dr 

— a 

Due to Lemma 1531 and Theorem lll.il 



< 



C_ 

\/a t TT 



(266) 



P«{2t + 2a,\) J /(t) cos(Ar)rfr Il a {X + iO) XE -W J f(f) cos(Ar)dr, (267) 
o o 
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where the possible singular component of da is supported on the Borel set E, and 
Xe<= is the characteristic function of the complement to E. The convergence is 
understood in L 2,2<T (i?) sense. Here we also used Remark after the Lemma | 
Generalized Minkowski inequality and (|266[) yield 



A(2t + 2s)P{2t + 2s,X) / f{\T\)e~ lAT dTds 



< C 



\A(2t + 2s) |ds 



P(x,X)\' 
sup / ^— ; — da 

x >o J A 2 + 1 



2.CT 

1/2 



The second factor is bounded due to Lemma [8751 Function A(x) € L 2 (M. + ), 
fore the first factor tends to as t — * 00. 



there- 
□ 



Remark. We not only proved the existence of the wave operators, but also 
deduced the formula for them, the right-hand side of (|267p . Notice that this map is 
isometry. Part of the arguments above are well-known in the theory of polynomials 
orthogonal on the unit circle 26, 66 . 

Very interesting effect can be observed in the case da g(Szego), and A(r) — > at 
infinity in some sense (say, A(r) £ L p (R + ),p < 00). As was discussed before (see the 
paragraph after Theorem l8.2p , the limit P*(r n , A) is not necessarily uniquely defined 
and might depend upon the choice of the subsequence r n — > 00. One can easily 
show that the proof of the Theorem above can be adjusted to this situation with 
the exception that the limit lim^oo e** n2) e - ** n:Do depends upon the choice of time 
sequence t n and the limiting operators will actually differ only by the unimodular 
factor. However, due to Lemma 18.61 this phenomena cannot be observed for real- 
valued A{r) (i.e. when b(r) = 0). 

At this point, we need to mention that in OPUC theory, the free CMV matrix is 
unitarily equivalent to the shift in £ 2 (Z) space. Infact, the same is true about the 
evolution for free Dirac operator. Indeed, consider the following operator 

. d 



Do = Z- lr DoZ,Z 



1 
71 



,D 



dr 





. d 
dr 



The new domain of definition is [i/ 1 (R + )] 2 with additional condition /i(0) = 4/2(0). 
Now, if one maps the Hilbert space / = (/i, / 2 )* £ [L 2 (M+)] 2 to £ 2 (R) by 



f(r) — ► g(x) 



fi(x), x>0 
ih{-x), x < 



then the operator CD happens to be unitarily equivalent to the selfadjoint operator 

Lq = —i— on L 2 (IR) with domain of definition H 1 (M.). Its free evolution is just the 
dx 

shift: exp(it£jQ)g(x) — g(x +t),t £ R. Since we performed only the unitary trans- 
formations, the evolution exp(ztDo) (see Lemma 117.111 in Appendix) is unitarily 
equivalent to the shift in L 2 (R). Notice that potential is transformed to 



Z^QZ 





2A{2r) 



2A{2r) 




(268) 
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The unitary transformation of [L 2 (M + )] 2 to K. will result in nonlocal perturbation 
of £jq: 

[Hg](x) = -i%L+A(x)[&g](x) 
where Hermitian function A(x) and operator § are given by 



T/ \ J -2A(2x), x > c f \ i \ 
A (x)=*{ 2A(-2x), x<0 > § ^) = s(-z) 

The analysis of £> is non-trivial (rather than in the case § = 7) and is equivalent 
to analysis of the original Dirac operator T) or Krein system. Notice that in the 
Fourier space, this operator can be formally written as 

oo 

-A/(A)+ / V(X + t)f(t)dt 



with real- valued V(X) being the Fourier transform of A(x). The analysis we have 
done before implies the corresponding properties of this operator. 

Let us consider the scattering problem for Dirac operator and relate scattering 
parameters to the parameters of the corresponding Krein system. Consider, for 
simplicity, operator D with finitely supported coefficients a and b. Then, there 
is the so-called Jost solution F(r, A) : T>F = XF defined by the asymptotics at 
infinity: F(r, A) = (/i(r, A), f2(r, A)' = exp(iXr)(i, 1)* for r large enough. Let us 
introduce the scattering data for the Dirac operator: 

A d (X) = (/ 2 (0, A) - i/i(0, A))/2, B d (X) = (/ 2 (0, A) + i/ x (0, A))/2, 

Tk(A) = Af{X),R d (X) = B d (X)/A d (X) 

Coefficient T d (X) is called the transmission coefficient, R d (X) is the reflection coef- 
ficient, /2(C), A) is Jost function. These notations are quite natural. If one extends 
a and b to the negative half-line as zero, then 

F(r, A) - J 4 d (A)exp( i Ar)( i ,l) t + S d (A)exp(- i Ar)(-i,l) t ,r < 0, 

or 

T d (X)F{r, A) = exp(zAr)(z, l) 4 + R d (X) exp(-zAr)(-i, 1)*, r < 0, 
Lemma 13.4. The following relations are true 

A d (X)=Vl(X),B d (X) = <B(X),R d {X) = f(A), / 2 (0, A) = n Q (A), 

tt|/2(0,A)| 2 7r|n(A)| 2 
f(A)- Schur function of the Krein system, functions 21(A), 53(A), II Q ( A) are taken 
from consideration of A(r) f= L 2 (R + ) or A(r) € i 1 (R + ) cases. 

Proof. The proof is a straightforward calculation. □ 

The last formula in (|269[) is of great importance. It gives a factorization of the 
spectral measure density via some function analytic in the upper half-plain. Formu- 
las of that sort have analogs in the scattering problems for some PDE [T9l [18] . The 
equivalence of reflection coefficient from the scattering theory of quantum mechan- 
ics and Schur function is a remarkable fact, which, perhaps, was not completely 
understood and used by both mathematical physicists and analysts. 
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The general spectral theory allows to get new natural interpretation for various 
quantities considered before. One example is provided by the following Lemma. 

Lemma 13.5. For any Xq G C + . the operator Im(D — An) -1 has matrix-valued 
kernel G(x, y, An) and 



2 Im A 



\P*(2r, A) | s 



da(X) = Tr G(r,r, A ) 



|A-A | 2 

— OO 

Proof. The spectral representation for the resolvent yields (see (|259j> ) 

OO 

1 



A -An 



v(x,X)(p(y,X) ip(x,X)ip(y,X) 
^(x,X)ip(y,X) tp(x,X)ip(y,X) 



da d (X) = (D-Ao)- 1 ^) (270) 



Therefore, 

Tr lmG(x,y, A ) = 



Im A 



|A-A | 2 

Now, the Lemma is straightforward. 

This Lemma allows to control the integral 

" P,(r,X)\ 2 
A 2 + l 



[<p(x, X)(p(y, X) + ip(x, X)ip(y, A)] da d (X) 



□ 



da(X) 



by using the standard tools of, say, perturbation theory. In particular, if A — > in 
some sense, then this integral tends to 1/2, the value for unperturbed case. 
Recall that the CMV matrix corresponding to Verblunsky coefficients a n is 



e = 



* * * 

* * * 
* * * * 
* * * * 
* * 



A 
A x 



A k 



A = 



0-2kP2k-\ 
P2kP2k-l 



«0 

Po 



aiPo 
-Sido 



-Cl2ka2k-1 
-P2k0-2k~\ 

PlPtt 

-piao 



Cl2k+lP2k 
— 0,2k+lO,2k 



P2k+lP2k 
— p2k+\0-2k 



and pk = (1 — Iftfcl 2 ) 1 / 2 - To show how CMV matrix corresponds to Dirac operators, 
we prefer to write it in the equivalent way (by introducing two Hilbert spaces ^ 2 (Z + ) 
corresponding to even and odd indices): 



e = 



Gii C12 
C21 C22 
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with 



Cn = 



e 2 i 





PoPi 





—aia% 




















Po 


-picto 





-p2Ctl 























(02/03 



-0304 (04/05 





-P4»3 









-p 5 a 4 
* 




Poai 
















P2a 3 














P304 


P4fl 5 














Psa & 


* 














* 


* 



,e 2 



-a ai 
P1P2 











-a 2 a3 

p 3 p 4 -a 4 a 5 

p 5 a e 





It is well known that the formal discretization of the continuous Schrodinger 
operator produces a discrete Schrodinger operator, a particular case of the Jacobi 
matrix. For the Dirac operator, the situation is a little bit different because it 
is self-adjoint and the CMV matrix (an analog of Jacobi matrix in this case) is 
unitary. Notice that D is unitarily equivalent to 



-i— -2iA(2r) 
d 

2iA(2r) i — 
dr 



1 

71 



1 i 
1 -i 



(271) 



with an appropriate boundary condition at zero. Conjugating D by the matrix 



i 
-1 



we get that r 1 Dr has the same diagonal but the off-diagonal elements have the 
form (|268[) . Consider the following formal discretization of D 



A = -ih-\L-I) -2U(2r n ) 
h 2iA{2r n ) ih~ l {I-R) 

where R is the right shift, L - the left shift. Then, consider 



I + ihVf 



L 2hA(2r n ) 
-2hA{2r n ) R 

Ah) 



Taking the Verblunsky coefficients {a k } : a\„ = 0, a^n-i — 2hA(2r n ), we get 
/ + ihDh = CMV + d(h). Thus, formally, CMV matrices and discretization of 
Dirac operators are related via this very simple identity. Since L — I and I — R are 
of order h when acted on smooth functions, we could also say that 



exp 



L 

-2hA(2r n ) 



2hA(2r n ) 
R 



o(h) 



We do not pursue the goal of making any accurate statements regarding these 
discretizations but that can be done. 
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Remarks and historical notes. The one-to-one correspondence between 
Krein systems and Dirac operators was discovered by M.G. Krein in his semi- 
nal paper [44 . Unfortunately, no proofs were given. The determinantal formulas 
obtained in this section are new to the best of our knowledge. In the meantime, 
analogous results for differential equations were obtained earlier (e.g. [38]). These 
ideas were also used quite recently [41]. If a(x),b(x) £ i 1 (R + ), the existence of 
wave operators follows from trace-class perturbation argument |59j . In the case 
a(x),b(x) E L P (R + ), 1 < p < 2, the wave operators were studied by Christ and 
Kiselev [TU] . The analysis was based on establishing the asymptotics of generalized 
eigenfunctions (essentially, asymptotics of P*(r, A)) for Lebesgue almost any value 
of A. If p > 2, one can use results from [43] to construct examples with no abso- 
lutely continuous spectrum. Thus in this case the wave operators might not exist 
at all. The proof of Theorem 113.31 is taken from [17] . Independently, Barry Simon 
obtained analogous results for CMV matrices in the Szego case. 
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14. SCHRODINGER OPERATORS 



Let us consider Dirac operator (|246|) with b(r) = and absolutely continuous 
a(r). For the corresponding Krein system, we have A(r) G R and from Lemma I7T41 
we learn that if (a;) is real and continuous on R, the measure da is even. Operator 
D takes form 

d/dr — a 

—d/dr — a 

It has the following domain of definition {/i(r), /2M G £ 2 (R + ) x L 2 (R+)}, / X ( 2 ) are 
absolutely continuous, / 2 — a/2, /{ + a/i G L 2 (R + ), /2(0) = 0. Consider operator 

T> 2 = 



D 



0-Ci 
Jt 2 



where 



M i = -^2 + 9i. /i(0) + o(0)/i(0) = 



^ = -^ + 92, / 2 (0) = 



(272) 



and potentials are 

91 = a 2 — a', 52 = a 2 + a ' 

Obviously, Jfj and Jf 2 are different only by the order in factorization: JCj 
0*0, 3-f 2 = 00*, where = —d/dr — a is formal differential expression. 
The following is true 

d 2 <p 
dr 2 
d 2 ^ 
dr 2 



+ qi <p = AV p(0 I A) = iy(0,A)=-o(0M0 l A) 
gaV = A 2 ^, V>(0,A) = 0>'(0,A) = A 



That means if and ^ are generalized eigenfunctions for "K\ and IK2, respectively. 
Since D(<^(r, Q),ijj(r, 0))* = 0, we get (compare with Lemma r4.5|) 



ip(r, 0) = 0, (p(r, 0) = exp 



a(t)dt 



(273) 



Definition 1.5. The non-decreasing function p^)(E) is called a spectral mea- 
sure for the general Schrodinger operator OXi = —d 2 /dr 2 +q with mixed boundary 
condition /'(0) = hf(0) if the following is true. For any f(r) G L 2 (R + ), we have 



\f(r)\ 2 dr = 



f(r)<P(h)(r,E)dr 



dp(h)(E) 



where <P(M(r, B) is the generalized eigenfunction, i.e. for any E € R 
-<P(/0 + =Eip [h) , tp[ h )(0,E) = l,ip[ h) (Q,E) =h 
If = 0, we get the Neumann boundary condition. 



(274) 
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Definition 14.1. The non-decreasing function p/^^E) is called a spectral measure 
for the general Schrodinger operator %2 with the Dirichlet boundary condition if 
the following is true. For any f(r) G L 2 (R + ), we have 



oo oo 



\f(r)\ 2 dr 



/M<P(oo)(r-, E)dr 



dp(oo)(E) 



where <^(oo) ( r J E) is the generalized eigenfunction, i.e. 

-V('oo) + W(oo) = E( P(oo) , <P(oo) (0, E) = 0, <// (oo) (0, E) 



1 



The spectral measure for the Schrodinger operator with locally integrable po- 
tential always exists. But it is not necessarily unique [56l I49j . 

One can use Lemma 113.11 to prove the following Theorem. We consider the 
usual normalization of the measure da by saying that the function a is odd and 
a(X) = (a(A-0)+ a(X + 0))/2. 



Theorem 14.1. If p\{2) are spectral measures for the operators !Hi(2), then 

PiW = { VZI , P*W = < (275) 

0, A < 



Proof. Indeed, take any function f{x) € L 2 (R + ). Let /(— x) = f(x),x > and 
consider f(x) on the whole line. We have 

oo oo 

= J 8.(x,X)f{x)dx = 2 J <p(x,X)f{x)dx 

-oo 

Therefore, the first formula of (|275p is straightforward due to Lemma Tl3. II To get 
an expression for /02(A), one should take the odd continuation of function f(x). □ 

Notice that both /9i(2)(A) are constants for A < 0. That means UC 1 ( 2 ) are both 
nonnegative operators. That is not surprising since D 2 > and D 2 is decoupled 
into the direct sum of JCi and "Ki- Formula for /02(A) shows that not only P2W 
has no jump at (no zero eigenvalue for IK2), but also it decays at in a certain 
way. 

Consider the so-called "free" case, i.e. a(r) = 0. Then, A(r) = 0, da = dA/(27r), 
<7i = q 2 = 0, (p(x, A) = sin(Aa;), tjj(x, A) = cos(Ax). Moreover, pi(X) = 2\ 1 / 2 /n 
on K + is the standard spectral measure for the free Schrodinger operator on R + 
with the Neumann boundary condition and /52(A) = 2A 3//2 /(37r) on R + is the spec- 
tral measure for the free Schrodinger operator on R + with the Dirichlet boundary 
condition. 

Let us assume we are given nonnegative self-adjoint Schrodinger operator 3^.^ = 
—d 2 /dr 2 +q with mixed boundary condition /'(0) = hf(0) at zero and q G L 1 1 oc (K + ). 
Denote its spectral measure by dpi l (E). Then, there is the unique Dirac operator 
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and Krein system that generate this Schrodinger operator in a way described above. 
Indeed, define the functions 

a(r) = -^,4(r)=a(r/2)/2 (276) 

where ifi(r) is solution to the equation —ip" + qip = 0, ^'(0) = h, ip{0) = 1. We 
have ip"(r) € L 1 1 oc (IR + ) and ip(r) > by the oscillation theory for Sturm-Liouville 
operators ([77j. p. 218). Therefore, a(r) is absolutely continuous on R + and a(r) 
satisfies the Riccati equation 

q = a 2 -a' (277) 

Let A(r) — o(r/2)/2. Notice that A(r) is absolutely continuous. Krein system with 
coefficient A(r) generates "Kh- Let us assume that there are two different Ai(r) and 
A2(r) that generate the same 3~Ch- Then, due to Theorem 1 14. 11 measures o~\ and a% 
of these Krein systems are the same. Then, A\ — A-x- Notice that A(0) = —h/2. 
Thus, we proved 

Lemma 14.1. The operator 3~Ch is generated by the Krein systems if and only if 
Kh is nonnegative. The Krein system is unique and A{r) is given by \270\j . 

In general, Riccati equation (|277| has many solutions. For instance, if q = 0, 
the general solution is given by 

a h (r) = -h(l + rh)-\h > (278) 

But a,h(0) = —h so they can be all distinguished by the value at zero. If one 
considers q = and boundary condition f'(0) = hf(0), then the corresponding 
operator is nonnegative for h > only. In this case, ip( r ) = hr +1 and the formula 
(|276p gives exactly a^(r). 

The case of a Dirichlet boundary conditions is a bit subtle. Consider positive 
IKqo = —d 2 /dr 2 +q with Dirichlet boundary condition at zero and locally summable 
potential q. The problem here is that we don't know a(0), the initial condition for 
solving Riccati equation 

q = a 2 + a 1 

Interestingly enough, we might have many Krein systems that correspond to the 
same Schrodinger operator "Hoc. For instance, all Ah(r) — — a/ l (r/2)/2 (formula 
(|278|) ) generate the same Schrodinger operator with Dirichlet boundary condition 
and q = 0. Different Ah have different ah- Due to (|275[) . these measures are different 
by the jump at zero only. Let us consider cr(A) = <xo(A) + h8(X)/2, where 9(X) is 
Heaviside function, h/2 > is the jump at zero. Then, the corresponding H (x) = 
h/2 for all x. Given any r > 0, solution of equation (|35f is T r (s,t) = h(2 + hr)^ 1 . 
Then, A(r) = h{2 + hr)- 1 = A h (r). 

Thus, the natural questions are when is IKqo generated by Krein system and how 
to describe all Krein systems that give rise to JCqo? 



Lemma 14.2. The operator IHoo is generated by some Krein system if and only if 
there is some ip(r) > for r > such that —tp" + qip = 0. Moreover, if A(r) is 
coefficient of Krein systems generating JCqo , then 

AM = a(r/2)/2, where a(r) = ( 279 ) 

ip(r) 

with some ip( r ) satisfying the properties given above. 
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Proof. If there is some positive tp{r) satisfying equation, then the Krein system can 
be easily constructed by letting 

a(r) = ||^, A(r)=a(r/2)/2 

Conversely, assume that there is at least one A(r) generating 3ico. Then, the dual 
system with coefficient —A(r) will generate Dirac operator D_. Let ip_ (r, X),ip-(r, A) 
be the corresponding generalized eigenfunctions. Then, (|273|) yields 

a(r) = tff_(r,0)/<p-(r,0),A(r) = o(r/2)/2 
and ip(r, 0) satisfies —(p"(r, 0) + q(r)(p(r, 0) = 0. □ 



There are several different ways to reformulate this criteria. For instance, solu- 
tion i/)(r, 0) mentioned above exists if and only if the corresponding operator with 
the boundary condition f'(Q) = hof(0) is non-negative for some ho- Notice that if 
that is true for ho, then it must be true for any h > ho- One can also easily state 
this criteria in terms of the spectral measure p2- If er, obtained from the formula 
(|275p . generates a Krein systems, then this Krein system generates Moo- 

In view of these two Lemmas, one can suggest the following reduction of Schrodinger 
operator to the Krein system. Let "Kh be any Schrodinger operator bounded from 
below. Add large positive number 7 to "Kh so that it becomes strictly positive. 
This transformation simply moves the spectrum to the right and does not change 
the spectral types. For %h + 7, the Lemma Tl4. II is applicable. For Moo, the algo- 
rithm is the same but one has to apply Lemma 114.21 

Now, let us briefly discuss the solution to the inverse problem for Shrodinger 
operators. For 3~Ch bounded from below, the problem can be reduced to the inverse 
problem for Krein system which we know how to solve (just follow the construc- 
tion in the first Sections). Assume that we are given the spectral measure p of 
Schrodinger operator !H bounded from below. Assume also that the potential q we 
want to find is continuous. Then, the asymptotics of p at infinity [49] is 

p(X) = 2A 3/2 /(3tt) +o(A 1/2 ), fortheDirichletb.c. 
p(X) = 2X 1 / 2 /tt- h + o{l), for mixed b.c./'(0) = hf(0) 

From this asymptotics, we can find the corresponding boundary condition and apply 
one of the algorithms discussed above to find measure a for one of the Krein systems, 
generating !H. Once a is known, we can find A, then a, and, finally, q. Notice 
also that this method gives a one-to-one correspondence between all Schrodinger 
operators, bounded from below, and spectral measures that yield accelerant: H G 
C m+1 (R + ) iff the potential q <E C m (R + ), m is an integer. An accelerant H is a.c. 
onM+ iff qe Ll oc (R+). 

There is a direct way of solving the inverse spectral problem for Schrodinger 
operators. This method is due to Gelfand and Levitan [25] • Let us discuss this 
method and compare it to Krein's approach. Consider the operator (|274| with a 
continuous potential. Assume that its spectral measure pth) is given. Define 



/3(A) 



p ( , l) (A)-2A 1 /2/^ j A>0 
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and 

F(x) = lim / cos(A 1/2 a;)(i/3(A), x > 



It turns out that the limit exists and F(x) is continuously differentiable in x. Con- 
sider 

F { x,y) = nX + y)+ 2 F{X ~ V) 

and an integral equation 

X 

K(x,y)+F(x,y) + J K(x,t)F(t,y)dt = 0, < y < x < oo (280) 
o 

One can prove that the solution K(x,y) exists and is unique. Then, the following 
relations solve the inverse problem. 

h = K(0, 0) = -F(0, 0), q(x) = 2-^-Klx, x) (281) 

ax 

Now, let us make an assumption that the operator 5£ is nonncgative. Then, we 
can find the unique Krein system that generates "K. Take <r(A) = p(A 2 )/4. For an 
accelerant, we have the formal representation 

oo oo 

H(x) = J cos{\x)d (<r(A) - A/(2tt)) = i J cos{ v qix)d (p(jj,) - 2^ 1/2 /vr) = F(x)/2 

-oo 

One can check that 

- K(x, y) = Y 2x {x + y, 0) + T 2x {x - y, 0) (282) 
Indeed, from (f35]) , we have the following identities 

X X 

F 2 x{x + y,0)+ / H(y-u)T 2x (x + u,0)du+ / H{y + u)T 2 x(x- u,0)du = H{x + y) 



T 2 x(x-y,Q) + J H(-y-u)T 2x (x+u,0)du+ J H(u-y)T 2x (x-u,0)du = H(x-y) 

o o 
Since A(r) is real-valued, H(x) is a real-valued, even function. So, by adding the 
last two formulas, we get 

X 

[T 2x (x+y, 0)+T 2x (x-y, 0)]+ J [H(y -u) + H(y + u)} [T 2x (x -u,0) + T 2x (x + u, 0)] du 

o 

= H(x + y)+H(x-y) 

Therefore, we have ([282]) , and K(x,x) = -[T 2x (2x,0) + r 2a (0,0)]. Taking the 
derivative, we obtain 

2-^-Kix, x) = 2(-A'(2x) + 2A 2 (2x)) = a 2 - a' = q 
dx 

and 

/'(0)//(0) - -o(0) = -2A(Q) = K(0,0) 
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Thus, Krein's approach gives the same answer and these two methods are essentially 
identical. The difference is that imposing condition on ()280[) to have the unique 
solution is weaker than saying that F(x)/2 is an accelerant. That allowed authors 
of [25 to deal with a more general situation. Notice that the inverse problem for 
the Krein system is in fact the problem of the factorization for integral operators. 
Indeed, given measure a, we construct the accelerant. For any r > 0, operator 
1 + "K r > 0. Therefore, T r (x,y) exists. To find it, we need to factorize (1 + !K r ) _1 . 
Once we do that, A{r) = r r (0, r) = -V+(0, r) by 1(13]). 

Let us consider the scattering theory for Schrodinger operator 3-C2 defined by 
(|272[) . We assume that a has a compact support that belongs to, say, [0, R]. That 
means potential qi has a compact support too. Consider the Jost solution F(r, A) : 
'K2F = X 2 F defined by its asymptotics at infinity: F(r,X) — exp(iXr),r > R. Let 
us introduce the scattering data 

A S (X) = (F'(0, X) + iXF(0, X))/(2iX),B s (X) = {iXF(0, A) - F'(0, A))/(2tA), 

T S (X) = A S (X)-\R S (X) =B S (X)/A S (X) 

Function F(0, A) is called the Jost function for Schrodinger operator. The next 
Lemma relates scattering and spectral data for Schrodinger operator with Dirichlct 
boundary conditions to the corresponding parameters of Krein system. One should 
remember that we deal not with arbitrary Schrodinger but with the one generated 
by Krein system with compactly supported absolutely continuous coefficient. 

Lemma 14.3. The following relations hold true 

mx) = ^A) + ^mmim MX) = B(A) _ aMm^m, 

2iX 2iX 
Rs{X) = f(A) ~ 2a + am\%X)Y m A) = ^ s(A) + Bs[X) = m(A) + S(A) = H(A) ' 

Proof. The proof is a direct corollary from Lemma Tl 3. 41 Notice that if a(0) = 0, 
the data coincide with main parameters in Krein system. □ 

Simple calculations show that |^1 S (A)| 2 = 1 + \B S (X)\ 2 if A e R\{0}. That follows 
from the identity | 21(A) | 2 = 1+ | 03(A) | 2 which holds for any Ael. These formulas 
also show that the analytical properties of A S (X) and R S (X) are worse than of the 
analogous functions for Dirac operator or Krein system. For instance, j4 s (A) has 
pole at zero iff <x(0) 7^ 0. Moreover, 

-o(0) 



MX) - 5g> 



iX 



+ l+F(X) 



where F(X) = Il(A)n _1 (A) is Weyl-Titchmarsh function for Krein system. Thus, 
we see that A S (X) might also have zeroes in C + . These zeroes must be purely 
imaginary. Indeed, if Ao is such a zero, then Ag is eigenvalue for the Schrodinger 
operator considered on the whole line, a self adjoint operator. Therefore, A§ < 
and Ao is purely imaginary. 

Remarks and historical notes. 

For the first time, the inverse spectral problem for the Schrodinger operator was 
solved by Gelfand and Levitan [25]. Their approach is applicable to any operator, 
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not necessarily bounded from below. In the recent paper [67], Simon essentially 
introduced an "accelerant" directly for the Schrodinger operator. Different factor- 
izations of Schrodinger operators and applications were discussed in many papers 
(see, for instance, [13]). These methods allow one to insert eigenvalues below the 
essential spectrum, for instance. 
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15. Scattering theory for Krein systems 



In this section, we will consider the scattering theory for Krein systems and Dirac 
operators from slightly different perspective. The strategy is close to what is best 
known as approach by Marchenko and Agranovich to solution of inverse scattering 
problem for Sturm-Liouville operators 2 . We will try to emphasize the algebraic 
aspect of this argument, i.e. why Hankel operators appear and how their inversion 
is related to scattering data. Let us start with Krein systems. For simplicity, 
assume that the coefficient A(r) is finitely supported within, say, interval [0, R] and 
is continuous on this interval. Then, clearly, there is a unique solution X sc (r,\) 
such that X sc (r, A) = Xo(r, A), if r > R where 

s lAr 



X (r,X) 







1 



is the fundamental solution for A(r) = 0. This solution X sc (r, A) is normalized at 
infinity by its asymptotical behavior. We will study this solution and the scattering 
data it defines. Then, we will find its relation to spectral data and show how to 
solve an inverse scattering problem. This construction will be valid for more general 
case A(r) € In the meantime, let us assume first that A{r) is compactly 

supported and do some preliminary calculations. Clearly, if X sc is given by its value 
at r = R (i.e., X SC (R) — X Q (Rj), we might try to study it by solving the Krein 
system backwards, from R to 0. That is equivalent to dealing with "mirrored" 
coefficient A^(r) = A(R — r) ■ X[o,R]( r )- Therefore, we can study the following 
problem: consider the Krein system with coefficient A\(r), it will be later taken 
equal to —A^ R \r) but so far it is an arbitrary function with support within [0, R]. 
For this A\(r), consider the following solution 

-iXR q 



Y(r) = X(r) 







1 



where X is the fundamental solution for A\(r) normalized by X(0) = I. Notice 
that 

Xi(r,X) =Y(R-r,-X) (283) 
satisfies the following properties: X[(r, A) = V^ x _ A (R)yXi(r, A) and X\{R, A) = 

Xq(R, A). Thus, X±(r, A) is the scattering solution for —A[ R ' and so if one wants 
to study the scattering problem for A(r), we just need to study Y(r) associated to 
Ai(r) = -A( H \r) and then use the formula ([555]) . 

Now, let us start with obtaining the formulas for elements of the matrix Y(r, A). 
To simplify the calculations, take 



Y 1 (r,X) = Y(r,X) 



e iX(R-r) o 

1 



X(r) 







We have0 



2l(r,A) 

-i\r 



»(r,A) 



<B(r,A) 2t(r,A) 



Yi(r, A) = 

Let us write the following identity which follows from (|108[) 

-e iXr <B{r, A) ' 



(284) 



Yi(r,X) 



a(r, A) 



r «B(r,A) 2l(r,A) 



4 In the calculations below we assume that A is real. 
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or 



Yx(r,X) 



1 s r (A) 

S r (X) 1 



2t _1 (r,A) 



2l(r, A) 



(285) 



where s r (A) = - e - lXrs B{r, A) %L~ l (r, A). Let 

r 







Y 1 (r,X) = I+l K^(s,0) 
Q3(r, A) 



o 



o 



ds,K^(s,0) 



K[ r \ s ,0) Ki r \s,0) 



C r {x)e lAx dx,s r (X) 



C r (x+r)e l x dx,J r {x) — —C r (r—x)-X[o. r ]( x ) 



Subtract the identity matrix from both sides of (|285|) . take adjoint of the matrices 
on both sides, and act by the operator 

3>+ 
T_ 



from the left. We recall that f± are projections from L 2 (EL) onto H 2 (M.) and H 2 (I 
Then, we have the following integral equation for matrix K^(x,0), x S [0,oo) 



9r I 



^ (r) *(-,0) + 



J r (x) 
Jr(x) 







(286) 



and an operator 



[5rf](x) 



J r (x + u)f(u)du 



is acting in L 2 (R+). Notice that K^(x,0) = if x > r. Consider an operator 30 r ) 
in [L 2 (R+)} 2 such that 



o 5; 

Sr o 



(I + X (r >)=I 



Operator Sr is contractive Hankel operator since 23 (r, A) / 2l(r, A) is analytic contrac- 
tion. Therefore, 0C r exists, is self-adjoint, and has matrix- valued kernel K^ r \x,y) 
such that K^ r \x,Qi) is exactly the solution of (|286[) . Now, let us compare these 
equations for different r e [0,R]. Notice that 23(i?, A)/2l(i?, A) is the Schur func- 
tion for Ai(r). Denote 



C(x)=C R (x),J(x) = -C(R-x)- X [ ,R](x) 



(287) 



As it follows from (fT72"]) and (fT75|) . we have C r (x) = C(x) for x e [0,r]. There- 
fore, 

J r (x) = J(x + (R-r)) (288) 
andK (r) (ir,0) can be obtained by inverting the matrix-valued Hankel operator. On 
the other hand, it can be expressed through kernels T, T by formula (|284|) . We then 
have 



K^ r \x,0) 



T r (x, 0) 4- T r (x, 0) T r (r, x) — T r (r, x) 

r r (x,r) — r r (x,r) r r (o, x) + r r (o, x) 
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Theis formula allows to find A\(r) as off-diagonal elements in the matrix 

r r (o,o) + f r (o,o) 



tfM(0,0) = -i 



2A 1 (r) 

2A 1 (r) r r (o,o) + r r (o,o) 



Now, we are ready to translate these calculations to the original setting for coef- 
ficient A(r). Take Ai(r) = -A^ir) and the n use the formula (|283[) . We then 
have 

e i\r/2 o 

e"^/ 2 



X sc (r,X) = Y l (R-r,-X) 
For Z(r, Afl 

JXr/2 Q 



iAr/2 _ i\r/2 



Z(r,X) 



Z(r,X) 





e *Ar/2 o 

e ~ lAr / 2 



ft-r/2 



+ / X^- r )(s-r/2,0) 

r/2 
oo 

iAs 



e iXs o 

e" lAs 



ds 



+ / L^ 2 )( S ,r/2) 



e 



-iAs 



r/2 



where L^ r / 2 \x, y) = K^ R r *>(x — r/2,y — r/2); x, y> r/2 satisfies 



r/2 



Jn- r (x + U — r) 

Jri- r (x + u — r) 



L^ r ^(u,y)du 
— 0; x,y>r/2 



J R - r (x + y-r) 

J R - r (x + y-r) 
Now, let us obtain the formula for Jfj_ r (x + y — r) . Notice that the Schur function 
for the coefficient Ai(r) = -A^ R \r) is equal to —*8(R, —X)/QL*(R, -A) (Lemma 
1731) . If 

A) 



= / D(x)e tXx dx 

JUL 



A(R, A) 

then Jfl_ r (a; + y - r) = D(x + y) by (|2"57)l and (|2"ggfi Therefore, 7>/ 2 )(>, y) is an 
integral kernel of the operator £/ r / 2 ) in [L 2 [r/2, +oo)l 2 given by 

CD* 



7 



D r/2 



r/2 



and 



[Dr/a/Ka?) = / D(x + u)f(u)du 

r/2 



Since \ < S(R,X)/A(R,X)\ < f, ||D P || < f for any p > 0. Therefore, 



D 

r/2 



r/2 



is invertible. Also, D r/2 = if r > R so £( r / 2 ) = for r > i?. 



^The motivation to introduce Z(r,X) comes from the fact that it is the scattering solution for 
Dirac operator, that will be made clear later. 
6 Notice also that D(x) = for x > R. 
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Now, we want to take R — > oo. For this, we need some regularity at infinity, just 
like we needed regularity of, say accelerant near zero in the previous constructions. 
The simplest and quite natural class of functions to consider is L 1 (R + ) n Co(R + ). 

Theorem 15.1. Let A(r) G nC*o(R + ) be the coefficient in Krein's system. 

Then, there exists the scattering solution X sc (r, A) such that 

X sc (r,X) = X (r,X)+o(l) (289) 

as A S K, r — > oo. This solution can be obtained as follows 



X sc (r,X) = e^ r ^Z(r, X),Z(r, A) = 



e iAr/2 






e -iAr/2 



L {r/2 \s,r/2) 





— i\s 



r/2 



where L^ r / 2 \x,y);x,y > r/2 is the kernel of the operator £( r / 2 ) : 

D 



D 



r/2 



D 



r/2 



r/2 





,[©r/2/](x) 



[/ + s>r/ 2 ] (/+i: (r/2) )- 

arid D(x) is given by 

oo 

g(X)=m== J D(x)e iXx dx,D(x)eL 1 (R)nC ( 

— oo 

The coefficient A(r) can be found from the following identity 

A(r)=4 1 /2 \r/2,r/2) 



D(x+u)f(u)du 



r/2 



(290) 



Proof. As it follows from the section on the Baxter Theorem, Levy- Wiener The- 
orem, and Lemma [12721 D(x) e L 1 (R + ) n C (M + ). Also, D p is contraction in 
L 2 [p, oo). By the general theory of Hankel operators, D p is compact in any of 
L p [p, oo),oo > p > 1 and any eigenvalue in L x \p, oo) is also an eigenvalue in 
L 2 [p, oo). This is because -D(x) g L°°(R). In particular, that means I + D p is 
invertible in L^p, oo) and, therefore, L^ p '(x,p) 6 L x [p, oo) n Co(M + ). Moreover, 
||L( p )(a;, p)||ii[p !00 ) — > as p — > oo which can be checked by simple iterations. Thus, 
X sc (r, A) is well-defined and satisfies (|289[) . Now, let us show that it is actually a 
solution. 

Consider R- any positive number and let Aji(r) = A(r) ■ X[o,B]( r )- It has finite 
support and therefore the statement of the Theorem follows from the calculations 
given above. If Dr(x) is the corresponding function, then ||.Dr(x) — Z?(x)||i — > 
as R — ► oo. That follows from the proof of Baxter's Theorem. Consequently, 
Lft\x, r) — > L^ r )(x, r) and X^(r, A) — > A sc (r, A) as i? — * oo and r is fixed. On the 
other hand, 

Then, 



X*(r, A) = X(r, X)X~ 1 (R, X)X (R, A), r e [0, i2] 
X£(r,X)=X(r,X) j^^ V >re[0,fl] 



For fixed 



2l(i?, A) -SB(fl, A ) 
-®(.H,A) 2l(i?,A) 

X«(r,A)^X(r,A) 



21(A) -95(A) 



-33(A) 21(A) 
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as R. — ► oo . Therefore, 



X sc (r, A) = X(r, A) 



21(A) 
-18(A) 



-23(A ) 
21(A) 



Consequently, X sc (r, A) is a solution. Relation (|290|) is true for truncated A(r) and 
therefore holds after taking R — > oo. □ 



Notice that 



*sc(0,A) - 



21(A) 
-05(A) 



-»(A ) 
21(A) 



This matrix can be regarded as scattering data for Krein's system. On the other 
hand, function -D(x) can also be regarded as scattering data. What is the relation 
between these functions? If f(A) or 23(A) are given, then | 21(A) | can be found from 
|2l(A)|~ 2 = 1 - |f(A)| 2 or |2l(A)| 2 = 1 + |23(A)| 2 , respectively. Then, 21(A) can 
be found from |2l(A)| because it is outer. Function 23(A) can be obtained from 
23(A) = f(A) 21(A) and D(x) may be recovered by taking inverse Fourier transform 
of 23(A)/2l(A). The converse result can be obtained from the following calculations 
which is in the core of the method. Assume that A(r) S L 1 (R + ). Then, 



21(A) -23(A ) 
-23(A) 21(A) 



1 3(A) 
3(A) 1 



a _1 (A) 
2n x (A) 



(291) 



All elements in the matrices are from the Wiener algebra. Let 



CXJ CXJ 

w=i + /„(,)fta ( Ai=/«4 



/A.) 



dx 



o o 
Writing down the integral equation for coefficients and taking suitable projections, 
we get 

1 -D* 1 [ «0) 00 1 _ [ ^) 
-D 1 [ /?(•) a(-) J " [ D(x) 



(292) 



where D is the corresponding Hankel operator. Of course, this equation can be 
used to find 23(A) and 21(A). Consequently, if D(x) is generated by Krein system, 
then it defines this system uniquely. The converse is also true. 

Theorem 15.2. Assume that D(x) is a given function from L 1 (IR + ) n C*o(K + ). 
// the Hankel operator T> is contraction in L 2 (W + ), then there is a unique Krein 
system that generates this D(x). Moreover, the coefficient A(r) e L 1 (R + )nCo(K + ). 
Conversely, any A(r) € L 1 (K + ) n Co(R + ) generates D(x) with these properties. 

Proof. The second part follows from the arguments given above. Now, let us start 
with D(x). Consider {291]). Let 21(A) = 1 + (f>(X),g(X) = 5 i(A) + g 2 (X), 2T 1 (A) = 
1 + -0(A), where 0(A), g 1{2 ) (A), ip(X) € W+(R). Then, (p9T|) is equivalent to the 
algebraic system 

4> = i> + 23(5i + g 2 ), 23 = gi + g 2 + ~$(g x + g 2 ) 
That can be rewritten as a system for cj> and 23: 

4> = 3>+(®0i), 23 = g x + y+$gi) (293) 

with two more equations 



(12 



:,V ) = -23ff2 - 3 5 +(235T) 



(294) 
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that determine gi and ip consecutively. The system (|293[) is equivalent to (|292[) . 
Since D is contraction, (1292p has the unique solution. Now, that <j> and 23 are found, 
$2 and if) can be found by (1294p . Notice carefully that we do not know yet that 
■0(A) + 1 is indeed an inverse to 1 + 0(A). It will be clear in a second. We have 

Lemma 15.1. Matrix 



21(A) 93(A) 
03(A) 21(A) 

obtained in this way is J -unitary. 
Proof. We only need to show that | 21 1 2 = 1 



| 23 | 2 . The following is true 



I21I 2 = 1+123 



!23 



-H 2 O = ?+(| 23 | 2 -|0| 2 ) ^ = 3>+(23 23- 



Plug in the expressions for 23 and <f> from (|293|) into the last formula. Thus we just 
need to check that 

= 3>+(®<7i) + ?+(»?+ (.910)) - 3> + (0P+G;i»)) 
Due to the first identity in (I293|) . the last equality can be written as 

but this is always true since the left-hand side is 

T+^T+Cg^)) - ®O>_( 5l 0)) = 5>+( ff iW) 

and the right-hand side is 

3> + (#P+(<7i®)) = P+fe® - ^P-foiB)) = T+{giW) 

□ 

As a corollary, we get (1 +4>)(1 + tp) = 1. Indeed, we have 5 = 23 21 1 from (|293|) 
and (1294]) . Then, 1 + = 2t - 23 g = (| 21 1 2 - | 23 | 2 ) 2l _1 = (l + 0)" 1 by Lemma 
above. 

Now, we can say that the function f = 23 2P 1 is analytic contraction for which 
the conditions of the Baxter Theorem hold true. Therefore, it generates the Krein 
system with A(r) S L 1 (M + ). Moreover, a(x),(3(x) £ C (K+). Therefore, C(x) S 
Co(K + ). By Lemma Tl2. 21 we have A(r) is also from Co(M + ). Clearly, the function 
D(x) corresponds to this Krein system. 

□ 

The class of D(x) we considered was the simplest one. In principle, the Hankel 
operator 2) is bounded under much weaker conditions (e.g. BMOA space for gi, 
see [ST]). 

Let us make a remark regarding the Dirac operator. Notice that the matrix- 
function Z(r, A) introduced above satisfies the following properties 



-iJ 



dr 



-iA(r) 
iA(r) 



Z(r,X) 



£ i\r/2 







e -i\r/2 



6(1) 



Thus, Z(r, A) is the scattering solution for the Dirac operator (written in a slightly 
different way, see (|271|l ). All results from the section can be easily translated to 
this case. 



CONTINUOUS ANALOGS OF POLYNOMIALS 



105 



Remarks and historical notes. 

The main results of this section are not new. They are essentially contained in 
the papers by Krein and Melik-Adamyan [35 11521155] , The authors were motivated 
by certain problem of continuation related to Hankel operators (and, essentially, 
coming from the complete solution of Nehari problem, see [I] [57]). Also, Melik- 
Adamyan considers slightly different canonical system. Using the approach of this 
section, one can provide another proof for continuous analog of Baxter's Theorem. 
We tried to make the argument almost purely algebraic. 

Also, one can use upper (lower)-triangular factorization of operator J) r to rep- 
resent its determinant via certain integral of A(r). Analogous calculations will be 
done in the next section for truncated Wiener-Hopf operators. This is a way to 
obtain formula similar to the so-called Borodin-Okounkov identity [9]. 
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16. Truncated Wiener-Hopf operators. The Strong Szego Theorem 



Consider the continuous accelerant H . In this section, we obtain an important 
formula for the Fredholm determinant of / + "K r in terms of the coefficient A(r) of 
the associated Krein system and prove continuous analog of the so-called Strong 
Szego Theorem [66]. We start with the following well-known result [5]. We omit 
the proof which can be obtained, e.g., by upper (lower)-triangular factorization of 
1 + JC r . 

Theorem 16.1. If H(x) is a continuous accelerant on R, then 



det(l + JC r ) = exp 



r u (0,0)du 



Lo 



Now, assume that we are given a continuous accelerant H(x). We have a relation 



and 



-fr u (o,o) = -|A( u )| 2 

au 



r 

r r (0,0) = H(Q) - J \A(s)\ 2 ds 



det(l + % r ) = exp[F(0)r] exp 



(r-s)\A{s)\ 2 ds 



and 



Notice, that 



det 2(1 + M r ) = exp 



- J (r-s)\A{s)\ 2 ds 




(295) 



(296) 



d_ 

rir 



r t 

lndet(l + 5£ r ) = fl'(0)- J \A{s)\ 2 ds, J- In det 2 (1 + J{ r ) = - j \A{s)\ 2 ds 



The fact that the both sides have limits for A(r) £ L 2 (M. + ) can be regarded as the 
weak Szego theorem. Notice that an approximation argument yields (|296[) without 
continuity assumption on H{x). The regularity H(x) € Lj 2 oc (K) is enough and the 
weak Szego theorem reads: In det 2(1 + JK r ) is nonincreasing and 

00 

Urn A lndet a (l + 5Cr) = - f \A(s)\ 2 ds 

r~ too dr J 



The strong Szego asymptotics is as follows 

det(J + 5C r ) det 2 (/ + ^r) 



exp 



r(H(0)- J \A(s)\ 2 ds) exp -r J \A(s)\ 2 ds 



exp 



s\A(s)\ 2 ds 
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and the limit of left-hand side exists iff 

oo 

G = J r\A{r)\ 2 dr < oo 
o 

Notice that under this condition one has 

r oo 

r J \A(s)\ 2 ds = r J \A{s)\ 2 ds + o(l) 
o o 
and assuming that A(r) € L 2 (R + ), we get 



(297) 



T 

J. r 



det(/ + J£ r ) 



det 2 (I + 0< r ) 



exp 



r(ff(0)- f \A(s)\ 2 ds) 



exp 



-r J \A{s)\ z ds 
o 



exp / v r (s)s\A(s)\ ds 



(298) 



where v r (x) = 1 for x £ [0,7*] and v r (x) = rx^ 1 for a; > r. Since v r (x) is non- 
negative, monotone in r, and v r {x) < 1, the left-hand side is also increasing and 
has a finite limit iff (|297p holds. The following result also gives a characterization 
of this case in terms of the spectral measure. That is also a continuous analog of 
Ibragimov's and Golinskii-Ibragimov's Theorems [27l [35l [66] . Some parts of the 
arguments below are borrowed from the discrete case [66] . 

Theorem 16.2. Assume A(r) £ L 2 (R + ). Then, the following statements are 
equivalent 

(i) G < oo 

(ii) The measure da is purely a.c. and 



ln(27rcr'(A)) = / l(x)aq>(i\x)dx € H 1/2 ( 



(299) 



(hi) T r is bounded 
Moreover, 



G = T^ = L = I 



wher$\ 



L 



x\l{x)\ 2 dx, I 



^aIUA) 



c+ 



n Q (A) 



d z \ 



where 1I Q (A) is given by (|206D . 



(300) 



Proof. The equivalence of (i) and (iii) as well as G — Too follow from the argument 
above (see (|298jl ). The rest of the proof is divided into several Lemmas. 



Lemma 16.1. If G < oo, then da s = 0. 



7 Recall that if A(r) £ L 2 (R), then In [27rcr'(A)] e L L (R) + L 2 (R) (Corollary [TOt . Also, 
notice that condition A{r) £ L 2 (R) can be expressed purely in spectral terms, i.e. through da. 
In particular, the Theorem 111.51 and Lemma ll7,12l from Appendix imply that conditions da s = 
and ln(27rcr'(A)) G # 1/2 (R) guarantee A(r) G L 2 (K+). 
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Proof. For real A, 



|P*(r,A)| >exp 



- \A(s)\ds 



> exp[— CVln r] 



That follows from Lemma l4~6l From (I224[) . we get 

oo oo 

|n(r,A)| 2 



Therefore, 



Thus da s = 0. 



/ ^V + l <cj \A(s)\ 2 ds < CV" 1 



da s (X) ^ Cexp(CVm7) 

—3 < > 0, as r — > 00 

A 2 + 1 ~ r 



□ 



Lemma 16.2. Assume g(X) is outer in N(C + ), g(ioo) = 1, and In |<?(A)| G L 2 (R) + 
Z-^R). J7ien, PTT]) below holds. 

Proof. Using multiplicative representation for outer functions and normalization at 
zoo, we have 



lng(A) = 2 / l(x) exp(iXx)dx 



with l(x) G L 2 (R + ) + W{\ 



Simple calculations show that 



xe- 2tx \l{x)\ 2 dx = i^tt)" 1 / |a A ln.g(A)| 2 d 2 A = {An)- 1 

Im A>e Im A>e 



dxg(X) 



5(A) 



Therefore, we always have 

x\l(x)\ 2 dx = (Att)- 1 



dxg(X) 



.9(A) 



d 2 X 



even though the both quantities can be infinite. 



d 2 X 



(301) 



□ 



For A(r) G L 2 (R+), we have ]n(2na'(\)) G L 2 (K) +i 1 (R), U a (ioo) = 1 and the 
Lemma can be applied to g(X) = LI" 2 (A) and thus we get L = I in (|300[) . 
We need the following auxiliary 



Lemma 16.3. For any R, we have the following inequality 

R 



r\A(r)\ 2 dr > n^ 1 



d x P*(R, A) 



P*{R, A) 



d 2 X 



(302) 



Proof. For any T > 0, consider n T = {X G C+, |A| < T}. Let us show that 

R 



r\A(r)\ 2 dr > n~ 



d x P*{R, A) 



P*(R, A) 



d 2 A 



(303) 
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for any T > 0. Then the general statement follows upon taking T — > oo. Fix any 
T. Then, it is sufficient to prove (|303|) assuming that A(r) £ C[0,R]. Indeed, any 
A(r) G L 2 [0, R] can be approximated by continuous functions in L 2 [0, R] norm and 
the both sides of (]303[) are continuous in A{r) with respect to L 2 [0, R] metric. 

Then, we just need to use the suitable formula from the discrete case and an 
approximation result given by Corollary 19.11 Consider large n, the discretization 
step h — R/n, Verblunsky parameters given by (|175|) . and the corresponding monic 
orthogonal polynomials Pk(z) and P£(z). Then, we have the following formula 
([55], Theorem 2.1.4) 



In 



IK 1 

3=0 



With our choice of Verblunsky parameters, 



In 



IK 1 

3=0 



-3-1 



r\A{r)\ 2 dr, n 



d 2 z 



(304) 



oo 



by the Riemann sum approximation. Over fi^, we have P*(e tXh ) 

ihe iXh d z P*{e iXh ) -+ d x P*(R,\) 
Since e lXh — > 1 uniformly over f^T, we also have 

ihd z PZ(e iXh )^d x P*(R,\) 



(R, A) and 



(305) 



Therefore, making the change of variables and using (|304p and (|305[) . we get (130311 . 

□ 



For A(r) e L 2 (M + ), we have P*(R, A) — > IL a {X) uniformly in {Im A > 6} for any 
5 > 0. Therefore, we always have an estimate G > I and consequently (i) implies 
(ii) due to Lemmas 116. 1[ 116.21 

Now, we are left with proving that G < I. Let us assume we have purely a.c. 
measure with density satisfying ln(27rcr / (A)) G iJ 1 / 2 (M) so that L < oo for the 
corresponding function l(x). We need to show G < L for the associated Krein 
system. 

We know that A(r) G L 2 (M+). Then, 2tt(j'{\) - 1 = ^(x) with ^(x) G £ 2 (R) by 
Lemma [T7T21 Formulas ([2"2T)1) and P2T|l imply 

H (x) = 2tt^(x) (306) 

where t/» is Fourier transform of ip. 

Let us take Hermitian Ir(x) such that Zij(a;) is continuous with compact support 
within [— R, R] and 

oo 
— oo 

as i? — > oo. Formula (|306|) . Lemma [17.12[ and Lemma [5.21 show that the corre- 
sponding Aft(r) — > A(r) in L 2 [0,T] for any T > 0. For each R, we can apply the 
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Theorem 19.41 For the corresponding sequence of Verblunsky parameters, we have 
[66] . Chapter 6: 



IK* 

3=0 



,(n)|2s-i'-l 



exp 



(")|2 



.fc=l 



(307) 



where 



Clearly, 



ln/4(0) = £ L 



(n) ikS 
k e 



k= — oo 



oo " 

fe=l ^ 



exp 



as n — ► oo. On the other hand, for any 6 > 0, we have 

R 

II (l-lafl 2 )- 5 '- 1 

<5<j7i<n 

as it follows from the Theorem 19.41 Therefore, 

R R 

r\A R (r)\ 2 < [ x\l R (x)\ 2 dx 



r\A R {r)\ 2 dr 



o 



and then G < I since A R {r) — * A(r) in i! 2 oc ( 
finished. 



"). So, G = I and the proof is 

□ 



The proof we used essentially utilized the strong Szego Theorem for Toeplitz 
matrices and the approximation of continuous orthogonal system by the sequence of 
discrete ones. In the meantime, one could have adjusted the various proofs directly 
to continuous case (see, e.g. [9] for continuous analog of Borodin-Okounkov identity 
which can probably be used for this purpose). Notice that if dak, k — 1,2 satisfy 
conditions of Theorem 116.21 then da = [eri] 7 [o^] 1 7 dX, 7 G [0, 1] also satisfies these 
conditions. 

Formula (|295[) is very important for many applications. In the theory of random 
matrices, one needs to calculate asymptotics of Fredholm determinants for some 
specific accelerants. Assume that coefficient A, corresponding to a given accelerant 
H, tends to zero fast enough. Then, solving the inverse scattering problem by 
methods of the last Section, one can obtain an asymptotics of A(r) at infinity. 
Assume that A(r) decays at infinity fast enough such that 



A(r) = 



In . C N+ i(r) 



„N+1 ' 



n— 2 



holds with constants j n and CV+i(0 G L°°(R + ), N is arbitrary. Formula (|295[) 
and 

r 00 00 00 00 

J(r- s)\A{s)\ 2 ds = r J \A{s)\ 2 ds-r J \A{s)\ 2 ds- J s\A{s)\ 2 ds + J s\A(s)\ 2 ds 

r r 
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oo oo 

shows that as long as all j n , J \ A(s)\ 2 ds, and / s|^4(s)| 2 <is are known, we can 

o o 

compute complete asymptotic of dct(l + "K r ) as r — > oo. But these two integrals 
can be explicitly expressed via the spectral data. The idea of using the inverse 
scattering theory to compute asymptotics of Fredholm determinants was pioneered 
by Dyson 22J. 



Remarks and historical notes. 

Various generalizations of the strong Szego formula to continuous case were ob- 
tained in 5, 591170] . Our version is optimal and new to our knowledge. On applica- 
tion of inverse scattering to random matrices see [SSJ . In [75] , the theory of Krein 
systems was used to study an asymptotics of the certain Toeplitz determinants. 
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17. Appendix 

In this Appendix, we collected the general results that we used in the main text. 

17.1. For section [2l The proofs of the following two Lemmas are given in [69] . 
p. 71 and p.99. 

Lemma 17.1. IfY r (x,y) is the resolvent kernel for K(x,y), which is continuous 
on [0,r] 2 (see Section 2), then 

S r (x,y) 



where 



T r (x,y) 



1 f [ K x & ... & 

nlj —J \y ti ... £» 

o o 



d£i . . . di n 



■\ i / k { | Wi+^y Jk( i 1 1 



u 



■kl-l K (a I ::: 1 >" £ l - 







I 1 6 ^ ) =det{J)T(e fc , »»)},! <fc,Z< 



»?i m ■■• »7i 

T/ie series converges absolutely. 

The next Lemma gives Carleman-Hilbert determinantal representation for resol- 
vent. We recommend an excellent book [3T] (Theorem 2.2, p. 207) for the modern 
presentation of that subject. It also contains the discussion of integral operators 
with discontinuity on the diagonal. 

Lemma 17.2. IfT r {x,y) is the resolvent kernel for K(x,y) £ C*([0,r] 2 ) (see Sec- 
tion 2), then 

r,(x,y) = ^f^G(7([0,r] 2 ) 
o r 

where 8 r (x,y) and S r are defined as S r (x,y) and 5 r given above, but relative to the 
modified kernel K(x,y) = K(x,y) if x ^ y and K(x,y) — on the diagonal. 

17.2. For section [4J We begin with some simple facts about the linear spaces 
with indefinite metric. Let [x,y] = (Jar, y) be indefinite inner product (in C 2 ). 
For any matrix M, introduce M c — JM*J. Then, [Mx,y] — [x,M c y\. Clearly, 
(AB) C = B C A C , (A )- 1 = {A- 1 ) if A is invertible. The following result is well- 
known [33 \8\ 

Lemma 17.3. If M is J -unitary, then |detM| = 1, andM~ x ,M* are J -unitary 
too. If M is J- contraction then M* is J- contraction also. 
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Proof. Taking determinant of M*JM = J, we get |detM| = 1. It is straightfor- 
ward that M is J- unitary if and only if M^ 1 is J-unitary. 
Then, clearly, M is J-unitary if and only if 



but that means 



so 



M*JMJ = I 



M*J = (MJ)- 1 



MJM*J = I 



which is the same as saying that M* is J- unitary. 

Assume that M is J-contraction and 1 is not an eigenvalue. Then, we have a 
general algebraic formula 

(M— /) ~ 1 (MM C — I) (M c —I) ~ 1 = {M C -I)- 1 (M C M-I)(M-I)- 1 = I+{M-I)- 1 +{M C -I)~ 1 

which can be rewritten as 

MJM* - J = JQ*(M*JM - J)QJ, Q = (M- J) _1 (M C - 7) 

Since Q is invertible, we have M*JM < J iff MJM* < J. In other words, M is J- 
contractive iff M* is J-contractive. If 1 is an eigenvalue, multiply M by unimodular 
scalar factor such that 1 is not an eigenvalue of the resulting operator and apply 
the argument above. □ 

Lemma 17.4. Assume that 

r 

p(X) — 1+1 f(x) exp(i\x)dx 



where f(x) € L 2 [0,r] and p(X) ^ for A S M. Then, p(X) is uniquely determined 
by 

1 



-r,r] 



i 



Proof. By Levy- Wiener theorem, 

oo 

. , . .„ — 1 = / h(x) exp(i\x)dx 

\pW\ J 



(308) 



where h(x) € L 1 (R) and h(—x) = h(x). We know h(x) on [— r, r] and need to find 
f(x) on [0, r]. As in Lemma [3~8l 



p{\) 1+ / h(s) exp(i\s)ds - 1 



which gives 



f(x) + h(x) + h(x- t)f(t)dt = 0,0<x<r 
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This equation determines f(x) on [0, r] uniquely from h(x) on [— r, r] since I+Oi > 
where 



•Hg(x) = h(x- t)g{t)dt 



The positivity oil + "K follows from (|308f . □ 
17.3. For section [H 

Theorem 17.1. Let dp,{\) be finite ncmnegative measure defined on the whole 
line. Consider the linear manifold L consisting of the finite linear combinations of 
exponents exp(iAr) with r > 0. Then, L is not dense in L 2 {dp) iff 

l ^-dX > -oo (309) 
1 + A 2 

Proof. Denote the closure of L in L 2 (dp.) by L^. Consider the function 

oo 

/(A) = / exp(— x + i\x)dx 



i + A 



o 



We can find sequence {/ n } in L which is uniformly bounded on R in L°° norm and 
/„ — ► / in the uniform norm on any fixed compact in R. For instance, 



„ 2 -i ^ +1 j n 
/™(A) = ^2 exp(iAj/n) / exp(-x)da 



Therefore, / G and g = (A — i)(X + i)^ 1 = 1 — 2/ S £ M . Let us show that L M 
is invariant under the multiplication by g. Because g is an elementary Blaschke 
factor, \g\ = 1 on the real line. Therefore, we only need to show that gip £ L M for 
any t/j £ L. But that clearly follows from 

(i) ij) is finite linear combination of exponents exp(iAr) for different r > 0, 

(ii) space is invariant under the multiplication by exp(iAr) for any r > 0, 

(iii) g itself belongs to . 

Now, once we know that is invariant under the multiplication by g, we know 
that g n £ L M for all n £ Z + . Let us consider the standard conformal map of A £ C + 
onto the unit disc: w £ J},w — (X ~ i)(X + i)^ 1 . Under this map, measure d\i goes 
into a new finite measure dr on the unit circle, which generates a new Hilbert space 
L 2 (dr). The subspace goes into the subspace L r . Moreover, since all g n £ L^, 
w" £ L T for any n £ Z, w £ T. Let us consider the subspace Z T obtained by closure 
of Z = Spanjl, w, ... , w n , . . .} in the L 2 {dr). Clearly Z T C L T . Let us show that 
actually Z T = L T . To do that, it is enough to prove that any function 



h(w) = exp 



w + 1 



,r > 0,w £ T 



1 

(which is an image of exp(iAr) under the conformal map) can be approximated by 
analytic polynomials in w in the L 2 (dr) metric. Consider functions 

pw + 1 
pw — 1 



h p {w) = exp 



,W £ T 
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for p < 1. Since d/j, is finite on R, we have r(— e, e) — * as e — > 0. In other words, 
r has no mass point at w = 1. Function /i(it>) is bounded in D and is continuous 
there except for the point k; = 1. Therefore, h p — > h in L?{drr) as p — ► 1. At the 
same time, each /i p is analytic in small neighborhood of D and therefore can be 
approximated by polynomials uniformly on the unit circle. Thus, L T = Z T . The 
Szego theorem ([55], Chapter 2) says that Z T is not dense in L 2 (dr) iff 



7T 

/" In r'{9)d9 > -oo 



Clearly the last condition is equivalent to (|309|) . 
Notice that due to 

oo 

/ dfi(X) < oo 

— oo 

we have that (|309[) is equivalent to 



□ 



In" //(A) 
1 + A 2 



dX > — oo 



(310) 



This Theorem has interpretation in the theory of Gaussian stationary processes 
with continuous time and that is very useful point of view on the whole theory 
of Krein systems. In the meantime, there are the so-called Krein strings [21j . 
differentiable operators more suitable to deal with stationary processes. 

Theorem 17.2. Assume that da is a measure on the real-line such that 

° rfcr(A) 



1 + A 2 



< oo 



Consider the linear manifold X of functions 



/(A) = / exp(i\x)f(x)dx, < n < r 2 
o 

where f(x) € C 1 [ri,r2] and is zero outside [ri,ra] C [0, oo). Then, X is not dense 
in L 2 (da) iff 



lnV(A) 
1 + A 2 



dX > — oo 



Moreover, let Xq G C + . Then 
1 



Dist 



A-Ao 



1 



L 2 (da) 



V2ImA, 



: exp 



ImAo f ln(27rcr'(A)) 
~2^~ J |A-A | 2 



dX 



(3H) 



(312) 



Proof. Consider a new measure dp: = da/(l + A 2 ) which is finite on the real line. 
Denote by Y the linear manifold of functions of the following form (A + i)f(X), f £ 
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X. Let Y/t be the closure of Y in L 2 (d[i). We only need to show that Y^ ^ L 2 (dfi) 



iff 



ln^(A) 
1 + A 2 



d\ > — co 



(313) 



Let be the space of functions from the proof of the Theorem 117. 1[ i.e. the 
closure in L 2 {d^i) of finite linear combinations of exponents exp(iAr),r > 0. It is 
not difficult to show that exp(iAr) 6 Y^ for any r > 0. That follows from the 
representation 



exp(iAr) = — i exp(r)(A + i) / exp(— x) exp(iXx)dx 



So, CYfj,. At the same time, each function 



(A + i) I exp(iXx)f(x)dx — i(f(ri) exp(iAri) — /(V2) exp(iAr2) 



+ / exp(iXx)[f'(x) + f(x)]dx 



can be approximated in L 2 (d[i) by the finite linear combinations of exponents 
exp(iAr). One should replace the integral by the Riemann sum and use conti- 
nuity of the functions /, /' to estimate the error. Thus — Y" M and one can use 
Theorem 117.11 to finish the proof of the first statement of the Theorem. 

Now, let us obtain the formula for the distance. For simplicity, consider Ao = i. 
The general case can be treated in the same way. We have 



inf 

fex 



1 



da = inf 



l-^(A + i)/(A) 



da 



1 + X 2 



inf 



dfi(X) = inf 



dfi(X) 



inf / II — wv(w)\ drtw) 



where the measure dr(w) was obtained from d/z(A) by mapping C + onto B via 
w = (A — i)(X + l) ■ Here we also used an approximation result from the proof of 
Theorem ll7.ll For the last inf, we can use the Szego formula [66], i.e. 



Dist(l,wZ T ) L2(dT) = exp 
and the proof is finished. 



2;r 



— f ln(2nT'(6))d9 
4?r J 



V2 



exp 



1 f ln(27rcr'(A)) 



2tt 



1 + A 2 



dX 



□ 
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17.4. For section [H 

Lemma 17.5. IfC(x) £ L 1 [0,i?], C(x) is continuous at zero and 

R 

fW — / C(x) exp(iXx)dx 



then 



1 



v 



lim - / sf(is)ds = (7(0) 

y->oo y J 


Proof. The proof follows from the standard estimates: 



[CO) - C(0)}exp(-sx)dx 



v y 
lim - / sf(is)ds = C*(0) + lim - [ s 

y^°° y J y J 



because the second term before the lim can be bounded by 



ds = (7(0) 



(7 



\C\\ 1 + \C(0)\ 

y 



V 

J s exp(~sS)ds 



□ 



where ws(C) — sup^grg^ \C{x) — (7(0) | — > as S — > 

17.5. For section 1101 The following Lemma controls the zeroes of the continuous 
orthogonal polynomial 

Lemma 17.6. Let 

r 



p(A) = 1 - / 7(x)e~ lAx dx 
o 

where j(x) £ (7 [0,r], 7(7') 7^ 0. //A n are zeroes ofp(X) and |Ai| < |Aa| < iften 
A n = A^ + o(l), n — > 00, where A„ =i„ + iy n and 



^ + 2/^ = ItMI 2 exp(2rj/„), x„ = r 1 [tt/2 + 7m + <p] , n £ Z, 



(314) 



Here, j(r) = \j{r)\e lv . 

Proof. From Lebesgue-Riemann Lemma, A„ £ C + for large n. Integrating by parts, 
we get 



1 



1 



1 



p(X) = 1 + — [i(r)e~ iXr - 7(0)] - ^ [i{r)e~^ - 7 '(0)] + - / j"(s) e - iXs ds 



iX l ' v '~ /x ~ ;i A 2 

Therefore, the equation p(X) = can be rewritten as 

Y(r) 1 



zA 



1 , 7(0) f (0) 
zA A 2 



7(r) 



1 



7(r) 



OdAp 1 ) 



7"(s)e jA(r - s) ds 



By Rouche's Theorem, A„ will be approaching the roots A„ of equation 

e~ lXr 1 



iX 
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Then, the first equation in (|314j) easily follows upon taking the absolute value 
squared. The second one can be obtained by taking the real part of identity 
ry(r}e~ lr ( x+lyS> = —i(x + iy), which yields cos(rx — if) = ye -ra |7(r)| -1 — > 0, as 
y — > oo. The last equation yields the needed quantization for x n . □ 

The following result is due to Widom (see, e.g. [66], Lemma 8.1.9) 

Lemma 17.7. (Widom's lemma). Let F, D be disjoint compact sets in C and C\F- 
connected. Then there is m such that for any Zi, z%, . . . , Z m £ D there is a monic 
polynomial Q m (z) of degree m, such that 



sup 



1 

< - 

- 2 



n 3 =i(*-fj) 

The following result is the mean- values formula for analytic functions of a special 
type. 

17.6. For section [TIT 

Lemma 17.8. Assume that g(\) G B(C + ), g(X) G L 1 (M), g(X) = 5(1) as |A| oo 
and g(iy) — o(y _1 ) as y — * +oo. Then, 

In 11 + ff(A)|d\ = 



Proof. Since Re(l + g) > in C + and g G B(C + ), the function 1 + g is outer from 
iV(C + ). Therefore, 

OO 

s, y f \n\l + g(X)\ „ 



A 2 + y 2 

— oo 

Multiply the last identity by y and take y — » +oo. □ 

Lemma 17.9. ///(A) is suc/i ttat (A + i^fiX) G H 2 (C + ) and f(X) G £ 2 (M) ; 
t^en /(A) G H 2 (C+). 

Proof. Since /(A) = (A + £) 5 (A) with ,g(A) G H 2 (C + ), we have /(A) G JV(C+). 
Then, the statement of the Lemma follows, for example, from the multiplicative 
representation of iV(C + ). □ 

17.7. For section 1121 The following considerations are used in the discussion re- 
garding the case A(r) G L 1 (M + ). We borrow the notations, statements, and proofs 
from [66] ■ Chapter 5. For the reader's convenience, we decided to include this 
material. 

Let X be a Banach space, C- linear bounded operator, and P+- projection (i.e. 
linear bounded operator such that P\ = P+). Notice that P_ = I — P+ is also a 
projection. 

Definition 17.1. The Toeplitz operator is an operator acting in P + (X) by the 
formula T = P + GP + . 

Definition 17.2. A linear bounded operator XL is called upper triangular if P_UP+ = 
and L is lower triangular is P+LP- = 0. 
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Definition 17.3. A linear bounded operator C is a Wiener-Hopf operator if C = £11 
where £,11- invertible, XL, XL -1 - lower triangular, IX, IX -1 - upper triangular. 

Theorem 17.3. (Wiener-Hopf Theorem) . Let G be a Wiener-Hopf operator. Then, 
the corresponding Toeplitz operator T = P+GP+ is invertible and 

T- 1 = {P + U- 1 P + ){P+L- l P + ) 

Assume that Q, R are projections and 

QR = RQ = 0,Q + R = P+ (315) 

Theorem 17.4. (Baxter's Lemma). Let G be a Wiener-Hopf operator so that 
G = £11 = U£. Consider Q, R obeying H315\) . Assume that 

P£Q = RL~ l Q = QUR = QU^R = 

and 

WP-L^KUW < l/2,||mt- 1 P_<C|| < 1/2 
Then, Jq = QGQ is invertible and 

IIT^n < IIL-Hr 1 !! + 2max(||ir 1 ||, ||£- 1 ||)(||p_xl- 1 || + iipu- 1 !!) 

Corollary 17.1. Let G be a Wiener-Hopf operator and {Q„, R n }~ sequence of 
projections obeying 11315]) with Q n x — > x for any x £ P + (X). If they also satisfy 
conditions of Baxter's Lemma, then 

(Q n GQ n r 1 Q n x - Q„7~ l x -> 0, x e P+{X) 

Now, let us apply these mainly algebraic results to the concrete situation. Let 
H(x) £ L 1 (M)- Hermitian function and l + i/(X) > 0, where v{\)~ Fourier transform 
of H. Let X be L 1 (R), Gf = f + H * /, [P+f] (x) = Xn+(x)f(x). The function 
v(X) S W(K) and 1 + v(\) > 0. Therefore, by general result from the Wiener 
algebra theory, we have g(X) = ln(l + v (A)) G W(R) so 



g= J g(x) exp(i\x)dx — J g(x) exp(i\x)dx + J g(x) exp(i\x)dx = g- + g+ 

-co — oo 

Notice that u = exp(g + ) - 1 e W + and I = exp(g_) - 1 £ W_. Therefore, G = £11, 
where £/ = / + 1* f, Uf = f + u*f. Both operators U and £ are invertible and one 
can easily check that tl, II -1 are upper triangular, £, £ _1 - lower triangular (notice 
that at the moment the definition of upper(lower) triangular operator is different 
from what we used in the section on factorization of integral operators). Therefore, 
the Wiener-Hopf theorem is applicable to the operator I + Oi^ = P + GP + . Let 
r(a;) = (J + Jt^)" 1 ^). 

Then, consider the following projections [Q n f] (x) = X[o,n]( x )f( x ) an d [R n f] ( x ) = 
X[n,oo] ( x )fi x )- The result below is what we use in the proof of continuous analog 
of Baxter's Theorem for OPUC. Recall that I + 3~C r is given by ([2"5]) and can be 
regarded as an operator from L 1 [0,r] to L 1 [0,r] due to Young's inequality. 

Corollary 17.2. If n > n , then ||(7 + IK n ) — 1 ||x,i[o 1 n] 1 i 1 [o,n] < @ . Moreover, 
||r re (0,x) - r(x)||x,i[ ,n] -> 0, where r n (0,ar) = {I + JX„) _1 [x[o. n ]( x ) ■ H( x )] ■ 
Proof. Indeed, for n large enough, all conditions of Baxter's lemma are satisfied 
which yields the necessary estimates on the norms. Then, in the Corollary 117. 1[ 
take x = H(t). All conditions of Corollary 1 1 7 . 1 1 are satisfied and we get the second 
statement on convergence. □ 



120 



SERGEY A. DENISOV 



17.8. For section Il3l The next Lemma shows that the spectral measure for Dirac 
operator is uniquely defined. 

Lemma 17.10. The spectral measure dad for Dirac operator D is unique. 
Proof. Indeed, if r(A) is another spectral measure, then Lemma \l 3 . 2 1 yields that 



fW= / f(x)exp(i\x)dx € L 2 {R,dr) 



(316) 



for any f(x) G L 2 [—r,r], r > 0. By taking f(x) — X[o,B]( x ) 1 exp(— x) (R is large 
we have 

dr(X) 



. < oo 
A 2 + 1 

Moreover, from the definition of the spectral measure, 

|/(A)| 2 ^(A)-r(A))=0 



(317) 



Using (I317|) . we can approximate X[a,b) by functions (j316[) in both L 2 (R, dad(X)) 
and L 2 (R,dr(\)). We have 



d(a d (\) - r(A)) = 
for any a and &. That implies dr — da^. 



□ 



The following result is quite elementary. It is used in the proof of existence of 
wave operators for Dirac operator with square summable potential. 

Lemma 17.11. For the unperturbed operator Do, the action of the group e~ %tT>Q 
is given by the formulas 



WD 


' / " 


1 


f(x + t) + f(x-t) 







~2 


-i(f(x-t)- f(x + t)) 



where f{x) G L' 2 (M. + ) is extended to the whole line as an even function. 



itD 


' " 


1 


' -i(f(x + t)- f(x-t)) ' 




. / _ 


~2 


f(x-t) + f(x + t) 



f{x) is extended to K as an odd function. 

Proof. One can use the definition of exp(— UDq) to verify this statement directly. 
Another way to see that is to use the spectral resolution for Dq. 



h 



oo oo 

F(A) = J fi(x) cos(\x)dx + J f2(x)s'm(Xx)dx, (318) 



1 



fi{x) = - / F(X) cos(\x)d\, f 2 (x) = - / F(\)sm(\x)d\ 
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Therefore, 

(e- u ' Do f) 1 (x) = - I e~ itx F(X) cos(Xx)dX, 



OO 

1 



OC 



— oo 

It now suffices to apply the Fourier inversion formula. □ 

f 7.9. For section 1161 The following Lemma proves one simple property of the 
exponential map on the H 1 ' 2 (M.) functions. 

Lemma 17.12. If-y(x) £ i? 1 / 2 (R), then e<^ - 1 6 L 2 (R) and this map is contin- 
uous. 



Proof. We have 

M*) 1 + J2 ||7"(z)l|2 = ||7 * . . . *7||2 



71=1 



( 2 - v \ (2 - p)/{2p) 
IItHUp < C —4 H7llffV*(R), K P < 2 



\2p- 2, 

by Holder's inequality. By Young's inequality we now have 

||7* • • • * 7II2 < ||7llp„,Pn - 2n(2n - I)" 1 

So, ||7™(a;)||2 < Cn n ' 2 ||7||^ 1/2 , R ,. The Stirling formula for factorial yields conver- 
gence of the series and continuity of the exponential map. □ 
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